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IONIZING POTENTIALS OF MULTIATOMIC GASES 
By C. A. MacKay 


ABSTRACT 


Ionizing potentials for helium and fifteen compound inorganic gases were 
measured by accelerating photo-electrons through a gauze into a chamber 
where the positive ions produced were drawn to a fine Pt electrode, made small 
so as to eliminate effects due to radiation. The current-voltage curves ob- 
tained consist of straight lines except at the breaks. Corrections for initial 
velocity, etc., were made by obtaining on the same curve both a break for 
the unknown gas and one for mercury (10.4 volts). The values found are, 
in volts: helium, 24.5; hydrogen, 15.8; nitrogen, 16.3; oxygen, 12.5, 16.1; HCl, 
13.8; HBr, 13.2; HI, 12.8; water, 13.2; NHs, 11.1; Cle, 13.2; Bra, 12.8; Iz, 10.0; 
NO, 9.4; COz2, 14.3; CO, 14.1, 15.6; H,S, 10.4. Comparison is made with results 
of other observers. Theoretical interpretation. An attempt is made to test 
the thermochemical method of estimating ionizing potentials in those cases 
(H2S, NHs, H:O) where data are available. No process of ionization involving 
molecular dissociation seems consistent with the experimental results. This 
and other evidence suggests that in those cases (HCI, HBr, HI) in which the 
thermochemical calculations yield accurate results, the agreement may be 
fortuitous, and that the lowest ionization potentials in each case corresponds 
to ionization without dissociation. 


INTRODUCTION 


N recent years the study of the critical potentials of monatomic gases 

and vapors has been of great importance. The experimental data 
available at present are extensive and accurate, and fruitful correlations 
between electrical and spectroscopic measurements have been found 
which have increased our knowledge of the structure of the atom and of 
the process of radiation and ionization. In the case of multiatomic mole- 
cules, not nearly so many measurements have been made and the ob- 
served values vary from one another over a much greater range than can 
be accounted for by experimental error. In much of the earlier work 
in this field, the distinction between radiating and ionizing potentials 
was not clearly drawn,! which, of course, added to the uncertainty of 
the results. The practically unknown structure of the molecule has made 


1 Davis and Goucher, Phys. Rev. 10, 101 (1917) 
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the interpretation of the experimental results very difficult. The greatest 
progress has been along thermochemical lines,? especially in the case 
of the halogens and some of their salts. The application of the positive 
ray method to molecular ionization*® has also proved valuable, since it 
gives definite experimental knowledge of the various products of ioniza- 
tion, which have to be assumed in any other method of measurement. 
The present investigation was undertaken for two reasons; to measure 
the ionizing potentials of as large a number of multiatomic gases as 
possible, so that more data would be available; and to test the application 
of thermochemical reasoning to additional gases. It was thought that 
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Fig. 1. Diagram of apparatus: A, nickel target; B, platinum gauze; C, fine wire 
electrode; D, electrometer; VV, voltmeters 

interesting results might be found, especially in compounds which are 

non-polar, and that, in connection with other lines of work, our knowledge 

of the process of molecular ionization might be increased. 


METHOD 


The experimental tube and the electrical connections are shown in 
Fig. 1. Light from a quartz mercury arc fell on the nickel target A. The 
electrons emitted by the target were accelerated towards the gauze B by a 
variable difference of potential. lJIonization took place in the region 

? A summary is given in Foote and Mohler, Origin of Spectra, Chapter VIII. 


*H. D. Smyth, Phys. Rev. 14, 409 (1919); Proc. Roy. Soc. A 102, 283 (1922); 104, 
121 (1923); 105, 117 (1924) 
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beyond B and the positive ions were collected on the electrode C by 
means of a small constant difference of potential between B and C of 
either 4 or 6 volts. C was made of fine platinum wire, so that its area was 
reduced to such an extent that the effect of radiation on C was practically 
eliminated. For more complete details and for the method of correcting 
for the initial velocity of the electrons when they leave the surface of the 
target, the reader is referred to a previous paper.‘ 


OXYGEN 


The results have already been given for nine gases‘ and it has been 
shown that in every case except oxygen the agreement with the more re- 
cent determinations has been good. The work on oxygen has been re- 
peated. using the same experimental tube, but preparing the oxygen in 
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Fig. 2. Typical current-voltage curves for oxygen, with nickel target as a source of 
electrons (A); with a hot filament as a source of electrons (B). 


two ways; first, by the electrolysis of sulphuric acid and, later, by heating 
potassium chlorate. Hydrogen was removed by red hot copper oxide and 
a phosphorous pentoxide tube was used to dry the gas. A second series 
of runs was made in which the target A was replaced by a small molyb- 
denum filament, which was electrically heated to a dull red. The latter 
method eliminated the possible presence of ozone which might have been 
produced by the action of the ultraviolet light on the oxygen. 

Fig. 2 shows typical current-voltage curves in which the necessary 
correction in the voltage has been made. In these curves, as in all the 
others, the current scale is purely arbitrary. The reciprocal of the rate of 


‘ Mackay, Phil. Mag. 46, 828 (1923) 
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deflection of the electrometer measured the current and it was plotted on 
a convenient scale. The lower break at 12.5 volts-is present whether the 
source of electrons is the nickel target or the hot filament. In addition 
there is a clearly marked discontinuity at 16.0 volts. Both these values 
must be attributed to molecular oxygen, since the moderately heated 
filament could have produced no appreciable dissociation. 16.1 volts 
has been found by Lockrow® to be the lowest value at which an arc can 
be maintained. 


HYDROGEN SULPHIDE 


The observed value for the ionizing potential was 10.4 volts. Fig. 3 
shows a typical curve B, in which the accelerating voltage was varied 





Fig. 3. Curves for nitric oxide (A); hydrogen sulphide (B); bromine (C); chlorine (D); 
carbon dioxide (E). 


from 0 to 17 volts. The absence of any break at 15.8 volts shows that 
there was very little if any dissociation of the hydrogen sulphide by the 
ultraviolet light, so that the discontinuity at 10.4 volts is the result of 
ionizing the molecule of the compound. The single break is also an 
indication that no large proportion of the molecules is dissociated by 


electron impact. This point will be discussed in more detail in a later 
section. 


CARBON DioxipE, Nitric OXIDE AND CARBON MONOXIDE 


Carbon dioxide and nitric oxide gave curves very similar to that for 
hydrogen sulphide. They are shown in Fig.-3, curves E and A respec- 
tively. Carbon monoxide has two ionizing potentials, both of which are 


§ Lockrow, Abstract of paper presented.at New York Meeting Am. Phys. Sorc. 
February, 1924, Phys. Rev. 23, 553 (1924) 
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clearly marked in the typical curve B in Fig. 4. There is no evidence of 
dissociation. It is interesting to note that carbon monoxide is one of the 
few molecules which occurs with two positive charges in positive ray 
analysis.° The present method gives no direct information regarding the 
products of ionization, but in this case it is at least probable that the two 
values of the ionizing potential correspond to the loss of one and two 
electrons from the molecule. 
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Fig.4 Curves for hydrogen iodide (A) and for carbon monoxide (B) 


CHLORINE, BROMINE AND IODINE 


The vapor pressures of chlorine and bromine are high enough to permit 
direct measurement of their ionizing potentials, while for iodine an indirect 
method was employed. The type of curve obtained for chlorine is shown 
in Fig. 3 curve D. Bromine is very similar, as appears from curve C. 

To obtain the ionizing potential of iodine, hydrogen iodide vapor was 
placed in the experimental tube. Fig. 4 curve A shows the result between 
0 and 17 volts. Four breaks are quite apparent. To show the higher 
breaks more clearly, after the accelerating potential had been increased 
from 0 to 14 volts, the intensity of illumination on the target (Fig. 1) was 
reduced by inserting copper gauze screens and the run was continued-from 
14 volts on, and curve B was obtained. It was observed that when 
readings were made immediately after the light was admitted to the 
tube, the first ionization occurred in the neighbourhood of 13.0 volts. 
Prolonged exposure of the gas gave the type of curve shown, with the 
first ionization at 10.0 volts. Guided by this fact and by the results of 


¢ J. J. Thomson, Positive Rays, pp. 84 and 218 
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previous measurements’ on hydrogen iodide and iodine, the interpretation 
of the four breaks is clear. The lowest is due to the pressure of iodine 
vapor, produced by the dissociation of hydrogen iodine, while that at 
12.8 volts is given by the undissociated molecule. The other two are very 
probably to be attributed to hydrogen, and they agree fairly well with 
the results of other investigations.*® 

For the convenience of the reader the values for all the gases thus far 
measured are given in Table I. Some of the more recent measurements 
by other observers have been added, but no attempt has been made to 


TABLE I 


Ionizing potentials of multiatomic gases 





Gas Observed Values observed References 





values by others 
H: 15.8 16.0 Olmstead® 
15.8 Davis and Goucher! 
16.0 Foote and Mohler, Origin of Spectra, p. 68. 
N: 16.3 16.2 Duffendack, Phys. Rev., 20, 665 (1922). 
16.9 Foote and Mohler, J. Opt. Soc. Am. 4, 49 (1920). 
17.0 Brandt, Zeit. f. — 8, 32 (1921). 
Oz 2.5, 6.1 15.5 io Mohler, Bur. Stds. Sci. Paper, No. 400 
1920 
14.0 Boucher, Phys. Rev. 19, 189 (1922). 
15.5 Smyth, Proc. Roy. Soc. A105, 116 (1923). 
HCl 13.8 13.7 Knipping’ 
Franck, Zeit. f. Phys. 11, 160 (1920). 
HBr 13.2 13.3 Franck, loc. cit. 
HI 12.8 12.7 Franck, loc. cit. 
a + 13.0 Foote and Mohler, Phys. Rev. 17, 394 (1921). 
3 5t. 
He 24.5 24.5 Franck and Knipping, Zeit. f. Phys. 1, 320 (1920) 
Also Franck, loc. cit. 
24.5 Horton and Davies, Roy. Soc. Proc. 97, 23 
(1920). 
- ae Also Davies, Phil. Mag. 45, 797 (1923). 
4 ‘ 
Bre 12.8 11.3 Mohler (unpublished) 
I; 10.0 9.4 Smyth and Compton’? 
10.5 ronson Mohler, Bur. Stds. Sci. Paper No. 400, 
1920). 
NO 9.4 9.3 Hughes and Dixon, Phys. Rev. 10, 495 (1917). 
CO, 14.3 
co . MAA, 8.6 14.3 Foote and Mohler, Phys. Rev. 17, 394 (1921). 
14.5 Stead and Gessling, Phil. Mag. 40, 413 (1920). 
15.1 Found, Phys. Rev. 16, 41 (1920). 
H,S 10.4 (9. 3?) Hughes and Dixon, Phys. Rev. 10, 495 (1917). 





7 Foote and Mohler, Jour. Amer. Chem. Soc. 43, 1832 (1920); 
Knipping, Zeits. f. Phys. 7, 328 (1921); 
Smyth and Compton, Phys. Rev. 16, 501 (1920). 

* Olmstead, Phys. Rev. 20, 613 (1922) 
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give a complete summary of the data. A complete summary of experi- 
mental results and their interpretation is soon to appear in a Bulletin of 
the National Research Council, written by K. T. Compton and F. L. 
Mohler. 


DISCUSSION OF THE RESULTS 


The quantum theory has made possible very many correlations be- 
tween spectroscopic lines and critical potentials for the atoms of many 
elements. When the band spectrum of the multiatomic molecule is con- 
sidered, the problem becomes greatly complicated, since we may have 
energy transfers in the electronic, vibratory, or rotational quantized 
states of the molecule. While the critical potential method’ of attack 
has been employed to interpret radiating potentials, for ionizing poten- 
tials a consideration of the thermochemical relations has led to more 
direct results, at least in the case of certain polar compounds. Calcula- 
tions by means of this method are given by Foote and Mohler” for the 
three halogen acids and the theoretical values agree well with the experi- 
mental data." For other types of compounds the agreement is not so 
satisfactory. Born and Bormann” have predicted 31.8 volts for the com- 
plete ionization of hydrogen sulphide, the final products being assumed to 
be 2(H)*+ and S--. The equation they used may be written, according to 
Foote and Mohler," in the following way: 


(H2S]+Q(H2S) +2D(H:) +2J(H) +5(S) — E(S)>2(H]*+S~ ~ 


where square brackets denote the gaseous phase, Q(H,2S) is the heat of 
formation of H.S, D(H:2) is the heat of dissociation of Hz per 1/2 gram 
molecule, J is the ionizing potential, S is the heat of sublimation per 
gram molecule, E is the electron affinity of one gram atom. From this 
equation we get at once 


J (HS) = Q(H2S) + 2D(H2) + 2J(H) +5(S) — E(S) 


which gives J(H2S) = 31.8 volts if we use Born and Bormann’s values for 
quantities on the right side of the equation. This is far above the lowest 
observed ionizing potential so that it is probable that the first ionization 
does not involve the loss of both electrons. 


* Foote and Mohler, Origin of Spectra, p. 190. 

‘° Foote and Mohler, Origin of Spectra, p. 186. 

1 Mackay, loc. cit.* 

1? Born and Bormann, Zeits. f. Phys. 1, 250 (1920) 
'8’ Foote and Mohler, Origin of Spectra, pp. 180-186. 
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There are several possibilities, for which sufficient data are available 
to get at least approximate results. The process may be simply 


[H.S]+J(H.S)—[H.S}* +e 


where e denotes a free electron, and [H2S]*+ is a molecule of hydrogen 
sulphide with a single positive charge. No thermochemical quantities 
are involved here since we have no dissociation. Or, we may assume that 
we have dissociation but not double ionization and write an equation thus 


(H2S]+ Q(H.S) + S(S)+2D(Hs)+J(H) — E(S)~—Q(H+S~)-H*++ (HS), 


which gives 


J(H2S) =Q(H2S) +S(S)+2D(H2) +J(H) — E(S)- -Q(H+S~) 


In this equation the quantities are the same as in Born and Bormann’s 
with the exception of the last two. E(S)~ means the electron affinity of 
sulphur for one electron and Q(H+S7~) is the heat of formation of the 
hypothetical compound (HS)~ from neutral hydrogen atoms and singly 
charged sulphur atoms. As to E(S)-, it is probable that it has a value 
between E(S)-~ and 0 and nearer the former than the latter; that is, 
nearer 2.0 volts than 0. For an approximate value, let us assume E(S)~ 
is 1.5 volts. Q(H+S~) is also unknown. On Bohr’s theory, S~ would have 
the same number of electrons in the outer shell as the atom of chlorine, 
but probably the tendency to complete the shell would be greaterin thecase 
of chlorine than for S~. This gives us a means to approximate the value 
of Q(H+S~) by making it equal to or less than Q(HCI). The value given 
in the tables of Landolt-Bérnstein for the latter is 5.24 volts, so that 
Q(H+S~) is probably less than 5.2 volts. Given the two quantities E(S)~ 
and Q(H+S7~) we can substitute all the values in the equation given above, 
and obtain 


J(H2S) =.2+2.9+3.6+13.5 — 1.5 —5.2 =13.5 volts. 


ven if the approximations made were considerably in error, we have a 
difference of 3.1 volts between the calculated and the observed values of 
the ionizing potential when we use the most favorable estimate of E(S)~ 
and Q(H+S~—). Probably the discrepancy between calculation and experi- 
ment should be greater. In exactly the same way it may be shown that 
any of the following products of ionization, H.+S+ or H+t+H+5S, or 
H,*++S~ lead to an ionizing potential still higher than 13.5 volts. From 
these considerations one concludes that the lowest observed ionization 
(10.4) consists in the removal of an electron from the hydrogen sulphide 
molecule and that dissociation does not take place. As we have already 
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seen, this is supported by the experimental curve, which shows only a 
single ionization up to 17 volts. 

Both ammonia and water may be treated in the same way as hydrogen 
sulphide. For the former we may write NH;+J(NH;)—NH;++e as the 
simplest possible equation. If, on the other hand, we assume that N, H: 
and H* are formed we have 


NH;+Q(NH;)+D(N)+D(H)+J(H)--N+H2+Ht+e . 


We do not know D(N), but Langmuir’s'* work shows that it is greater 
than D(H); that is, greater than 1.8 volts. This gives us the equation: 


J(NH3)>.4+1.8+1.8+13.5= 17.5 volts, 


where Q(NH;) is taken from Landolt-Bérnstein tables. In the same way 
it may be shown that any assumed process of ionization which does not 
involve the electron affinity of any atom or molecule will give an ionizing 
potential several volts above the experimental result. From the chemical 
properties of hydrogen and nitrogen, it seems unlikely that there is any 
molecule or atom derived from ammonia for which the electron affinity is 
more than a fraction of a volt. We are left with the same conclusion as 
in the case of hydrogen sulphide, that the first ionization is due to the 
expulsion of an electron and that the positive ion is (NHs3)*. 

Turning our attention to water vapor, we find that the data here are 
not sufficiently well known to do much with thermochemical equations. 
The oxygen atom is known to have an electron affinity but its magnitude 
is still undetermined. The same holds true for the OH group which may 
be one of the ions resulting from ionization. We are limited by these 
unknown quantities to consideration of those cases in which the electron 
is supposed to be liberated. A possible equation may be written thus 


H20+0(H20)+J(H2)- Ht +0+e , 
and from this we get 
J( HO) =Q(H,O)+J(He) =2.3+15.8=18.1 volts. 


If the oxygen atom has an electron affinity of 5.0 volts, then the calcu- 
lated result would be reduced by that amount and would be in agreement 
with the observed value, 13.0 volts. Another possible equation is 


H.0+0(H:0)+2D(H:)+J(H)-H++H+0-+e , 
whence 
J( HO) =2.3+3.6+13.5=19.4 volts . 


4 Langmuir, Jour. Amer. Chem. Soc. 34, 860 (1912); 37, 417 (1915) 
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Neither of these equations seems to hold since both give too high an 
ionizing potential and there is no evidence of any ionization of free oxygen 
or hydrogen showing in the curve for water vapor. The more probable 
supposition is that water vapor behaves in the same manner as hydrogen 
sulphide. 

In addition to the negative evidence obtained from considerations 
similar to those just outlined, there is positive evidence as to the nature 
of the first ionization. Smyth® has studied hydrogen, oxygen, and nitrogen 
by means of positive ray analysis. He has found that for all three gases 
the lowest ionization is not accompanied by dissociation, unless as a 
secondary process occurring when the ionized or excited molecule later 
collides with a neutral molecule and either itself dissociates or effects the 
dissociation of the impacted molecule. He states that in view of his 
results, the validity of thermochemical reasoning, when applied to non- 
polar compounds, is extremely doubtful. 

In the case of polar compounds like the halogen acids and salt vapors, 
the case is also doubtful. On the one hand we have the satisfactory agree- 
ment of experiment with thermochemical calculations in the case of the 
halogen acids" and also the recent detection of I~ and Cl- in the products 
of ionization of HgI,and HgCl, in the work of Kondratjeff and Semenoff,’® 
who used Smyth’s positive ray method of analysis of products of ioniza- 
tion. These authors thus detected I~ and Cl- at voltages which agreed, 
within the limits of experimental error, with the values predicted from 
thermochemical considerations. The possible error, however, was several 
volts, and, furthermore, it remains to be discovered whether these nega- 
tive ions were primary products of ionization or secondary products 
arising from subsequent collisions with neutral molecules. This point is 
now being investigated by Smyth and Barton for several halogen com- 
pounds. 

In this connection it may be noted that Duffendack,'’ in intense arcs 
in HCl vapor up to 70 volts, has failed to observe the hydrogen spectrum. 
This is very surprising if the products of ionization are H+ and Cl-, as 
assumed in making thermochemical calculations. 

In conclusion it may fairly be stated that the experimental evidence is 
strongly against dissociation with ionization as a primary process at the 
lowest ionizing potential of non-polar compounds and that, in the case of 
polar compounds, there is still a question as to whether the close agree- 


1 Foote and Mohler, loc. cit.!° 
1° Kondratjeff and Semenoff, Zeit. f. Phys. 22, 1 (1924) 
17 Duffendack, unpublished 
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ment between experiment and thermochemical calculation is not fortui- 
tous, and whether such dissociation as has been observed may not be a 
secondary process. 

The writer takes this opportunity to thank Professor K. T. Compton 
for his most valuable suggestions and unfailing interest throughout the 
course of this work. 

PALMER PHYSICAL LABORATORY, 


PRINCETON UNIVERSITY, 
June 17, 1924. 
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PART I. SOME EXTENSIONS OF THE CORRESPONDENCE 


J. H. VAN VLECK 


THE ABSORPTION OF RADIATION BY MULTIPLY 
PERIODIC ORBITS, AND ITS RELATION TO 
THE CORRESPONDENCE PRINCIPLE AND 

THE RAYLEIGH-JEANS LAW. 


PRINCIPLE 
By J. H. VAN VLECK 


ABSTRACT 


This part deals with the quantum theory aspects of the problem. In the 
absence of external radiation fields the distortion in the shape of the orbit 
is essentially the same in both the classical and quantum theories provided in 
the former we retain only one particular term 71, 72, 73 in the multiple Fourier 
expansion of the force 2e*#/3c* on the electron due to its own radiation. The 
term to be retained is, of course, the combination overtone asymptotically 
connected to the particular quantum transition under consideration. Then 
the changes AJ;, AJz2, AJ; in the momenta J, which fix the orbits and which 
in the stationary states satisfy the relations J, =m,h, are in the ratios of the 
integers 71, Tz, 73 in both the classical and quantum theories, making the 
character of the distortion the same in both even though the speed of the 
alterations may differ. One particular term in the classical radiation force 
is thus competent to bring an orbit from one stationary state to another. 

The correspondence principle is then extended so as to include absorption 
as well as the spontaneous emission ordinarily considered. Commencing 
always with a given orbit it is possible to pair together the upward and down- 
ward transitions in such a way that in each pair the upward and downward 
optical frequencies (determined by the hy relation) are nearly equal for large 
quantum numbers (usually long wave-lengths). That is, if s denotes the 
initial orbit there exist levels r and ¢ such that the ratio (W,—W,)/(W,— Ws) 
OF Yp_/¥s¢ approaches unity when the quantum numbers become large. We shall 
define as the differential absorption the excess of positive absorption due to 
the upward transition s—r, over the negative absorption (induced emission) for 
the corresponding downward transition st. It is proved that for large 
quantum numbers the classical theory value for the ratio of absorption to 
emission approaches asymptotically the quantum theory expression for the 
ratio of the differential absorption to the spontaneous emission. Consequently 
a correspondence principle which makes the numerical values of the emission 
in the two theories agree asymptotically, of necessity achieves a similar 
connection for the absorption. 

The correspondence principle basis for a dispersion formula proposed by 
Kramers, which assumes the dispersion to be due not to the actual orbits but 
to Slater’s ‘virtual’ or “ghost’’ oscillators having the spectroscopic rather 
than orbital frequencies, is then presented. Kramers’ formula has both 
positive and negative terms and the differential dispersion may be defined in a 
manner analogous to the differential absorption. It is shown that the quantum 
differential dispersion approaches asymptotically the dispersion which on the 
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classical theory would come from the actual multiply periodic orbit found in 
the stationary states. This asymptotic connection for the general non-degen- 
erate multiply periodic orbit must be regarded as an important argument for 
Kramers’ formula. 


1. INTRODUCTION 


CCORDING to the first postulate of the Bohr theory of atomic 

structure the electrons can move only in certain particular quantized 
non-radiating orbits or stationary states. In order that the quantum 
conditions may be applicable it is necessary for the motion to be of the 
so-called ‘“‘multiply periodic’ type, which can be represented by multiple 
Fourier series of the form 


x= LYrirsrsX (71,72,73) cos [2r(r 1 + T2We + T3W3)t + Vrizars' *)] (1) 


with similar expansions by y and z. The constants w;, we, ws are the 
intrinsic orbital frequencies, and the summation is to be extended over all 
possible positive and negative integral values of the integers 71, re, 73 
subject to the restriction that 7,w;-+72w2+73w3 be positive. For simplicity 
in notation we have assumed that there are only three w’s, but in general 
the number of such frequencies can be equal to or less than the number of 
degrees of freedom of the atomic system to which the electron under 
consideration belongs. (The slight modifications necessary to extend 
the results of the present paper to systems with more than three frequen- 
cies will be discussed in section 17). We shall also suppose that the 
system is ‘‘non-degenerate”’ so that the number of degrees of freedom is 
equal to rather than greater than the number of frequencies. This means 
that in the case of three w’s, the problem may be thought of as one of a 
single electron moving in an asymmetrical three-dimensional static force 
field. It will also be assumed that the orbits conform to the classical 
mechanics and that all solutions of the differential equations of motion 
can be represented by series of the form (1), thus making the complete 
dynamical system multiply periodic, rather than merely certain particular 
families of orbits. 

The quantum conditions for determining the size of the stationary 
states consist in equating to integral multiples of h a set of orbital con- 
stants J:, J2, J; defined by the relations 


we = OW/aTK(k = 1,2,3); 27 = Jyor + Jews + Jews; (2) 


where W is the total energy and T is the average kinetic energy. We can 
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therefore set J, =n,h (k =1, 2, 3), where m,, m2, m3 are the integers or quan- 
tum numbers characteristic of a stationary state.' 

According to the second postulate of the Bohr theory an electron may 
pass from one allowed orbit to another. The frequency »,, of the quantum 
of light thus radiated or absorbed is determined by the familiar relation 


hv,, = W, — W,, (3) 


where W, denotes the energy of the stationary state 7 (with r=,, s, etc.). 
In the case of emission, 7 is the initial state and s the final state, while for 
absorption the significance of these symbols is just reversed, as r repre- 
sents the higher energy level. 

Considerable information concerning the probabilities of the various 
transitions between different orbits is furnished by Einstein’s derivation? 
of the Planck radiation formula, p(v) =8mvhc-?/(e”/*?—1), for the 
specific energy density of black body radiation. Einstein assumed: 

(I) The number of atoms JN, in the state 7 is given by the statistical 
mechanics formula 

Ne = Ne~W@/k?/ Sre-W/kT (4) 
where N is the total number of atoms, and the sum is to be extended over 
all possible states.* 

(II) The amount of energy emitted in a time A? by transitions from 
the state r to the state s is represented by a formula of the type 


AE, = hires N; = + Bo: p(vrs)| Ad = (5) 


(IIL) The amount of energy absorbed by transitions from the state s 
to the state 7 is given by 


AE.—r -_ hrrs N,B.—, p(vr,)At . (6) 


(IV) There is to be statistical equilibrium (i.e. AE,,,=AE,—,) when 
the energy density has the characteristic black body distribution given by 
the Planck formula. Using Eqs. (3) and (4), one can verify that this 
equilibrium condition will be fulfilled provided the probability coefficients 


1 The quantum conditions have been stated above in what may be termed the corre- 
spondence principle form. Soime readers may be more familiar with the formulation 
used by Sommerfeld and others, which consists in equating certain phase or “quantum” 
integrals to integral multiples of h, but the two methods can readily be shown to yield 
the same results in non-degenerate systems. See, for instance, appendix 7 of Sommer- 
feld’s ““Atombau.” 

? Einstein, Phys. Zeit. 18, 121 (1917) ; 

* The “‘a priori probabilities” p, which ordinarily appear in front of the exponentials 
in Eq. (4) have been omitted, as we are concerned with non-degenerate systems where 
all states have the same a priori probabilities, making the p,’s in numerator and de- 
nominator cancel. 
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A,-45, B,+,, and B,,, which are independent of the temperature and 
energy density, satisfy the relations 


| = Bouse = (c3/Sahvr,*)A re ° (7) 


Einstein’s criterion of statistical equilibrium under black body radia- 
tion thus tells us how emission and absorption vary with the density, 
Eqs. (5) and (6), and how they are related to each other, Eq. (7), but still 
leaves undetermined the magnitude of the coefficient A,,, in (7). To 
evaluate approximately the latter it is customary to resort to the corres- 
pondence principle, as explained below. 

It is well known that according to classical electrodynamics an acceler- 
ated electron radiates energy continuously and that if in particular the 
orbit is of the form (1) then the light thus emitted should be resolved by 
spectroscopes into frequencies which are combination overtones 7,0;+ 
T2wW2+73w3 Of the orbital frequencies w1, we, w3. According to the quantum 
relation (3) there is, however, no such immediate connection between the 
spectroscopic frequency and actual frequencies of motion, but it can 
easily be proved that if the quantum numbers 7, me, 3 of the states r 
and s differ from each other by 71, 72, 73 units respectively, so that 
AJ; =7,h, then in the region of high quantum numbers (usually also 
long wave-lengths) the optical frequency v approaches asymptotically 
the combination overtone 11w,;+7T2w2+73w3. Since the classical and 
quantum mechanisms thus then give nearly the same numerical values 
for optical frequencies even though irreconcilably different in character, 
it is natural to assume that they give in the case of high quantum numbers 
the same numerical results for the relative intensities of different spectral 
lines, at least when there is no radiation field (i.e. p(v) =0). It is, however, 
to be clearly understood that the asymptotic connection of frequencies is 
a necessary mathematical consequence of the quantum conditions, and 
is hence, following Ehrenfest, best termed the correspondence theorem for 
frequencies.‘ On the other hand the existence of an analogous relation for 
the intensities of lines radiated when p(v) =O must be regarded as an addi- 
tional hypothesis, which we shall call the correspondence principle for 
emission. The latter is generally accepted not only because of its in- 
herent reasonableness but also because of its excellent experimental veri- 
fication in the ‘‘selection principle” and in the researches of Kramers and 
others on the intensities of Stark effect components. 


‘To prove the correspondence theorem for frequencies we need simply note that 
hives = W,—W,=AW. Nearly consecutive orbits of large quantum numbers differ but 
little from each other in relative size and we can then without great error replace the 
increment AW of the energy by the differential dW = 2*,(5W/5J,)AJ%. Using Eq. (2) and 
the relation AJ, = 74h, we thus get the result that in the limit »/(1,0:+ 7203+) =1. 
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To formulate analytically the correspondence principle for emission 
we need simply note that from the classical expression (2e?/3c*)#? for the 
rate of radiation from an electron having a vector acceleration , it 
follows that the amount of energy radiated as light of frequency 71; + 
T3wW9-+73ws in the time interval At is® 

= (1624e?/3c?) (Ty04 oa T2We ot T3w3) * [D(r, > T2, Ts)|? At (8) 
where D? = X2+ Y?+ Z*._ If we multiply by N, to take into account the 
radiation from N, electrons, then on comparing this classical expression 
with the quantum Eq. (5) for p(v)=0, and using the approximation 
v=7,01+72W2+73ws, we see that the correspondence principle for emission 
requires that for high quantum numbers the probability coefficient A,-,, 
must have a value close to 

Ary = (162'e?y*,,/3hc*) [Dr(r1 » Ms T3)|?. (9) 
wines D' denotes D evaluated for an orbit of the same size and shape as 
that found in the stationary state r. 

An equation of the precise form (9) can be expected to hold only 
asymptotically. At ordinary wave-lengths it is more probable that 
A,—,, depends not on the orbital frequencies and amplitudes evaluated for 
one particular stationary state but instead involves these quantities 
averaged in some manner or other over the continuous succession of 
‘‘dis-allowed”’ orbits intermediate between the initial and final states. We 
might still have (9) valid in the limit if we assume, for instance, that a 
more exact expression is (10) 


+1 
Aw—».= (16xt'94/3hct)(- =f (rt rat T 33) m{D(r1, T2, ro}tean)” 


Here m is some integer, while \ is an auxiliary parameter such that 
Jy=(nmet+r7Ad)h, where n, and n,+7;, (k=1, 2,3) denote respectively the 
quantum numbers of the states sandr. Either the + or — sign must be 
consistently used throughout. Ordinarily the + sign is taken, and then 
the special cases of m = 0,4 have been studied in some detail by F. C. Hoyt.*® 


2. A CORRESPONDENCE PRINCIPLE FOR ORBITAL DISTORTIONS 


In the review of the corresponderce principle for emission given in 
section 1 it must be remembered that actually a classical electron radiates 
simultaneously all the combination overtones in the multiple Fourier 


‘For simplicity in printing, the arguments 17, 72, 73 of the amplitudes X, Y, Z, D 
will often be omitted; also the subscripts are often dropped from py. 

*F. C. Hoyt, Phil. Mag. 46, 135 (1923); 47, 826 (1924). Eq. (10), of course, gives by 
no means all the possible expressions for A(,_,,). In fact from certain viewpoints the 
best arguments appear to be for a certain type of logarithmic average first introduced 
by Kramers (The Intensities of Spectral Lines, Dan. Acad. Memoires, 1919, p. 330) 
and also studied by Hoyt. Only two out of six formulas studied by Hoyt are included 
in Eq. (10). 
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expansion rather than just one harmonic vibration component 1r,;+ 
T2We+73w3. One can, however, formally avoid this difficulty by intro- 
ducing what we shall term an “abridged” radiation force obtained by 
retaining only the one term 7, 72, 73, in the multiple Fourier expansion 
of the complete radiation force 2e0/3c*. It can then be verified that the 
total radiation is on the average given by Eq. 8.7 We can then imagine the 
emitted light in the classical theory to be monochromatic, and its fre- 
quency will approach asymptotically the optical frequency » in the 
quantum theory, while an asymptotic connection of the average rates of 
spontaneous radiation can be secured by using some such equation as (10). 

A question which naturally arises at this point is whether the abridged 
radiation force with only one term produces just the same distortion in 
the size and shape of the orbit in the classical theory as the electron ex- 
periences in the quantum theory when it passes from one stationary state 
to another. In the classical theory an orbit having originally the same 
size and shape as the initial quantized stationary state will obviously 
after a properly chosen lapse of time be sufficiently damped by the 
abridged radiation force to make the energy the same as that of the final 
stationary state in the quantum theory. However, we cannot imme- 
diately infer that after this interval of time has expired the classical orbit 
will have the same shape as the final quantized orbit unless the system has 
only one degree of freedom, for then only does the energy determine 
uniquely the shape of the orbit. In the case, for instance, of an elliptical 
trajectory modified by a relativity precession to make the system non- 
degenerate (in two dimensions) the classical and quantum orbits might 
have initially the same semi-major axis and eccentricity; but when the 
classical orbit has been sufficiently damped by the abridged radiation 
force to make its energy (and hence approximately the semi-major axis) 
the same as that of another smaller stationary state, we cannot predict 
off-hand that then its eccentricity will become identical with that of the 
latter. 

The writer is not aware of any specific statement in the literature as 
to whether these orbital distortions are the same in the two theories, 
except for a brief allusion to this question in an interesting article by 
H. A. Senftleben, which has just appeared.* However, an examination 

7 The average value of the work — FvAt done in the time At against the abridged 
radiation force F has the value (8) even though the complete expansion of v (but not 
of # in F) is retained, for “‘cross-terms” involving products of different combination 
overtones cancel out on the time average. For greater detail see section 16 of part II. 

8H. A. Senftleben, Zeit. f. Phys. 22, 127 (1924). Following Eq. (50) of his article 


Senftleben states the necessity of having the average time rate of change of the J’s 
approach each other asymptotically in the two theories, and this implies an asymptotic 
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of this point is not difficult. Ina system with 3 degrees of freedom the size 
and shape of the orbit are completely determined if we know the values 
of the expressions J;, J2, Js defined in (2), for the other arbitrary con- 
stants enter only as epoch angles. Now in the quantum theory the J’s 
always change by integral multiples of #, and their alterations in going 
from one state to another are in the ratio 


AJ, :AJ2:AJ3 = 71:72:73. (11) 


- 


On the other hand it is proved in section 16 of part II that if in the 
classical radiation force we keep only the term involving the combination 
overtone 73w1+72 we+73w3, then the J’s are also, in the classical theory, 
altered in the ratio given above for all time intervals. Therefore by 
including only the abridged force in the classical theory, the distortion is 
the same in the two theories. Another way of saying this is that one 
particular term in the classical radiation force damps the orbit in such a 
way that if at the start it coincides with a given quantized stationary 
state it tends to pass through a succession of stationary states of smaller 
energy, whose quantum numbers differ from those of the initial state by 
mr,, MT2, Mrz units respectively, where m is an integer. 

The results quoted in Eq. (11) hold even at ordinary wave-lengths 
(small quantum numbers) for it is shown in section 16 that according to 
the classical theory the J’s always tend to change instantaneously in the 
ratio (11). Consequently the AJ’s are in the ratio of integers even 
though the abridged radiation force be acting for such a long time that 
the alterations in the J’s are of the order of magnitude of their initial 
values, making the radiated energy comparable with the total energy. 
The relations are also valid for systems with more than three J’s (see 
section 17). 

If we form a three-dimensional space for plotting values of J1, J2, Js, 
the deformations of the orbit produced by one particular term in the 
classical radiation force may be represented by a straight line (neglecting 
small periodic fluctuations which cancel out on the average). The equa- 
tions of these lines may be written J, =J,°+r,h, where \ is a parameter. 





connection of the orbital distortions. No detailed proof of this, however, is contained 
in part I of Senftleben’s article (the only part available at time of writing) but an in- 
dication is given of how it might be obtained from an asymptotic connection of the rates 
of radiation of angular momentum and energy in the two theories provided the system 
has only oneelectron. Both the method of proof and interpretation of the result (reached 
independently) in the present paper differ in many respects from the above, especially 
the ability to generalize the results to atoms with more than one electron (section 17), 
and the emphasis on the validity of the relation AJ;:AJ2:AJ = r,:r2:7y at ordinary wave- 
lengths. 
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Consequently the present considerations give a semi-theoretical basis for 
determining the probability coefficient A,, by averaging some function 
of the amplitudes and orbital combination frequencies along a straight line 
in the three dimensional space connecting the points corresponding to the 
initial and final states. This procedure, illustrated in Eq. (10), was 
adopted by Kramers and by Hoyt, but it is not apparent from the articles 
of these investigators whether they realized they were actually aver- 
aging along a path which would be traversed by the orbit under the 
influence of the abridged radiation forces. The various methods which 
have been suggested for determining A,—,, by averaging over the straight 
line path mentioned above are such as to make the mean free time in the 
initial state different from the actual time of transit between the two 
states under the influence of the abridged radiation force, except, of 
course, for an asymptotic connection of the two times for high quantum 
numbers. Neglecting higher powers of 1/c?, the quantum mean free time 
and the classical time of transit can be shown to be identical at ordinary 
wave-lengths or small quantum numbers provided in Eq. (10) we take 
the — sign in both exponents and take the integer m= —3.° This result 
seems worth noting, but it is questionable how much significance should 
be attached to this method of procedure for there is no obvious reason why 
the two times should be exactly equal rather than asymptotically con- 
nected.!° 


3. INSUFFICIENCY OF A CORRESPONDENCE PRINCIPLE RELATING ONLY 
TO EMISSION 


The correspondence principle for emission reviewed in section 1, cannot 
in a certain sense be regarded as entirely adequate because it establishes 
(or rather postulates) an asymptotic numerical connection of the classical 


* To get the result quoted we need simply integrate Eq. (8) between \=0 and A=1, 
noting that dE =(1:0:-+ r2w2+73w3)hddA by Eq. (2). This gives the classical transit 
time. The quantum mean free time is on the other hand simply the reciprocal of A;—.. 

10 [t is interesting to note that the classical time of transit from a 33 to a 2; orbit in 
hydrogen works out as 1.04X10-* sec. Wien’s canal ray experiments indicate a mean 
free time for Ha of 1.85 X10-8 sec. (Ann. der Phys. 73, 485, 1924). The latter, Wien 
shows, is almost exactly the reciprocal (1.87 X10-® sec.) of the logarithmic decre- 
ment for a linear oscillator of frequency v (not an actual orbit), but this very 
close agreement is probably only a coincidence. On the other hand the value 
1.04X10-* sec. must be corrected to allow for the transitions 3; 2;, 3,723, which also 
contribute to Ha and have longer mean free times than 3;-*>22,—the 3, state is, in fact, 
almost metastable—thus making the effective mean free time much larger. The agree- 
ment is as good as can be expected, for the correct formula for A,—, probably does not 
make the quantum mean free time exactly equal to the classical time of transit, but 
instead involves a different kind of average of amplitudes and frequencies than that 
obtained by putting m= —3 in Eq. (10). 
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and quantum theory for intensity only in the special case that there is no 
radiation field (i.e. p(v)=0). The term present in Eq. (5) when p(v) =0 
we shall call the spontaneous emission (Einstein’s ‘‘Ausstrahlung’’), 
while the second or remaining term proportional to the energy density 
we shall call the induced emission, although it is sometimes called the 
‘‘negative absorption” in distinction from the true or positive absorption 
given by Eq. (6). The correspondence principle for emission correlates 
quantum theory spontaneous emission with classical theory emission due 
to the radiation force 2e’#/3c*. On the other hand an electron can accord- 
ing to the classical theory absorb energy from a radiation field, and there 
must be some kind of an asymptotic connection between this absorption 
in the classical electrodynamics and the induced emission and absorption 
in the quantum theory. This will be discussed in sections 4 and 5. There 
can be no question of a correlation of the classical and positive quantum 
absorptions alone, as this would leave the induced emission unexplained. 
The existence of the induced emission term in the quantum theory may at 
first sight appear strange, but it is well known that this is qualitatively 
explained in that with the proper phase relations a classical electric wave 
may receive energy from an atomic system although on the average 
(i.e. integrating over all possible phase relations) it contributes more than 
it receives in exchange. It is therefore the excess of positive absorption 
over the induced emission which one must expect to find asymptotically 
connected to the net absorption in the classical theory. 


4. CORRESPONDENCE PRINCIPLE FOR ABSORPTION FOR A LINEAR 
OsCILLATOR 


Before seeking to develop a correspondence principle for absorption for 
the general case of an arbitrary multiply periodic orbit, we shall first 
for simplicity and clarity confine our attention to a one-dimensional 
linear oscillator. Here the multiple Fourier expansion reduces to 


x = Decos (2mwt + yy), (y =z = 0). 


As there is only one degree of freedom, there is just one quantum number 
n, and this can by the correspondence principle only change by one unit, 
as there are no harmonics of the fundamental frequency w in the Fourier 
expansion. Now the energy W of a linear oscillator of amplitude D, 
mass m, and frequency w is in general 27*w*mD?, Furthermore it is well 
known that for a linear oscillator the quantum conditions require that 
the energy W, of a state of quantum number m have the value nhw. 
Therefore the amplitude D, of an orbit of quantum number 2 is given by 


D,” = nh / (23*vm). (12) 
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In writing (12) we have utilized the familiar fact that for a linear oscilla- 
tor y=w, as here the spectroscopic and orbital frequencies are identical." 
We can now apply the correspondence principle for emission, as embodied 
in Eq. (9), to determine the approximate value of the probability coeffi- 
cient A,,,. Using the fact that in the present notation D"(r;,r2,73) is 
nothing but D,, while v,, is simply v, and taking r as an m quantum state, 
s as an (n—1) quantum state, we thus get A,_,, =82°e*v’n/(3c'm). From 
(7) it then follows that 


Brosinet) = Ben-y-—n = ne? / (3hmr). (12a) 


Hence by Eq. (6) we see that in a time At the energy which is removed 
from a radiation field of density p(v) by the positive absorption of quanta 
by JN oscillators all in (7—1) quantum states is 


AE (n-:)-sn = }nre?m= Nop(v)At . (13) 


On the other hand it is well known that according to the classical 
theory the average rate at which a linear oscillator absorbs energy in a 
field of radiation is independent of the amplitude and is given’ by 
tre’m—'p(v). Consequently in the interval At, N oscillators should 
absorb the energy 

AEgss = jre*m—! Np(v)Al . (14) 


Because of the presence of the factor m, the quantum theory expression 
(13) for the positive absorption differs increasingly from (14) as the 
quantum number becomes larger. This discrepancy is not surprising, for 
we have not taken into account the fact that oscillators in the state 
(w—1) may in the presence of a radiation field be induced to emit energy 
and pass to the state (7 —2) at a more rapid rate than when p(v) =0 and 
there are only spontaneous transitions."* Each of these excess or induced 
transitions may be thought of as returning to the ether or light wave the 
energy /yv,,, a sort of regenerative effect. Hence we can take as the net 
absorption of energy in the time At by JN oscillators in the state (m—1) 
the expression 


AF = [Bn—1)—+n = Bn-) 5 (n-2) ] hy Nop(v) ai. 


' The result y= is obvious from an inspection of Eq. (3), as in the present case 
W, =nhw, An= +1. 

12 Cf, Planck, Warmestrahlung, 4th Ed., Eqs. (260) and (159) 

18 If the electrons are in the orbit of lowest quantum number there can be only 
positive absorption because there are no still lower energy levels to which induced 
emission may take place. This case has been considered by Ladenburg (Zeit. f. Phys. 
4, 451, 1921); also by Ladenburg and Reiche (Naturwissenschaften, 27, 584, 1923). 
These writers note the necessity of having the factor m in comparing the classical and 
positive quantum absorption. 
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This way of measuring absorption we shall term the differential rate of 
absorption in contrast to the true positive rate given by Eq. (13). Now by 
(12a) we have,on lowering the integer” by one unit, the result Bin_1)—(n-2) = 
(n —1)2e?/(3mhv), and hence AF=AE,,. _ We thus see that in the limiting 
case of large quantum numbers, where Eq. (12a) is valid, the classical 
value (14) for the rate of absorption of energy is nothing but the differen- 
tial rate of absorption in the quantum theory. This connection of the 
classical and quantum differential absorption we shall term the corres- 
pondence principle for absorption; it is a purely mathematical conse- 
quence of the correspondence principle for emission, which was used in 
deriving (12a). Another way of stating the results is that because of 
Eqs. (9) and (12a) the ratio (Bin—1)n — Ben—1)—(n-2)) p(v)/A (n—1-+(n-2) of the 
differential absorption to the spontaneous emission for any given orbit 
has for large m very nearly the classical value c*p(v)/(162°mD,?v") for 
the ratio of absorption to emission." 


5. GENERALIZATION OF THE CORRESPONDENCE PRINCIPLE 
FOR ABSORPTION TO AN ARBITRARY NON-DEGENERATE 
MULTIPLY PERIopIC SYSTEM 


We must now seek to examine whether an analogous correspondence 
principle for absorption holds for an arbitrary non-degenerate multiply 


periodic system, whose orbits are given by Eq. (1) and which represent 
essentially the most general type of motion amenable to present quantum 
theory methods. Here we can no longer enjoy such simplifications as the 
existence of a single quantum number and the results An=+1, v=w 
characteristic of a linear oscillator. 

Let us consider a state s of quantum numbers m, m2, m3." Then when a 
quantum of energy is absorbed by an orbit of this type the electron will 
pass to some other state of quantum numbers 2,;+7, me+72, m3+73."° 
Starting with the state s (not r) there is no emissive transition to a state 
t such that the frequency »,; of the light thus emitted is just equal to the 
frequency v,, of that absorbed on passage to the state r. (This amounts 


Reference should be made while on the present subject to the discussion in section 
158 of Planck’s “‘Wiarmestrahlung” (4th Ed.) ,which establishes an asymptotic connec- 
tion between the classica] and quantum formulas for statistical equilibrium. A corre- 
spondence principle for absorption could doubtless be derived from such considerations, 
but no explicit mention is made of the asymptotic connection of the classical absorption 
and the differential absorption for a single orbit (where thermodynamic equilibrium 
need not be assumed) which is the primary concern of the present paper. 

 The’generalization to cases of more than three quantum numbers is given in section 
17. 

'® Not all the integers 71, 72, rs; need be positive but the state r must have a greater 
energy than s. 





SOME EXTENSIONS OF THE CORRESPONDENCE PRINCIPLE 341 


virtually to saying that in general energy levels are not evenly spaced, 
except, of course, for such an ideal case as the linear oscillator). However, 
if the state ¢ has quantum numbers 2,—71, m2—72, m3—73, then for large 
values of m,, m2, ms the ratio v,,/v,, approaches unity, for both these 
spectroscopic frequencies approach asymptotically the combination 
overtone 73#1+72W2+73w3 of the orbital frequencies. Hence it is natural 
to define as the differential rate of absorption the excess of light absorbed 
by transitions from the state s to the state r over the energy returned 
to the ether by induced emissive transitions from the state s to the state 
t, where, as mentioned above, the differences between the quantum 
numbers for the states s and ¢ are the same as the corresponding differen- 
ces for the states r and s. 

With this definition it follows that in the time At the differential 
absorption of light energy of frequency approximately 7,#,+Tewe+r3ws 
by N atoms in the state s amounts to 

AF = [hypp(vre) Ber — hyetp(ver) Bei] NAL. (15) 
If we admit the validity of the correspondence principle for emission 
then Eqs. (7) and (9) show that for large quantum numbers this expres- 
sion becomes approximately 

AF = 3x%e?h- [prep (Vrs) D,? — veep(ver) D,’| NAt 

where for brevity we have written D, for D’ (7, 72, 73), etc. Also for large 
quantum numbers the discrete succession of quantum orbits becomes so 
nearly consecutive that differences may be replaced by differentials, just 
as is done in the derivation of the correspondence theorem for frequen- 
cies.4 Consequently 


0 
AF = ane? h-! { 0)D.4, — Vet) +6 (*) D,(vrs = Vet) + Vetp(Ve4)6D" \ NAl 
v 


where 6 denotes the difference between an expression evaluated for the 
states s and ¢. It is, of course, to be understood that p(v) and dp/dv are 
to be evaluated with »v approximately equal to »,; (or v,.). Now since we 
are dealing with the large quantum number limit, the correspondence 
theorem for frequencies tells us that »,, may without sensible error be 
replaced by the combination overtone 171w:+72w2+73w3 of the orbital 
frequencies, and similarly v,,—¥v,; is to the desired approximation equal 
to 5(73w1+72w2+73w3). Now since the system is non-degenerate, the 
amplitude D as well as the w’s and the energy W, is a function of the 
constants J;, Jz, Js which according to the quantum conditions are 
equated to integral multiples of h, and consequently 


..... 421, , 4 Oats i( : + : + : ) MU, =rnh 
= nae = hf rs—+ r2-—- + 173— Joi, asdJp= Kh. 
a a — < —”CC 
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Similar relations hold for dw2, 5w;, and 6D?, The expression for the differ- 
ential absorption therefore — 


te] 
AF = $n°e? [ plo (n ", +9 oe + na) Gr + 


OS 2 " 
Op re) 
2 Ge(r poe +r a nae) v-| NAt (16) 
Ov "as, "iJ 


where we have written w, for 711 +72w2+73w3 and G for w,[D(r1, 72, 73)]?. 

It is shown in part II of the present paper that Eq. (16) also gives the 
amount of energy which according to purely classical mechanics is ab- 
sorbed in a time At by N systems, each with a multiply periodic orbit 
similar to the state s, when exposed to a radiation field of vanishing 
intensity except for frequencies in the vicinity of 713w:-+7T2w2e+73w3. That 
is to say Eq. (16) also gives the part of the classical absorption due to 
resonance of the impressed waves with the combination overtone 71, 72, 73. 
The total classical absorption is, of course, the sum of resonance effects 
for all possible overtones. 

We may therefore conclude that the correspondence principle for 
absorption holds even for an arbitrary multiply periodic orbit; i. e. 
assuming the validity of the correspondence principle for emission, the 
differential quantum absorption by a particular orbit approaches at high 
quantum numbers (usually long wave-lengths) the classical absorption of 
light of asymptotically corresponding frequency." 


" It is clearly to be understood that although in deriving (16) we have apparently 
utilized Eq. (9), the proof of (16) can be readily shown to be equally valid if (9) holds 
asymptotically, so that at ordinary wave-lengths A,—, might, for instance, be deter- 
mined by averaging the frequencies and amplitudes in the manner given in Eq. (10) 
In this connection it is interesting to note that it is not difficult to prove that for the 
special case of a linear oscillator the differential quantum absorption and the classical 
absorption are exactly equal to each other even for low quantum numbers provided 
A,—, is determined by taking the + signs in (10), with the integer m arbitrary. On 
the other hand it was found in section 2 that to make the mean free time of spontaneous 
emission in the quantum theory just equal the classical time of transit from one state 
to another in the absence of external radiation (i.e.p(v) =0), it is necessary to take the 
— signs in (10) and m= —3. 

Another way of formulating the correspondence principle for absorption is as follows. 
The quantity hvyA,—, plays in the quantum theory the réle of a coefficient of emission, 
and may be denoted by a. Using arguments similar to those employed in the derivation 
of (16), it is easy to show that the differential absorption has for large quantum numbers 
very approximately the value 


0 (] 
AF = (c3/8r) [ peor (ns arTs + 7, + te ~) (aws*) + 
1 3 


= : +? - + -.) | a 
wr — { r;— — — }wr < 
ary; 1 "Ts "3 Js . 
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In the present paper we have proved the correspondence principle for 
absorption after assuming the validity of the ordinary correspondence 
principle for emission. The procedure could equally well have been re- 
versed. That is to say, assuming an asymptotic connection of the differ- 
ential absorption in the quantum theory with the corresponding absorp- 
tion in the classical mechanism we could show that the spontaneous 
emission in the quantum theory must approach equality with the classical 
formula (8) for the spontaneous radiation of energy. 

The classical formula (8) for spontaneous emission involves rather 
more electrodynamics than the analogous classical formula for absorp- 
tion, given in Eq. (16), for the derivation of the former is based on 
retarded potentials and the validity of the field equations, while the latter 
requires simply that there be a certain incoherence in the radiation field 
and that the change in the energy be equal to the scalar product of the 
impressed force and the velocity of the particle, a very broad mechanical 
principle. As pointed out in the preceding paragraph, the asymptotic 
values of the probability coefficients A,,,, B,,, etc., can be derived from 
the absorption instead of the emission viewpoint. Because of this alterna- 
tive the correspondence principle for evaluating these coefficients is in a 
certain sense less electrodynamical and more purely mechanical in nature 
than it has sometimes been considered to be. From this standpoint the 
connection between the quantum theory and classical electrodynamics 
comes through the assumption of Wien’s law in obtaining Eq. (7), for in 
Einstein’s derivation of the Planck radiation formula considerations of 
statistical equilibrium without resorting to Wien’s law do not suffice to 
determine B,,,/A,—,. If we assumed both the correspondence principle 
for emission and that for absorption it would be possible to dispense with 
the need of assuming (rather than proving) Wien’s law. Virtually this 
procedure was adopted by Planck, who studied the special cases of the 
linear oscillator and rotating dipole.'® As pointed out to the writer by 
Prof. Ehrenfest, the content of the present paper may be regarded as 
placing such a procedure on a more general basis."® 





This is, however, nothing but the classical theory formula for the coefficient of ab- 
sorption in terms of the coefficient of emission as can be seen from (8) and (16). This 
way of stating the correspondence principle for absorption does not utilize (9) even 
asymptotically. 

18 Wiarmestrahlung, 4th Ed., section 158. Cf. footnote ™. 


19 The writer has just learned (Sept. 1924) that in unpublished computations made 
at Copenhagen, J. C. Slater has independently derived an absorption formula similar 
to Eq. (16) and has also noted the asymtotic connection of the classical and quantum 
absorption discussed in sections 4 and 5. 
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THE GENERAL CORRESPONDENCE PRINCIPLE BASIS FOR KRAMERS’ 
DISPERSION FORMULA 


In a recent note”® H. A. Kramers has proposed a formula for dispersion 
which is a modification of an equation previously developed by Laden- 
burg and Reiche,”! and which must be regarded as a distinct advance in 
the problem of reconciling dispersion with quantum phenomena. The 
physical principle underlying Kramer’s formula is that of “‘virtual’’ or 
“ghost”? resonators, first suggested by Slater” and elaborated by Bohr, 
Kramers, and Slater.» According to this viewpoint the dispersion is not 
to be calculated by considering the actual orbit (stationary state) as 
reacting classically to the impressed waves. Instead, the stationary states 
appear to be unaffected except for occasional quantum leaps, but the 
dispersion is to be computed as due to a set of hypothetical linear oscilla- 
tors whose frequencies are the spectroscopic ones rather than those of the 
orbits. The introduction of these virtual resonators is, to be sure, in some 
ways very artificial, but is nevertheless apparently the most satisfactory 
way of combining the elements of truth in both the classical and quantum 
theories. In particular this avoids the otherwise almost insuperable 
difficulty that it is the spectroscopic rather than the orbital frequencies, 
i.e. vy rather than 74w1+72w2+73w3, which figure in dispersion.” 

Kramer’s formula for the polarization due to N, electrons in the state 
sis” 


8 


r—s Let ae ~ 
P, = 324 [> new 2 = | N.€ . (17) 


Vrs” Cs waell4 as See vy”) 


*” H. A. Kramers Nature, May 10, 1924, p. 673; Aug. 30, 1924, p. 310. 

*! Ladenburg and Reiche, Naturwissenschaften 27, 584 (July 6, 1923). Ladenburg 
and Reiche’s formula differs from that of Kramers, given in (17), principally in that 
the negative terms are not included. 

*2 Slater, Nature, Mar. 1, 1924, p. 307 

** Bohr, Kramers and Slater, Phil. Mag. 47, 785 (1924) 

“ Reference should, however, be made to the alternative explanation proposed by 
Darwin (Proc. Nat. Acad. 9, 25, 1923). 

*® This equation differs by a factor 1/3 from the formula in Kramer's note to Nature, 
as (17) is written for the case where all orientations of the atom are equally probable, 
while Kramers’ original equation is for the case that the atoms are always so oriented 
that their free vibrations are parallel to the impressed electric field. The assumption 
that all orientations are equally probable may appear at first sight contrary to the 
assumption that the system is non-degenerate, as all three intrinsic frequencies (the 
criterion for non-degeneracy) can appear only when the atomic force field does not have 
a spacial symmetry, and with such a dissymmetry the quantum conditions allow only 
certain particular orientations. However, we may imagine the asymmetrical internal 
atomic force fields to have all possible orientations in the different atoms. This artifice, 
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where @ is the electric intensity and r is the frequency of the impressed 
wave. The summation with respect to r is to be taken over all states of 
energy content higher than s, while that for ¢ over all orbits of lower 
energy content than s. Each term in this equation represents the disper- 
sion due to a “‘virtual oscillator” and corresponds to a particular quantum 
transition starting from the state s. Each oscillator in the ¢ summation 
gives a negative dispersion, corresponding to Einstein’s negative absorp- 
tion, mentioned in section 3. 

Kramers states that (17) merges asymptotically for large quantum 
numbers into the classical formula for the polarization. It is our purpose 
to show that this is true not just when the quantized system is a linear 
oscillator, but also when it is the most general type of non-degenerate 
multiply periodic orbit.% Just as in section 5, the upward transitions 
(s-+r) and the downward ones (st) may be mated together in such a way 
that in each pair v,,/v., converges to unity for large quantum numbers. 
This is accomplished, just as before, by making the differences in quan- 
tum numbers between the states r and s the same as those between the 
states s and #, and we shall denote these differences by 7;, r2, 73. Each pair 
contributes a positive and a negative term in (17) and their net effect is 
a “differential dispersion’? very analogous to the differential absorption 


mentioned in sections 4 and 5. In precisely the same way that (16) is 





introduced for simplicity, enables one to have a random spacial distribution of orbits 
without having degenerate systems. 

In this connection it is interesting to note that if the axes of the asymmetrical atomic 
force fields had the same orientations in all atoms, as is the case when the dissymetry is 
due to a constant external magnetic or electric field, then we should expect the virtual 
resonators to have only certain particular orientations, as only a few quantized spacial 
positions can be assumed by the atoms. By analogy with the classical theory results 
mentioned after Eq. (41) in section 15 of part II, there might then be possible a polariza- 
tion in the y and z directions even though the electric intensity coincides with the x axis. 
The polarization and electric vectors would then have different directions, as in a crystal. 
Therefore the application of a weak constant external field would cause an abrupt dis- 
symmetry in the index of refraction just as it creates an abrupt outstanding polarization 
of the spontaneously radiated light (cf. N. Bohr, the Fundamental Postulates of the 
Quantum Theory, Supplement of Proc. Cambr. Phil. Soc., p. 27). These sudden dis- 
continuities are not found in the classical theory because the spacial distributions are 
continuous rather than quantized. 

*In his first note Kramers did not give his proof of the desired asymptotic 
connection. Inasmuch as the formula for the dispersion by an arbitrary multiply periodic 
orbit, to be developed in section 15 of part II, did not appear to have previously been 
given and was not mentioned in Kramers’ note, the writer supposed that Kramers’ prior 
demonstration was somewhat less general than that of the present paper. Since writing 
this paper, the author has learned, however, that Kramers’ unpublished computations 
are of the same scope and generality as those in § 15. Cf. Kramers’ second note”. 
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derived from (15) we find that for very large quantum numbers the 
contribution of one —— wee to the polarization becomes 


1 Ce) 
er ¥ [= ~ 9) id = * s+ rc); 
awe a a er 
J. a” a” Cae 
where G, and w, are defined as in Eq. (16). 

It is proved in section 15 that (18) also gives the polarization which 
according to the classical theory would be produced by resonance be- 
tween the impressed frequency and the combination overtone w, in a 
multiply periodic orbit of the same size and shape as the state s. The 
differential dispersion in Kramers’ theory therefore approaches asymptot- 
ically the classical dispersion for the actual orbit. The asymptotic con- 
nection is thus very similar to that previously encountered in the corre- 
spondence principle for absorption. 

It is particularly interesting to note that although both the positive and 
negative terms in the differential dispersion taken separately represent a 
type of dispersion characteristic of a linear oscillator (except that oscilla- 
tors with negatiye dispersion correspond to no ordinary physical reality, 
as they would have to possess a negative mass) the difference of the 
two terms approaches asymptotically a more complicated type of dis- 
persion appropriate to the general multiply periodic orbit. For Eq. (18) 
contains a term in (w,?—v*)~? as well as the familiar term in (w,?—v?)7! 
characteristic of the ideal resonator. The fact that the dispersion for the 
simple virtual oscillators thus merges asymptotically into the more 
complex classical dispersion for the actual orbits must be regarded as an 
important argument for the virtual resonator viewpoint, as such a con- 
nection is not an obvious outcome of the theory. Especially does this 
asymptotic connection increase our faith in the particular form of virtual 
oscillator theory embodied in Kramers’ dispersion formula, and makes 
it easier to accept the rather artificial negative dispersion terms thus 
involved.”’ 








7 It is to be noted that if A,—, is determined in accordance with Eq. (10), the + 
signs being taken, then the differential dispersion in the Kramers’ formula is for the 
special] case of a quantized linear oscillator easily proved identical even at low quantum 
numbers with the classical dispersion. It was pointed out in footnote ” that this kind 
of formula also makes the differential and classical absorption identical for the ideal 
oscillator without passing to the high quantum number limit. These facts taken together 
furnish a limited amount of evidence for the use of the average of the form (10) (with 
+ signs) rather than the alternative logarithmic average which Hoyt has shown explains 
observed intensities equally well in the instances studied by him.* 
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THE ABSORPTION OF RADIATION BY MULTIPLY 
PERIODIC ORBITS, AND ITS RELATION TO 
THE CORRESPONDENCE PRINCIPLE AND 

THE RAYLEIGH-JEANS LAW 


PART II. CALCULATION OF ABSORPTION BY MULTIPLY 
PERIODIC ORBITS 


By J. H. VAN VLECK 
ABSTRACT 


This part is largely dynamical and presents the mathematical analysis on 
which part I is largely based. It gives the calculation of the absorption by a 
multiply periodic orbit in a field of radiation and the corresponding perturba- 
tions in the orbits, also of the polarization produced by incident waves, causing 
dispersion, and the distortion of the orbit by its own radiation. Most of these 
calculations have previously been given only for special cases. 


7. INTRODUCTION. GENERAL DYNAMICAL THEORY OF MULTIPLY 
PERIODIC ORBITS 


It is the purpose of part II primarily to calculate according to the 
classical mechanics the average absorption and certain other properties of 
a multiply periodic orbit in a radiation field. Besides being of some 
interest taken by itself, this computation is necessary to furnish a 
mathematical basis for much of the material in parts I and III. Although 
multiply periodic orbits have been extensively studied by many investi- 
gators, the calculation of the absorption appears to have hitherto’’ 
been published for only a few special cases, such as the linear oscillator 
and one-dimensional rotating dipole whose rates of absorption have been 
examined by Planck and Fokker respectively.”* 

Let us consider an electron whose Cartesian coordinates are x, y, z, and 
which in the absence of a radiation field is subject to forces which can be 
derived from a static potential V(x, y, z). We thus assume that the 
system has only three degrees of freedom and hence but three intrinsic 
frequencies. These restrictions are made only in order to avoid excessively 
complex notation. (The appropriate modifications for the more general 
case of several electrons and hence more than 3 degrees of freedom are 
given in section 17.) Instead of using the ordinary three second order 


**M. Planck, Wairmestrahlung, Part IV, Ch. 1, 4th Ed. (Part IV, Ch. 2 of Engl 
trans. by Masius). A. Fokker, Ann. der Phys. 43, 812 (1914). 
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Newtonian equations (mx = —dV/dx, etc.), the motion can equally well 
be specified by the six first order Hamiltonian canonical equations 
dp: oH dx oH : 
= a eee, —— Ben. CC. (19) 
dt Ox dt Op. 
where the Hamiltonian function H is the total energy 7+ V, so that 
H = (1/2m) (p22 + py? + pe?) + Vix, y,2). (20) 


It is convenient to make a change of variables from x, y, z, pz, py, p, to 
new coordinates w;, We, Ws, J1, Je, Js. In a multiply periodic system it 
is possible to find a transformation possessing the three following proper- 
ties. 

(1) The six original coordinates and momenta, x, p., etc. are periodic 
in each of the w’s with modulus 27, and so can be expanded in series of 
the form?® 

% = Lr A(71,72,73) exp [27i(riwi + rewe + T3Ws) | (21) 


with similar equations for y, 2, pz, py, p:, except that different coefficients 
appear in front of the exponentials. The summation is to be extended 
over all positive and negative integral values of 71, 72, 73. By Fourier’s 
theorem 


2 2r le 
A(r 9T2 »T3) =(2n)-* f s f x exp[— 2ri(r1w1+72w2+73wW3) |X 
0 0 0 


dw,dwedw; (22) 
and consequently A(—7;,—72,—73) is the conjugate of A(r1,72,73), 
making the resulting series real despite the imaginary exponentials. The 
use of this complex variable notation makes the equations very much 
more compact and easier to handle than if sines and cosines are employed. 

(II) The expression p,dx+ p,dy+p,dz—J,dw; — J2dw2—J3dw; is the 
exact differential dF of some function F(x, y, 2, wi, We, W3) when the p’s 
and J’s are expressed in terms of x, y, 2, Wi, We, Ws by means of the trans- 
formation formulas. Such a transformation is often termed a contact 
transformation. We shall furthermore suppose the time average of 
the function F is zero, so that the expression ~,¢+),)+ ,z (which is 
twice the kinetic energy) is on the average equal to Jiw,+J2w2+Jsws. 

(III) The energy or Hamiltonian function H given in (20) is expressible 
in terms of J;, Jo, J; alone, (i.e., the w’s drop out when H is expressed 
in terms of the J’s and w’s by means of the transformation formulas (21)). 

In virtue of the fundamental theorem that a contact transformation 
preserves the Hamiltonian form,*® the differential equations for deter- 
mining how the w’s and J’s vary with the time may be written 


** For simplicity in printing we have denoted the exponential function e” by exp|id]. 
*° For proof of this theorem, see, for instance, appendix 4 of Sommerfeld’s Atombau. 
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dJ oH dw, 0H 
dt dwe dt ai 
where the subscript k may have the values 1, 2, or 3. This range of values 
for k is to be understood throughout, while the letter i will be reserved 
strictly for »/—1. Now by (III) we have 0H/dw,=0, and hence the J’s 
are constants independent of the time. Therefore the time derivatives of 
w’s, which by (23) are functions only of the J’s (since 7 involves only the 
J’s), are also constant, and consequently 
Wk = wet + €x (24) 
where the €, are arbitrary constants or phase angles, while the expres- 
SIONS W,W2,w3; are functions only of the J’s and are the intrinsic fre- 
quencies of the system. We shall suppose that the system is non-degen- 
erate, so that there is no commensurability relation of the form Jw,;+ 
mwe+nw3;=0, where 1, m, and n are integers.* 
By introducing the values of the w’s given in (24) and using the conju- 
gate property mentioned in connection with Eq. (22), it is seen that 
Eq. (21) may be replaced by a series of cosine terms of the form (1) 


(23) 


provided 
|X (71,72,73)]? = 4A(r1,72,73) A(—71, —T2, —T:s) , (25) 


Yann’ = 2r(r1€1 + T2€2 + T3€3) + arc A(11,72,73). 


According to the quantum conditions, the constants J;, J2, Js should be 
equated to integral multiples of 4, as it follows from (23), (24), and the 
latter part of postulate (II), that the J’s have the properties demanded 
in Eq. (2). However, this need not concern us immediately, as the 
present computation is based on purely classical mechanics without 
quantum modifications. 


8. DERIVATION OF THE EQUATIONS FOR THE PERTURBATIONS CAUSED BY 
ADDING A RADIATION FIELD 


In the preceding paragraphs we have reviewed briefly some of the 
general properties of multiply periodic dynamical systems which it will 
be necessary to use in the forthcoming analysis.** Let us now assume 


*t When there exists a commensurability relation of the form given above. it is 
possible to eliminate one of the three intrinsic frequencies by making one of the w’s 
zero through a transformation of coordinates. Cf. P. S. Epstein, Ann. der Phys. 51, 
179 (1916). 

82 In writing Eq. (2) W has been used in place of H, as these two letters can be used 
interchangeably because the Hamilton function defined in (20) is identical with the 
energy. - The factor 4 in Eq. (25) is due to the fact that the number of terms in the series 
(1) is half that in (21) because of the requirement in (1) that 711+ 72w2-+ T3ws be positive. 

33 For greater detail concerning the transformation theory outlined in §7 see Appendix 
4 of Sommerfeld’s Atombau or Ch. XT of Whittaker’s Analytical Dynamics. 
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that a field of radiation is present. The system will then no longer be 
conservative and our problem is to compute the rate at which work is 
done by the forces exerted by the light waves. 

As is usual in the study of analogous problems, we shall commence by 
assuming that the incident radiation consists entirely of monochromatic 
waves of frequency », and that the y and z components of the electric 
intensity are zero, while the peak value of the x component will be denoted 
by E,. It is clear that the Newtonian differential equations now become 


mx = —8V/dx—eE.cos[2x(vt+eo)]; my=—AV/dy; mz=—AV/dz; (25a) 
where ¢ is the absolute value of the electronic charge and €» is an arbitrary 
phase angle. 

If we were to study the perturbations in the motion produced by the 
incident wave purely with the aid of Eqs. (25a) it would be impossible to 
make further progress without specializing the form of the potential 
function V, as was done by Planck and Fokker. However, it is quite a 
different story when we seek to compute the perturbations not in x, y, z 
but in the “angle variables” w,, we, ws; and their conjugate momenta 
J, J2, J3, whose properties were sketched in section 7. Because of their 
intimate connection with the orbital frequencies w,, these coordinates 
possess important intrinsic properties. In fact by using them rather than 
x, y, 2, which is the essential feature of the present calculation, the 
periodic properties of the system come to light even without knowing the 
form of V, and it is thus possible to isolate the resonance between the 
frequency of the impressed wave and each combination overtone 7,0;+ 
T2ewWe-+73w3 Of the orbital vibrations. This resonance will be shown to 
depend, in the average, only on the associated amplitude X(r,, r2, 73). 

To get the equations for the perturbations of the J’s and w’s, we need 
simply note that the complete x component of force in Eq. (25a) may be 
written as the negative x derivative of the function V+xeE,cos [2(vt+ 
€9)]. Consequently the equations of motion in Cartesian coordinates 
may be thrown into the canonical form (19) provided in these equations 
we replace H by a new Hamiltonian function H’ which is no longer the 
sum of the kinetic and potential energies but instead has the value 


H' = H + xeE,z cos[2x(vt + €0)] (25b) 


where H is the Hamiltonian function in the absence of the radiation field, 
which is equal to the energy and given by (20). Since by postulate 
(II), the transformation from x, y, 2, pz, Py, Ps to the w’s and J’s is 
a contact transformation, and since it is a fundamental theorem that 
such a transformation preserves the Hamiltonian form even though 
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the Hamiltonian function contains the time explicitly, the equations for 
the w’s and J’s still have the canonical form (23). We must, however, 
utilize in place of H the new Hamiltonian function H’ defined above, but 
expressed now in terms of the w’s and J’s rather than Cartesian coordin- 
ates and momenta. We therefore have 

= = = oa = eE, cos|2x(vt + €0)]; 

dt OJ. OJ & 

dJ « OH Ox 

— = — — — —eE, cos[2n(vt + €)]. 

dt dw, dw, 
Now by postulate (III) concerning the transformation formulas, H is a 
function only of the J’s, so that dH/dw,=0. Consequently using the 
value for x given in (21), and writing cos [27(vt+e)] as the sum of two 
imaginary exponentials, we get 


dw, ‘dH 0A (11, T2, Ts) ( 

ete [ee jexpL2ri(riwr-+ race rates vt) 
texpleri(riertraeetrann—vt—e)]} | (26) 

dJ x : ( : 

rie —rieE, [ tT,A(71, T2, T3) jexPl2ni(rinrt ramet rast vl +er)] 

¢ 


+exp[2ri(r.wit+ ToWo + 73W3—vl— €o)| | ; (27) 


9. POWER SERIES INTEGRATION OF EQUATIONS FOR PERTURBATIONS IN 
w’s AND J’s 


If no impressed electric field were present, the J’s would have constant 
values, which we shall denote by Fa a P th while the w’s would be 
given by w ,t+e, (cf. Eq. 24). Our problem is essentially to compute 
the departures of the w’s and J’s from these values. This can be done by 
the standard method of determining the perturbations of the w’s and J’s 
as power series in the constant E, which represents the amplitude of the 
expressed wave, and which may be assumed to be small. For the results 
which we shall desire, it will only be necessary to determine the pertur- 
bations in the J’s to second powers of E,, and those in the w’s merely to 
the first power. 
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To determine the perturbations in the J’s to first powers of Ez we can 
substitute in the right hand side of Eq. (27) the unperturbed values 
w,t+e, for the coordinates w;, and in this equation we may similarly 
regard the amplitudes A (71,72,73) as having the constant values A °(r1,72,73) 
which result when their arguments J, Jo, J; are equated to t td -. 
Eq. (27) thus becomes 


dJ,./dt = — rick, 4 7A, exp[2ri(e, — €o) | exp[2ri(wt — ny (28) 


where for brevity we have written w, for r1w:+72w2+73ws3, €, for 71€,:+ 
Te€o+73€3 and A, for A%(r;,72,73). 2, denotes the summation with 
respect to 71, 72, 73, andit is to be clearly understood that this r summation 
is a triply infinite summation wherever it appears, in distinction from the 
k, | summations later introduced which involve only three terms. In 
(28) terms are wanting in (w,t+vt) and (e,+ 9), but instead of writing 
these terms explicitly we have inserted a prime above the summation 
sign to indicate that the entire summation is to be repeated with —» and 
— €) everywhere replaced by +yv and + ¢ 9. The symbol 2’ is used with this 
significance throughout, in order to halve the number of terms which 
must be written explicitly. 

We shall suppose the field is applied at ¢=0, so that at this instant 
the w’s and J’s have their unperturbed values. Taking account of these 
initial conditions we have on integration of (28) 





: 1—exp|27i(wrt—vi 
AJi=}eE:) > | Avexp[2ri(e,—e)] pl2mi( uM} (29) 
: T WO,—v 
where =AJx=Ji—J>°. 


Having just found the perturbations in the J’s to first order terms 
in E, we are now in a position to determine the w’s to this degree of 
approximation. We first observe that because of the fluctuations in the 
J’s, the expressions 0H/dJ; appearing in (26) do not have their original 
unperturbed values w;, but may, as they are functions only of the J’s, 
be expanded as power series in AJ;, AJ2, AJ3. We may therefore, to 
the order of approximation under consideration, replace 0H/dJ; by 
We FZ (10:1,2,3)(0W%/I0T1) AJ, where dw;,/dJ; denotes 6°H/dJ,0J; evaluated 
with each J; equal to J;,°. Similarly we may in the right hand side of 
Eq. (26) set wz =w,t+, and replace each expression 0A (11,72,73)/dJ% by 
its unperturbed value, which will be denoted by A ,/8J,. When these sub- 
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stitutions are made, and account is taken of the initial conditions, the 
integration of (26) gives 


. Ow 
Aw, =eEs )> {(2e-s207 ~*) A. exp|2ri(€r—€0)] x 
T l 


2x(w,—v)t—i{ 1 —exp|2ri(w,—»)é] } 


4n(w,—v)? 





OAr 1 — exp|2ri(w, —v)é] 
ij —— exp|2ri(e, — (30 
: OS * Pl padi 4r(w,—v) = 





where Aw, =w;,— wt — €,. 


10. COMPUTATION OF THE PERTURBATIONS AJ; TO 
2ND ORDER TERMS E, 


We can now commence the calculation of the J; to second powers of 
E,. We take the right hand side of (27) and expand the amplitudes as 
power series in the perturbations AJ,;, and similarly develop the factor 
exp [27i(r1w1+72wW2+73Ws)] as series in the Aw’s. We thus get to the 
approximation under consideration 


dJ x . : ; OAr 
eT imeE- 2; Zu- 15293) \r A,+2rir,AwiA, + AJ; ~*~ 
dt Osi 





exp[27i(€r—€0)] exp[2ni(e,—»)]} . (31) 


The values of the AJ’s and Aw’s given by (29) and (30) can now be 
substituted. We have then to deal with the product of two triply infinite 
series, which we shall assume can be multiplied together term by term. 
The only significant terms for our purposes are those which arise from the 
product of a term in 7, 72, 73 into one in —7;, —t2, —73. Other terms 
would involve a factor of the form exp [27i(e,+.,’)] and so would vanish 
when we average over an ensemble of systems having the same initial 
values of the J’s but haphazard values of the phase angles «1, €2, és. 
However, taking this average does not appear really necessary to make 
the other terms mentioned above disappear, as it can be shown that such 
terms vanish if we assume that instead of random waves applied only 
at ¢=0 there are a large number of waves of the same frequency whose 
times of application are haphazard.* 

“If the integration over all values of «, «2, ¢s is omitted (but the integration 
over €9 is retained) there would appear in Eq. (34) “Poisson terms” of the type form 


Rp(w,)tcos[24(w,+w,’)t +7]. This is on the assumption that the field is applied at ¢=0. 
If instead the field were applied at t= — 7, the change in a Poisson term from time 0 to 
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—r3), we may omit such products as involve the factor exp [+47eo], as 
these vanish when we average over a large number of impressed waves 
whose phase angles ¢) are haphazard. The first order terms in E, also 
vanish when such an average is taken. We thus get for the average value 


of dJ ;./dt. 


dJ i : eee | ; Ow, 
mao SB srie°E;” a“, . TA - 1X (i= 1.23)T? -——— x 
OJ 


Furthermore in multiplying a term in (71,72,73) by one in (—71;,—72, 











dt 1 
i| 1— exp|[2ri(w,—»)é] | —2r(w,—v)l exp[2ri(w, —»)é] 
(w,—v)? 
0A_, OAr 1—exp[2ri(w,—v)é] ) 
TDS (1= 2 3) TIA, TiA_.— ig 
OJ; oJ; @,—v ) 


where A_, denotes A(—71;,—72,—73). By combining terms in (w, Fv) 
real terms. (In doing this one must remember to reverse the signs of 
r, and 7, as well as w,.) Using the formula for X? given in (25) and 
performing the time integration we thus get 


with those in —w,+v, this expansion can be written as the sum of only 


" 1°#p 9° Pt r 9 dw, 
AJ «= §e° E722, 4 rk Xr’ ~ (1=15253) TI ounete x 
( OJ; 


| —2+2cos2z(w,— v)t]+22(w,—v)tsin2r(w,—v)t 
———ee + 
(w,—v)8 
0X,?\ 1—cos2r(w,—v)t ) 
tk Do=1.a)T2 : + (32) 
7) l (w,—v)? ) 


In this equation the summation 2,’ is to be taken not over all positive 
and negative values of 7, 72, 73, but over all such values subject to the 











t would be Rp(w,) (¢+T7)cos [24(w, +w,’) (t+T) + y]—Rp(w,)T cos[2x(w,+w,') T+ y] 
and so practically vanishes on the average when we allow the values of T to be 
uniformly distributed in an interval large compared to a period. Besides these Poisson 
terms, there would also be trigonometric and constant terms, but their effect during a 
long time interval is unimportant compared to the “‘secular” or linear variation with the 
time given in (34). It is to be clearly understood that integrating over the time of 
application —T is quite different from integrating over the phase angle « of the im- 
pressed wave, and the former integration is required over and above the latter in order 
to make the cross terms vanish without averaging over €1, €2, €3. 

%* The remarks concerning the vanishing of products may be summarized in the 
statement that the only significant terms result from the multiplication of a term 
wrFvin (31) by aterm —w,;+vin Awor AJ. In forming a term in —w, +v in Eqs. (29) 
and (30) care must be taken to reverse the sign of 7% (in 29) or 7; (in 30), as well as the 
signs of w,and». The sign of ws in (30), however, should not be reversed, as wy does not 
involve the 7’s. 
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restriction that w, be positive. This normalization, which halves the 
number of terms in the summation, is to be understood throughout the 
article whenever the equations involve the real amplitudes X,, D,, etc. 
instead of the complex coefficients A .. 


11. PASSAGE FROM MONOCHROMATIC TO CONTINUOUS RADIATION 


Hitherto we have assumed that the impressed wave is strictly mono- 
chromatic. Actually, however, we must deal with an aggregate of waves 
having a continuous distribution of frequencies. In our calculations the 
passage from the monochromatic to the continuous case will be almost 
identical with the procedure adopted by Planck in his computations for 
the absorption by an oscillator, and for a more detailed or rigorous justi- 
fication for the various steps the reader is referred to Planck’s original 
work,”® 

If p(v) denote the specific energy density of radiation, the amount of 
energy between frequencies vy and »v+dp is p(v)dv. On the other hand it 
is well known that in a homogeneous radiation field, each of the six com- 
ponents of the electric and magnetic intensity contribute equally to the 
energy, so that the mean value of (4?+J/°) /8z, the energy in the medium, 
is 3 62/42. We thus see that for the part &, of the electric intensity due to 
vibrations in the frequency range v, vy+dv, we have é2= $rp(v)dv. Nowin 
our problem ©,=E,cos 21(vt+eo), making é2 =} EZ. Therefore we see 
that E,? is to be identified with the element $2p(v)dv. 

Our problem is now to replace E by $7p(v)dy in Eq. (32) and integrate 
over all values of the frequency » of the impressed light wave. If x denote 
the expression (vy —w,)/, then it is clear that the contribution to AJ of one 
particular term 7), 72, 73 in the multiple Fourier expansion may be written 





ontnnie Ow, [2 —2cos(2rx) —2rxsin(27x) | 
(Mi), =dretn Xe(Zuaran7 t? 


—-— ———-—_—__— —— p(v)dv 
OJ; * 


3 
seis 0x,” [2—cos(2rx)] 
gmerr L(1=16%3)TI t p(v)dx, (33) 
OJ; x? 


where the upper limit in the integration is + ©, and the lower limit is 
strictly speaking —w,/, but if ¢ is sufficiently large may in the usual 
manner be without sensible error treated as — ©. We should in accord- 
ance with the definition of the symbol 2’ also consider in (33) terms which 
come from having —»v everywhere replaced by +», but if now x denotes 
(v+w,)t, these terms would lead to integrals of the same form as those 
given above except with the lower limit +w,f. This means that if ¢ is 
large the integration would be only over very large values of x, and 
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because of the factors x~* or x~* these integrals would be so small as to be 
negligible. This omission of terms in w;+ v is, in fact, customary, and 
simply another way of saying that the only appreciable contributions to 
the perturbations come when the factorsin the denominator approach 
zero, indicating resonance. 

In considering the second integral of (33), p(v) is a function of x, but 
may, following the usual procedure, be replaced by its value p(w;) at the 
resonance point y=w,.*® Since 


” 1—cos2rx 
PE tao steahenl dx =f (sin?rx) x "*dx = 21? 


—— 0 


the second integral in (33) becomes 27°p(w,). A partial integration con- 
verts the first integral in (33) into one identical with the second*’ except 
that p(v) is replaced by dp/dx, which is equal to f"dp/dv. The former 
integral is therefore 27°t—'dp/dv, where dp/dv is to be evaluated at »y=w 

Using the values just obtained for the integrals in (33), and summing 
over all the terms in the multiple Fourier expansion to get the total 
perturbation, we have 








dw,\ dp aD, ) 
AS, = 3 2etd leeD,? (2 15253) TL (20- 15253) T ) oe.) +(34) 
{ OJ iJ Ov OJ; ) 

In writing this equation we have everywhere replaced X,? by the 
expression D,?=X,?+Y,2+Z,2 in order to allow for the perturbations 
produced by the y and z components of the electric intensity. It is obvious 
from symmetry considerations that taken by itself the y component of 
the electric intensity produces perturbations in which X, is everywhere 
replaced by Y,, with an analogous substitution for the 2 component. 
It is also clear that the total average perturbation is the sum of the 
x, y, and z effects taken separately, for cross terms involving factors such 
as E,E, cancel out on the average if the phase relations of the impressed 
radiation are haphazard. Performing this sum amounts simply to the 
substitution of D,? for X,? 

In Eq. (34) it is to be understood that dp/dv and p(y) are both to be 
evaluated with the argument y=w,. This means that these expressions 
have different values for each term in the summation. A similar evalua- 


* Cf, Planck, Wirmestrahlung, 4th Ed., pp. 151 and 164. A Taylor’s expansion of 
p(v) in powers of »y — w, shows that the departure of p(v) from p(w,) may be disregarded in 
evaluating the 2nd (but not the 1st) integral of (33), as the terms thus neglected involve 
powers of ¢ in the denominator. 

37 The integrand of the Ist integral in (33) may be written as pdV/dx where 
V = —(1—cosx) /x? 
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tion of dp/dv and p(v), varying with each term, is also to be understood 
in all subsequent 7 summations. 


12. CALCULATION OF AJ;,AJ; 


The next step is to compute the average value of products of the type 
AJ; AJ), as a knowledge of this kind of expression to 2nd powers of E, is 
necessary for a determination of the absorption and also for the analysis 
in part III. To make the calculation we multiply together the series for 
AJ, and AJ, given in Eq. (29).°* For reasons similar to those explained 
in section 10, the only products which need be retained in taking the 
average are those which arise from the multiplication of a term in w,;}v 
by one in —w,+v. Following a procedure closely analogous to the 
derivation of (32) we thus find 

al y/ r 2 
AJ,AJi = =je°E,* Zz. yee’ ane iF 

The passage from a monochromatic plane wave to continuous homo- 
geneous radiation is identical in character with the transition from (32) 
to (34). In fact the above expression is similar in the form of variation 
with » to the second term of the right hand side of (32), and we can by 
inspection immediately write down the final result analogous to (34) 


1—cos2r(w,—v)é } 





(35) 


AJ,AJ = $ret, | 7.7 .D,*p(w,) (36) 
13. CALCULATION OF ENERGY ABSORBED 


It is now easy to compute the absorbed energy, which we shall denote 
by W—W,. The total energy (kinetic+potential) is identical with the 
Hamiltonian function defined in (20), and by postulate (III) concerning 
the J,w coordinate system, the latter can be expressed in terms of the 
J’s alone, without involving the w’s. Therefore 


W-—Wo= H(J} +A, J3+AJ2, J3+AJ3) — H(J}, J3, J$). 


Expanding in a power series in the AJ’s, and keeping terms of the order of 
magnitude of the square of the electric field, we have 


wes Owk 
W—-Wo= AJ;,)+3 37) 
. ) on ag) oe aR 1523) [ (5) ara] 


where, as previously, we have utilized the fact that w;, = (@0H/0J,)o, etc. 
Bars have been added to indicate that the average is to be taken over a 


38 Only the expressions for AJg and AJ} accurate to first powers of FE, need be em- 
ployed, as we desire their product to squares of Ez. 
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large number of waves. The values of AJ; and AJ,AJ; given in (34) and 
(36) may now be substituted. When this is done, and the summations 
with respect to / are expanded, the formula (37) becomes 





as a a a 
W—Wo= §r'etX,{ p(w,) | 11 — +72 +1r3—— JG, 
| OJ; 





OJ 2 OS 3 
4G ( Owr ‘ Owr 4 ~") } (38) 
==" a T2 is eee ‘ 
Ov ‘a 1 OS » ‘ad 3 J 


where we have written G, for D,w,;. This is the same result as was 


quoted in Eq. (16) of part I.*° 


It is readily verified that the general result given in (38) reduces to 
the values for the absorption derived by Planck and Fokker respectively 
should the system be a linear oscillator or a rotating dipole.*® 


14. SrupDy OF THE SPECIAL CASE WHERE 0p/dv=0 


A particularly interesting situation arises when the radiation field 
is such that the specific energy p(v) is independent of the frequency, 
making dp, dv=0, and enabling one to take p(w,) outside of the summa- 
tion sign in Eq. (38). 


” In Eq. (16) the effect of only the resonance to one particular term in the multiple 
Fourier expansion was considered, and a factor of N was inserted to give the absorption 
by N atoms rather than only one. Also in (16) the time interval was denoted by At 
rather than ¢. These facts immediately explain superficial differences between the 
forms of (16) and (38). 

For an oscillator there is just one frequency, so that terms in 72 and 7; do not enter, 
and we can only have 7;=1 as there are no harmonics. Also dw/dJ is zero as the fre- 
quency is independent of the size of the orbit, while G, is equal to J/2x*m, as is readily 
proved by reasoning similar to that used in the derivation of (12), except that mh is 
replaced by J. Therefore W—W = ze'p(v)t/3m, in agreement with Planck and others 
(cf. Eq. 14). 

For a dipole rotating in a plane we may as above take 7z2= 73 =0, while 7; must equal 
1. Also J=4r°Jw, where J is the moment of inertia. (This becomes obvious when we 
remember that the J which on the basis of the quantum theory is to be equated to nh 
is the angular momentum 277 multiplied by 27.) If] denotes the length of the dipole, 
the motion is of the form x=e~— J] cos 2rwt, y=e—l sin 27wt, so that X,;= Y,;=1/e, and 
G, =2el*e-? = JI? /2n%e21. Therefore by (38) W—W.=42xl2I—[p(w) +(dp/dv)]. This 
is the result obtained by Fokker (Ann. der Phys. 43, 818). (Fokker does not actually 
write down the absorption, but the latter is simply computed from his expression for the 
fluctuation in the angular momentum. In a comparison with the present article, 
allowance must be made for the fact that Fokker uses Heaviside-Lorentz units of charge). 

The preceding two examples illustrate how easily the results for special concrete 
cases can be derived from the general formula (38). This Eq. for the absorption might 
also have been computed as the time integral of —E,ex cos 24(vt+«@). This would be 
easier if we did not have to know the AJ's to second order terms for the work of part III. 
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It will be shown elsewhere by the writer*’ that J, =07/dw,, where T de- 
notes the time average of the kinetic energy, and where in the differentia- 
tion the frequencies w;, we, w; are to be regarded as variables, but all the 
amplitudes X,, etc., are held fast. The formula just cited for J; may be 
regarded as a method for defining the J’s which is alternative and equiva- 
lent to Eq. (2). Now since T =mrd,02D2 and w,=7101+Tewet 73s, 
we have 


Ji. = OT / de, = 2mx*X,} 7G; } (39) 
where G,=w,D,?.. Therefore when p(w,) is a constant independent of 
w,, (38) takes the very simple form 


W—Wo=ze'pt/m . (40) 


Thus we see that when dp/dv=0, the mean absorption of each electron is 
re’*pt/m, regardless of the shape and size of the orbit or the character 
of the dynamical system to which the electron belongs. That is to say, in 
a field of radiation whose specific energy does not vary with the fre- 
quency, we have the rather surprising result that the mean absorption is 
independent of the form of the force function V(x, y, s) which holds the 
electron in the atom.*® * 


| The Fundamental Concepts of the Quantum Theory of Line Spectra, Appendix 1, to 
appear as a Bull. Nat. Research Council. Eq. (39) may also be derived from a somewhat 
analogous relation given by Burgers (viz., Eq. 9, p. 47 of his dissertation, ‘‘Het Atoom- 
model van Rutherford-Bohr’’). 

42 The result (40) also holds for free electrons subject to a field for which dp/dv =0, 
regardless of whether or not free paths are stipulated. This can be seen from the analysis 
given by Jeans (Phil. Mag. 17, 773, 1909) if his expressions for the absorption (his 
Eq. (34) and equation at bottom of p. 790) are integrated over all frequencies, regarding 
p(v) as constant, and taking the dielectric constant as unity. In the case of free paths, 
the result may also be derived from an earlier paper by Sir J. J. Thompson (Phil. Mag. 
14, 223, 1907). If the orbital electron is constrained in such a way as to have only one 
or two degrees of freedom, we must clearly introduce factors of 1/3 or 2/3 into (40). The 
result thus obtained in the case of one degree of freedom is the same as the absorption by 
a linear oscillator in a variable field such that the density p(v) at the resonance point is 
just equal to the constant p employed above (cf. Eq. 14). 

8 It is tempting (but probably futile) to determine the probability coefficient A (r—s) 
in (5) in such a way that a relation similar to (40) is preserved when the quantum theory 
is substituted for classical electrodynamics. The quantum theory analogue of the total 
absorption is the sum of all the different kinds of differential absorption (cf. sections 
4 and 5) of various frequencies which may commence at a given orbit. If in Eq. (10) 
we use the + signs and take m=1, it can be verified without particular difficulty with 
the aid of (39) and (7) that even at ordinary wave-lengths or small quantum numbers we 
have in the quantum theory a relation of the form (40), as the right hand side of Eq. (15) 
becomes ze?pNAt/m when summed over all values of r and ¢. (Such a relation must 
hold asymptotically for large quantum numbers regardless of how A(r—s) is evaluated, 
provided only the correspondence principle for emission is valid.) However, such a 
method of determining A (r—s) is hard to reconcile with the work of F. C. Hoyt (loc. cit. 
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15. COMPUTATION OF POLARIZATION 


In order to furnish a basis for statements made in the discussion of 
dispersion given in section 6, it is necessary to calculate the polarization of 
the atom produced by the incident waves. It is clear from a power series 
expansion of Eq. (21), that to first order terms in the electric field the 
x component of the displacement of the electron brought about by the 
impressed waves is given by* 


Ax= {> (Su +2mir,AwiA )ex [2ri(wt+e,)] 
" “T) Lva=123 Vaz, ~ | a ee 

The values of the AJ’s and Aw’s given by Eqs. (29) and (30) can now be 
substituted. The displacement x is thus expressed as the product of two 
multiple Fourier series. In multiplying these series together term by 
term we need only retain the products which involve the time just through 
a factor exp [2zivt] (meaning that terms, for instance, in exp [27i(w,+ 
w,,+yv)t] etc. may be disregarded). The reason for this is that in con- 
sidering dispersion we are interested in the energy which is re-radiated 
by the atom as light of the same frequency as the impressed wave, and 
not in energy which is dissipated in light of the characteristic frequencies 
of the atom—a sort of fluorescence effect.” This means that the only 
products of significance for our purposes are those which arise from the 
multiplication of a term in w, by one in —w,, and furthermore in multi- 
plying together such a pair of terms we may disregard the constant and 
secular parts“ of AJ or Aw and retain only the portions of the right hand 





in footnote 6), who showed that much better agreement with a limited amount of x-ray 
data is obtained with m=4 than with m=0. The case m=1 would probably resemble 
more closely m=0 than m=4. Another objection to trying to get an exact quantum 
theory analogue of (40) is that in so doing it is necessary to sum over all + and — inte- 
gral values of AJ;,/h (k =1, 2, 3) and some of these would in general lead to transitions 
from positive to negative quantum numbers, which can scarcely correspond to any 
physical reality. Also it is doubtful how much real physical significance can be attached 
to radiation fields such that dp/dv=0. 

“ Higher order terms in the power series expansion would give a change in the 
dielectric constant with field strength which may be disregarded for weak fields. 

* Combinations w,+v of the natural and impressed frequencies also appear in the 
term by term multiplication and would be radiated. However these terms are not of 
direct interest in dispersion, as their frequency differs from ». Frequencies of the form 
w, are of course radiated even in the absence of an impressed field, and the total amount 
of energy emitted is only slightly modified by the addition of the latter. 

“ By “secular” terms are meant those which involve the time explicitly rather than 
through exponentials. Such terms are found in (30) but not (29). The effect of these 
terms in (30) is simply to change slightly the vibration frequencies of the system, rather 
than to introduce new vibrations resonant to the impressed wave. 
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sides of (29) and (30) which involve the time through exponentials. 
Remembering that the prime in the symbol 2’ signifies that the summa- 
tion is to be repeated with —v and —€, everywhere replaced by +» and 
+9, we thus get for the part of the displacement which is resonant to the 
impressed wave*’ 


Av=4eE 2,4 [ (= het ‘) 
V= 7h re, _— ~~ she age ye 
. Welioeaiees ™ oJ oJ 


exp|2ri(vt+ eo) | \ 


W.-Y 








Ow, 
+—-A A-(o,—»)!] 
OJ; 


By combining terms in w,}v with those in —w,+yv and using (25) we 
obtain by a procedure similar to that used in deriving (32) 


A __ e . 0X? )-1 Ow, x? -2 
e=be FD) Qusune | gy, rt + a7 Xe —”) 


where now the number of terms in the + summation is halved by the 
requirement that w, be positive, and where @ has been written for 
E, cos 2x(vt+ eo). 

To get the polarization we multiply Ax by the electronic charge —e 
and the number JN, of atoms per unit volume. On combining terms in 
w,—v and w,+v, and making some slight rearrangements, this gives 








Detexs 41) 
Pee T) oJ (w,X tT ) w,?X ,27;(dw,/dJ1) 
— MEQ) Qocnuee 2(w,?—v?) — (w,? — 2)? 


Besides P., which is the x component of the polarization due to the 
electric intensity &, there would also be y and z components P, and P,, 
involving products of the form A_,B,, A_,C,, etc. where B, and C, are 
coefficients analogous to A, except that they appéar respectively in the 
expansions of the y and z components of the unperturbed motion, rather 
than that for x. This would give an anisotropic effect similar to that 
found in a crystal, as & has by hypothesis the x direction, while P would 
be differently directed. However, this is on the assumption that all N, 
atoms are entirely similar,—i.e. have internal force fields with the same 
spacial orientations. If we assume that the axes fixing the force function 
V(x, y, 2) have random distribution,** then on the average X,?=4D,?, 
while A_,B, = A_,C,=0, making P,=P,=0 and eliminating the aniso- 
tropy. 


In forming a term in —w; in (29) or (30) care must be taken to reverse the sign of 
71 Or Tk as well as just w;. 
“8 Cf. footnote® 
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On making the substitution 4G,=40,D,? for w,X?, expanding the 
] summation, and taking only one term in the 7 summation so as to give 
just the resonance between the impressed frequency and one particular 
combination overtone w,, (41) becomes identical with Eq. (18) of part I. 
This is the desired result.*® 

By reasoning similar to that given in footnote” it may easily be verified 
that for linear oscillators whose axes are parallel to the electric intensity, 
(41) reduces to the well-known result (N,e&/42?m) (w?—v?)-'. For di- 
poles whose axes of rotation have random spacial orientations, the 
polarization is —(N,P?é/122°/) (v?+w*) (w?—v*)-? (notation as in foot- 
note 40).°° 


16. CALCULATION OF DISTORTION OF ORBIT BY FORCES DUE TO 
ELECTRON’S OWN RADIATION 


In the preceding sections we have calculated the perturbations in the 
orbits produced by impressed radiation fields. It is well known that 
there is also a force which is due to the electron’s own radiation. Neglect- 


ing higher powers of 1/c*, this force has*! the value 2e%/3c*, where v 
denotes the vector velocity. From Eqs. (21) and (24) it follows that the 


x component is 
F,=(—16rie?/3c3) 2, {w,*A ,exp|2ri(w,t+e,)| : (42) 


Strictly speaking the force will not have exactly this value, for in differ- 
entiating the expression for x three times to get 7, we have neglected the 


slow alterations in A, and w, produced by the damping due to the 


T 


49 The calculation which we have given of the polarization is the simplest possible 
one. A refinement which should be considered in a more extensive analysis is the effect 
of the damping due to the electron’s own radiation on the form of the dispersion curve 
near the resonance point. Also the modifications appropriate to degenerate systems 
would be of importance, as atomic systems are of this type in the absence of constant 
external fields. (We have avoided the degeneracy difficulty by introducing the artifice 
of internal force fields of random orientations). Considerations of such details are, how- 
ever, beyond the scope of the present paper. 

* It is to be clearly understood that in the computation of the polarization (as well 
as the rest of part IL) it is assumed all orbits have the same size and shape. This would 
not be the equilibrium distribution, as it is shown in part III that this is given by the 
Boltzmann formula, leading to a continuous distribution of sizes and shapes. The latter, 
however, would give rise to. continuous black body spectra rather than monochromatic 
phenomena, and would not explain the occurence of resonance frequencies in dispersion. 
Therefore in comparing the classical and quantum theories to establish correspondence 
principles, we must assume a discrete distribution even in the classical theory. An 
initial discrete distribution would according to the classical theory eventually become 
continuous due to the action of black body radiation and spontaneous emission. The 
present computations apply to a time interval long compared to a period but short 
enough so that the departures from the initial discrete distribution are not sensible. 

Cf. M. Abraham, Theorie der Elektrizitat (4th Ed.), Vol. II, Eq. (87). 
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electron’s own radiation. However introducing these corrections would 
clearly be a second order effect which can be neglected if the time interval 
is not too long.” 

It is clear that the mathematical analysis for the present case resembles 
that for an impressed monochromatic wave, already studied in sections 
8 and 9, except that now the force instead of being due to just one wave 
—eE,cos [27(vt+eo)] is due to the infinity of waves of different frequen- 
cies given by (42).°* Therefore in Eq. (27) or (28) for determining the per- 
turbations we must replace —3eE, { exp [2ri(vt+€o)] +exp | —2er(vt+eo)]} 
by the expression F, given above. We thus have in place of (28) 

dJx%/dt=2nriF ,=, TKA exp|2ri(w,t+er)| ? 

When the explicit value for F, is substituted and the series are multi- 
plied together, the ‘‘cross-terms’’ vanish on the time average, and the 
only terms which remain come from the multiplication of a term in w, by 
one in —w,. Consequently the time average of dJ;/dt is (—32m‘e?/3c*) x 
Y,7,w,3A,A_,. If we utilize (25) and halve the number of terms in the + 
summation by taking only such values of the integers 7;, r2, 7; as make 
the expression w, positive’ we have 


dJx/dt = (—16rte?/3c*)¥,{ rawr? Xr*} (43) 
t 


Integration of this relation gives 
AJ «= t(—169'e?/3c%) Er | riwr®Dr?} (44) 


In this equation we have inserted D, for X, in order to take into account 
also the y and z components of the radiation force. It is clear from sym- 
metry considerations that the effect of these components is simply to 
add terms in Y,? and Z,?. (Cross terms involving products such as 
X,Y, cannot arise in the first order approximation we are seeking, as the 
Hamiltonian function is H—xF,—yF,—2F,, where H is defined as in 
(20); and consequently the x, y, z effects enter separately unless squares 
of the perturbations are considered). 

From (44) it follows that if we considered only the radiation force 
arising from one particular term w, in the multiple Fourier expansion of 
F,, the AJ’s would be in the ratio 

AJ, 2: AJe: AJ3 = 112 72:73. (45) 


52 The interval must in a later part of the proof at the same time be assumed long 
compared with the natural periods of motion, but the addition of this restriction involves 
no contradiction. Cf. M. Planck, ‘‘Warmestrahlung,”’ p. 154 (4th Ed.). 

58 The transition from an equation of type (27) to one similar to (28) may be made 
just as in the monochromatic case, as the coordinates wz and the amplitudes may to 
the order of approximation under consideration be replaced by their unperturbed values 
wet + & and A_, just as was done in section 9. 

% Cf, derivation of Eq. (32). 
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This result has already been quoted in Eq. (11) of part 
The average rate of radiation is the negative of the rate of change of the 
system’s energy and is therefore on the average (cf. section 13) 


—dW/dt=  — wen = (162te?/3c) 2, w,4D,? (46) 
as 0H/dJ,=o% and w,=7,0,+7202+730; This is a well known result 
already utilized in Eq. (8) of part I. If only the combination 
overtone w, were operative in the radiation force given in (42) only 
the corresponding term would appear in the expression given above 
for the rate of radiation. Viewed from this standpoint, the result 
stated in (45) seems almost intuitive; for since dW /dt = (1~1,2,3) w:dJ;/dt, 
it is hard to see how the formula for the radiation of energy could be of the 
well known type given in the right of (46) unless the J’s were altered in 
the manner given in (43) (with X, replaced by D,) and hence in the ratios 
given in (45), when only one term is considered. 


17. GENERALIZATION TO CASE OF MORE THAN THREE DEGREES OF 
FREEDOM 


So far the results have for simplicity been stated only for the case of 
three J’s and hence three intrinsic frequencies w, we, ws. The case of more 
than three of these expressions can arise only when there are more than 
three degrees of freedom and consequently more than one electron or 
moving particle in the system. It is found without difficulty that in this 
case the Hamiltonian function in the presence of a monochromatic wave 
vibrating parallel to the x axis is no longer H+xeE, cos [27(vt+eo)] 
(cf. Eq. 25b) but instead H—2jx,;e;E, cos [27(vt+eo)] where the lsumma- 
tion is to be extended over all the particles in the system, and x; and e; 
denote the coordinate and charge of the jth particle. The analysis can be 
carried through just as for one particle except that —eA,,—eX,, etc., 


%® It is to be noted that the result (45) is not restricted to the case of the particular 
type of force obtained by taking one pair of conjugate terms in (42) but is easily proved 
equally valid for any other force of the type form Scos2rw,t+Tsin27w,t where S and T 
are arbitrary constant vectors. This immediately indicates the possibility of proving 
(45) even when higher order terms in 1/c* and the changes in w, and A, with time are 
considered, all of which were neglected in deriving (42). For when these effects are con- 
sidered and we keep just one combination overtone (terms in w, and — w,)\in the resulting 
expansion for F,, more rigorous than (42), we still have dJ,dt:dJ2/dt:dJs/dt = 11:72:13. 
Since the J’s always change instantaneously in this ratio (neglecting small periodic 
oscillations), we have (45) valid for an arbitrary time interval, even though (44) is an 
approximation valid only for intervals not long enough for the amount of damping to be 
appreciable. 
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must be replaced by the coefficients in the multiple Fourier expansion of 
2,e;x;, which is the x component of the electrical moment of the atom. 
The results of sections 11, 12, 13, 15 are therefore unaltered except that, 
of course, allowance must be made for the fact that there are uw rather 
than three degrees of freedom by adding terms in Ji, Js,. . . Jyu,T4,7s, 

» and in the final equations —eD, must signify an amplitude in 
the multiple Fourier expansion of the polarization of the entire atom 
rather than the displacement of just one charge. 

It is also verified without difficulty that the content of section 16 and 
the correspondence theorems stated in sections 2, 5, 6 of part II can be 
extended to include the general case if only the modifications indicated 
in the preceding paragraph are introduced. 


DEPARTMENT OF PHYSICS 
UNIVERSITY OF MINNESOTA, 
June 19, 1924. 


5% All of the generalizations given above are on the assumption that the dimensions 
of the atom are small compared to the wave-length of the impressed beam. This situa- 
tion has, in fact, been tacitly assumed throughout part II (even when there is only one 
electron), as it is required for the validity of the fundamental Eq. (25a). The case where 
the impressed wave-length is comparable with the dimensions of the orbits could be 
treated after the general methods of the present article, but the calculations would be 
probably rather laborious. 

In atoms with several electrons a case which is readily treated is that where each 
orbit is itself small compared to the impressed wave-length but the distance between 
orbits is large compared to the latter. The electrons then absorb independently whereas 
the generalizations given above are essentially for the case where they absorb in 
unison. The appropriate modification consists in introducing 2,e;D;* rather than 
(2je;D;)* in place of e®D*. This case of independent absorption could not well actually 
arise, as atomic dimensions are usually too small (except for x-rays) and also electrons 
travel all over the atom rather than being confined to small regions or cells. In asystem 
with several electrons the absorption per electron in a constant radiation field is re*pt/m 
(this statement is the generalization of section 14) only if the electrons absorb indepen- 
dently rather than in unison. 

It is interesting to note that when the wave-length radiated is large compared to 
atomic dimensions, the radiation force 2e*#/3c* affects not only the electron for which 
# is computed, but also all other electrons in the atom. For if each electron were af- 
fected by its own radiation force but not that of the other electrons, the work done by 
the forces would be proportional to 2e;?D;* rather than (Ze;D;)*, whereas it is well known 
that the loss of energy should involve the amplitudes in the latter rather than in the 
former fashion if the atom is small compared to the wave-length. 

The assumption that the wave-length is large compared to the orbit cannot be 
regarded as valid when we consider the “‘higher harmonics” or combination overtones of 
extremely large frequency. These terms, however, will not contribute appreciably to 
the total absorption or radiation if the series converges sufficiently rapidly. 

Because of the substitution of 2;(e;D;)* or 2je;*D;* for e*?D* in (9) or (10) it should be 
noted that in atoms with more than one electron the probability coefficient Ar—, 
depends in general on the amplitudes of all the different electrons rather than just one 
of them. 
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THE VARIATION OF THERMIONIC EMISSION WITH 
TEMPERATURE AND THE CONCENTRATION OF 
FREE ELECTRONS WITHIN CONDUCTORS* 


By A. T. WATERMAN 


ABSTRACT 


Equilibrium electron theory of thermionic emission.—The equation for the 
thermionic saturation current per unit area in the form J=AT*e~/T js derived 
according to Richardson’s thermodynamical method, assuming (1) the bound- 
ary condition that when y the free electron evaporation constant is zero, the 
free electron concentration inside and outside the hot body is the same, and 
(2) the interior free electron concentration given by the author’s equilibrium 
theory of conduction (Phys. Rev. 22, 259, 1923). a depends only upon the 
valence in the electron-forming reaction and always lies between 1/2 and 5/4. 
8 depends upon the valence and the two evaporation constants ¢o and yo, refer- 
ring respectively to a bound and a free electron. A is found to depend upon 
the valence, the electronic chemical constant and the concentration of electron- 
forming particles in the hot body, varying from 4.16(10)" for BaO to 8.13(10)" 
for Ni. Values of 8 calculated from the computed A and experimental data 
show excellent agreement with 6 observed for Pt, W, Mo, Ni, Ca, CaO, SrO, 
and BaO, the verification being more general than by use of Dushman’s 
formula. The experimental data however are not sufficiently accurate to decide 
between the two formulas. The failure of certain hot bodies to agree with 
either formula is discussed and possible explanations offered. Assuming the 
theory correct, the free electron concentrations in several metals, MoS, and CaO 
are calculated from experimental data. For metals at 0°C this concentration 
is of the order of magnitude of 10!” per cc, and decreases with rising tem- 
perature. A noteworthy consequence of this low concentration is that the 
contribution of the free electrons to the specific heat of the conductor remains 
negligible down to very low temperatures. 


HE variation of the thermionic saturation current density with 
temperature is given empirically by an expression of the form 


—b/T 


. 
I=AT e (1) 


where the exact value of \ matters little in obtaining agreement with 
experiment, provided it is of the order of magnitude of unity. The well- 
known theoretical expressions are! 


/kT 


~v 
IT=A,T'"e (2) 
* The main conclusions of this paper were presented before the American Physical 
Society on Dec. 27, 1923 (Phys. Rev. 23, 299, 1924) 
1 Richardson, Emission of Electricity from Hot Bodies, p. 36. y is used to denote the 


thermionic work function commonly denoted by ¢, since in this paper ¢ is reserved for 
the photo-electric work function. 
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and T=A,T? eve*T (3) 


derived respectively from kinetic theory and from thermodynamical 
considerations. The temperature functions in the two expressions become 
identical when the thermionic work function y is assumed to depend upon 
temperature according to the equation? y=yYo+(3/2)k7, also derived 
from thermodynamical reasoning. The same form of equation has been 
obtained by application of the quantum theory.’ Although the experi- 
mental evidence has thus far not been able to distinguish between Eqs. 
(2) and (3), Richardson and others have agreed that probably Eq. (3) 
deserves more confidence from the theoretical standpoint. 

The evaluation of the constant A is a matter of interest and importance. 
Differences in the form of this constant arise chiefly from the different 
methods by which the concentration of external electrons in equilibrium 
with the hot body is deduced. Thus we have 

vyskr 
Na = Npe (4) 
and 
Na =CTI%¢" ois (5) 
corresponding respectively to Eqs. (2) and (3), where m4 =equilibrium 
concentration of external electrons, mg =electron concentration within 
the hot body, and C involves a constant of integration. 

Dushman,‘ by application of the Nernst heat theorem, has recently 
derived an expression for the thermionic saturation current density 
of the form of Eq. (3), by evaluating the constant C in Eq. (5), with the 
interesting result that C, and hence A», should be universal constants. 
However the premises on which this latter derivation is based may pos- 
sibly be open to question, and furthermore the universality of A» does 
not appear to extend to all possible emitting substances. The author 
has submitted’ an expression for the thermionic saturation current 
density which is also consistent with Eq. (1). The previously omitted 
details of derivation are presented below, and the work has been extended 
to include evaluation of the integration constant appearing in the coeff- 
cient A (Eq. 1). This makes possible a numerical calculation of the value 
of A from data independent of thermionics, and thus permits a more 
substantial examination as to the validity of the formula proposed, as 


? Richardson, loc. cit,’ p. 35 

’ Richardson, loc. cit.’ p. 41. 

* Dushman, Phys. Rev. 21, 623 (1923) 

5 Waterman, Phys. Rev. 22, 259, 268 (1923) 
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well as allowing direct comparison with values of A derived by other 


methods. 


DERIVATION OF EXPRESSION FOR THERMIONIC SATURATION 
CURRENT DENSITY 


From the thermodynamic reasoning® 


dn/n,= (gs)er (6) 


where y is the energy required to liberate one free electron from the hot 
body ;® & is the gas constant per electron. Hence log ng =f (W/kT?)dT+ 
const. Using the boundary condition that, when y =0, the concentration 
of free electrons inside and outside the hot body is the same, or m4 =”,, 
we have 


Since? y=Yo+(3/2)kT 
_ —Wo/ kT 
na = npT?!e : (8) 
From the writer’s equilibrium theory of electrical conduction’ the 
internal equilibrium concentration of free electrons mg is given by 


np= BT —3¢/2(¥+1) ev (de—Po)/ (+1) kT (9) 


where ¢» and yp are respectively the energies (at 0°K) required to liberate 
a bound and a free electron from the conductor, whence ¢)— yo is seen to 
be the energy (at 0°K) required to change a bound electron into a free 
electron within the substance; 

v=valence in atomic reaction: Atom $ positive ion+y» electrons; 

B=c (yN) Vet) ; 

N=concentration of nuclei within conductor (or total concentration 
of particles in conductor from which electrons may be or have been 
emitted). 

The constant of integration c may be evaluated by noting that the inte- 
gration constant of Eq. (3) in the writer’s article’ is for the general atomic 
reaction equal to vi, where 7 is the chemical constant of the electron, 
reckoned in absolute units per electron and in terms of concentrations 
instead of partial pressures. 

Then c becomes e”/+), whence B = ert/@+1) (yN)1/0+1) 


* RicLardson, loc. cit.' p. 31, Eq. (6) 
? Waterman, loc. cit. p. 259 
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Substituting the value of ng from Eq. (9) in Eq. (8) 


F aie 3 vdotyvo 
ng=ertl (vN)*+1 T2041) @~ tk | (10) 





giving the concentration of external electrons in equilibrium with the 
hot body. This concentration is seen to depend upon the concentration 
of nuclei or electron-forming particles within the substance, the valence 
in the atomic reaction, and the chemical constant of the electron, in 
addition to the temperature and the evaporation constants. | 

Following the customary procedure, the thermionic saturation current 
density now becomes 


I=ATe-8/T (11) 
where 
A=er/(h/2xm) er/>+)(yN)V/04+) (€ being the electronic charge); 
a=(v+4)/2(v+1) ; 
B= (vpot¥o)/(U+1)k . 


For different valences v=1, 2, 3, the values of a are respectively a=1.25, 
1.00, 0.875, and thus the exponent of T for all valences is seen to lie 
between the usual 1/2 and 2. 

It should be observed that the constant 8 in Eq. (11) as experimentally 
determined does not necessarily yield directly the energy required to 
liberate a free electron from the substance (¢» as usually understood), as 
is true for } in the ordinary expression for the thermionic saturation 
current.6 Instead 8 involves both the energy required to liberate a free 
electron and the energy required to liberate a bound electron. It is clear 
that such must be the case if the external concentration of electrons 
depends upon the internal concentration of free electrons and if the latter 
in turn is determined by a chemical equilibrium between atoms, positive 
ions and free electrons. Thus an amount of energy (¢o9—yo) is required 
to change a bound electron into an internal free electron, and an amount 
of energy yo is then required to liberate this electron from the substance. 


EVALUATION OF CONSTANTS AND APPLICATION TO EXPERIMENTAL DATA 


Using values of constants as follows: 
«=4.774X (10)-"° e.s.u.; m (electron) =8.995 x (10)-*8gm. 
k =1.372 X(10)-"* ergs/deg. C; i= 34.98; 
we have log A = —4.129+15.19 »/(v+1)+[1/(v+1)] log rN. 
In view of the uncertainty attaching to the theoretical values of the 
chemical constant of the electron, the value of the chemical constant 








A. 


TABLE I 
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used is that selected by Tolman® from experimental data, expressed in 
terms of concentrations instead of partial pressures and in absolute units 
per electron. N, the concentration of nuclei within the conductor, may be 
determined by (density/atomic weight) Avogadro’s number.’ The 
valence vy is the only arbitrary constant appearing in the expression, and 
for most substances agreement with experiment is obtained for v=1 
(i.e. a univalent reaction liberating electrons within the conductor). Thus 
the value of the constant A may be calculated, after which an experimen- 
tal value for the thermionic saturation current density at any single 
temperature enables 8 to be calculated. This computed value of 8 may 
then be compared with the value of 8 obtained directly from the slope 
of the graph (log J=a log T) vs. 1/T. The following table shows the 
results of this procedure for a number of emitting substances. 








Dushman 
B(calc.) B(obs.) B(calc.) Blobs.) yo 





4 

q 
a 
$ 
y 





Observer a A (calc.) 
Davisson .25 7.16(10)" 57,300* 53,600* 51,950* 52,000*(volts) 
and Germer 56,800¢ 54,600t 
1.17 2.98(10)* 54,100* 54,000,* 4.66 
K. K. Smith 1.25 7.16(10)" 52,150 53,700 53,900 51,400 
Deininger a4 0.85 4.71(10)! 47,000 47,000 55,000 46,700 
Horton 2 1.00 3.57(10)" 41,100 41,000 43,800 38,800 3.54 
Stoeckle 1 1.25 7.35(10)" 52,800 52,300 47,850 51,200 4.51 
Schlichter 1 1.25 7.28(10)" 53,100 52,500 49,150 50,000 4.53 
Deininger 1 <3 e 51,100 50,000 4.31 
Schlichter 1 1.25 8.13(10)"% 51,600 51,500 47,900 49,000 4.44 
Jentsch 1 1.25 5.33(10)" 40,000 39,500 37,000 38,900 1.6 
Deininger * ss 40,750 42,700 (do =5. 3) 
BaO Jentsch .25 4.16(10)" 41,000 41,000 38,000 39,700 
” .25 4.24(10)" 44,250 44,200 41,100 42,200 





measured. 


* Worthing and Forsyth temperature scale; 
+ Langmuir temperature scale. 


limits of experimental error. 


throughout the temperature range. 


The value given for @(calc.) in the case of any substance is the mean of 
separate calculations of 8 for each temperature at which emission was 
Naturally, where agreement is close between #(calc.) and 
B(obs.), these separate values of B(calc.) showed a constancy within the 
Where a discrepancy occurs it is evident 
that the individual values of @(calc.) showed a progressive variation 


Where the accuracy of the experimental data seemed to warrant it, the 
current density was corrected for the expansion of the emitting surface 


§ Tolman, J. Am. Chem. Soc. 43, 1592 (1921) 
* In the case of the oxides the concentration of mole cules was used. 
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with temperature, since this may cause a difference in 8 (calc.) of 1 to 2 
per cent. 

It is seen from the table that a valence vy=1 gives very satisfactory 
agreement for all hot bodies listed except for W, Ca and Ta. For W exact 
agreement is secured with v=1.25, i.e., some of the atoms may ionize 
doubly, giving a valence between 1 and 2. Ca requires a valence of 2, 
which appears reasonable, and for Ta v apparently lies in the neighbor- 
hood of 3. The experimental data for Ta given by Deininger are not 
satisfactory however, in that the graph (log 7—a log 7) vs. 1/T is not a 
straight line for any value of a. It may furthermore be noted that the 
values of A which give excellent agreement between 8 (calc.) and B (obs.), 
vary from 8.13 (10)" for Ni to 4.16 (10)" for BaO, even though the 
valence is the same (unity). However, as far as the effect of this variation 
in A on the value of 8 (calc.) is concerned, A might be considered approxi- 
mately a universal constant for those cases where y=1, since taking the 
mean of A (calc.) fot such cases would change the calculated value of 8 
by only +2 per cent at the most. The obvious reason on this hypothesis 
for this approximate constancy of A for the same valence lies in the fact 
that the concentrations of nuclei in different elements and of molecules 
in compounds are closely of the same order of magnitude. However the 
use of the valence factor makes it possible to extend agreement to other 
substances, in which case a large change in A occurs. 


COMPARISON WITH OTHER FORMULAS 


For comparison the values of 8 (calc. and obs.) are included as com- 
puted according to Dushman’s expression,‘ in which a=2, and A is the 
universal constant 1.80X(10)" in the units employed in this table 
(e.s.u./cm (deg.)*). 

From the results contained in Table I it appears that the expression 
here proposed for the thermionic saturation current density compares 
favorably with that of Dushman, and in fact with any expression where 
A is taken to be a universal constant. It should be noted that Dushman’s 
value of A for best agreement is based upon the chemical constant of the 
electron computed from the Sackur-Tetrode expression. The values of 
A computed in this paper however are based upon Lewis’ and Gibson’s 
calculation from the experimental determination of the entropy of 
helium, following Tolman.* The use of the chemical constant derived 
from Lewis-Gibson-Latimer theory’? or the one selected by Tolman, when 


10 Lewis, Gibson, and Latimer, J. Am. Chem. Soc. 44, 1008 (1922) 
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used for calculation of A in Dushman’s expression, brings increasing 
disagreement between 8 (calc.) and 8 (obs.). As a matter of fact, the 
value of A in Dushman’s equation which gives the best general agree- 
ment for all substances in Table I (omitting Ta) is about 5x(10)" 
instead of 1.8 (10)". The use of the other two values for the chemical 
constant results in values of A of 1.53(10)" and 1.15 (10)" respec- 
tively. In the same units the value of A calculated from Richardson’s 
application of the quantum theory is 1.5 x (10)?°. 

In justice to Dushman’s equation it should be stated that, of the 
various substances listed, it shows its best agreement for the most reliable 
data. The expression derived in this paper on the other hand shows 
substantially the same agreement in these cases, but has a more extensive 
application. 

From the theoretical point of view, aside from the evidence mentioned 
regarding the choice of chemical constant, the proposed formula for the 
saturation current seems to rest upon a more satisfactory physical basis, 
if we leave out of account the particular expression used for the internal 
electron concentration. On the basic assumption employed in derivation, 
quantitative results cannot be obtained unless some expression for this 
internal concentration is assumed. This very fact however may be ad- 
vantageous in furnishing further evidence concerning the order of magni- 
tude and variation with temperature of the electron concentration in 
conductors. Whether the expression here used for this internal concen- 
tration of free electrons, based upon a chemical equilibrium relation, 
is justified or not, can only be decided by the success of its application to 
this and to other electrical phenomena. The measure of its success in 
explaining the general temperature variation of electrical conductivity 
for both good and poor conductors has already been discussed.® 

On the other hand the assumptions underlying Dushman’s thermo- 
dynamical derivation may be questioned. Thus the hypothesis that 
thermionic emission is completely analogous to evaporation from a solid 
or liquid is in all probability merely an approximation. As the electrons 
are held in the solid only by virtue of the solid atoms, the case would seem 
to be more strictly analogous to the emission of absorbed gas from a solid 
or of dissolved gas from a liquid. In the latter case the equilibrium pres- 
sure of the external gas is directly proportional to the concentration of 
gas dissolved in the liquid. This condition is satisfied for the analogous 
case of electron emission in the present derivation whereas according to 
Dushman’s conclusions the equilibrium pressure of an electron gas should 
be independent of the electron concentration within the emitting body. 
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Furthermore, in Dushman’s derivation, the specific heat of the internal 
electrons c, is assumed to be negligible as compared with that of the 
external electrons. However, the most probable reason why the internal 
electrons do not contribute appreciably to the specific heat of a metal, 
if there is any truth in the application of the kinetic theory to the internal 
electrons," is that this concentration is small compared with the con- 
centration of atoms. It is not beyond the bounds of possibility that the 
inclusion of the term involving c, might produce an appreciable correc- 
tion to the thermionic current as deduced by Dushman, even though the 
contribution of c, to the specific heat of the hot body appears to be 
negligible. In this way the constant A of the saturation current, Eq. (3), 
might be found to vary slightly for different substances, and thus agree- 
ment be obtained between the two expressions Eqs. (3) and (11). 


VALUES OF THE WORK FUNCTIONS @9 and po 


In the last column of Table I are given various values of the energy 
functions ¢o and yp which appear in 8. As stated above Wo represents the 
energy required (at 0°K) for emission of one free electron from the sub- 
stance, while ¢o is the energy necessary for emission of one bound electron. 
¢o has been identified’ with the photo-electric work function, and Wo with 


the work function entering into the usual theory of thermionic emission. 
In the case of metals calculation of ¢9 and ¥> from 8 may readily be made, 
since’ the conductivity evidence indicates that ¢o and yY» do not differ 


by more than 1/10 percent. Thus B= (vdo+wWo)/(v+1)k =Wo/k (or o0/k) 
for these substances. But for poor conductors like the oxides ¢o is sub- 
stantially greater than Yo. However, if from data on conduction, 


1 Note added July 19: It has been shown by Davisson and Germer (Phys. Rev. 20, 
329, 1922) and recently emphasized by Richardson (Proc. Roy. Soc. 105, 387, 1924) 
that in the accurate work of the former on pure tungsten the values of the thermionic 
work function y (in the present notation) as obtained simultaneously by calorimetric 
and by thermionic methods agree only when the thermal energy of the internal electrons 
involved in thermionic emission is negligible. This is naturally accepted as important 
evidence against the classical free electron theory. It should be pointed out, however, 
that the argument as presented may fail if the constant A of the thermionic equation is 
a function of the temperature, which presumably is the case if the concentration of 
internal free electrons depends upon temperature. In particular, using the function of 
temperature employed in this paper for the internal free electron concentration the 
above conclusion is exactly reversed. Thus ¥ = ¢) =4.66 v (cf. Table I) for W, whence 
¥=Yot (3/2)kT =4.664+-0.29 =4.95 v at 2270°K. On attributing to the internal elec- 
trons the classical value for their thermal energy this number should agree with Davisson 
and Germer’s calorimetric value of ¥, without correction, namely 4.91 v. The agree- 
ment is seen to be to less than 1 percent on the classical theory, while according to the 
non-classical point of view, ¥(thermionic) = 4.66 v and (calorimetric) =4.52 v, a dis- 
crepancy of 3 per cent. 
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b =v(do—Wo)/(v+1)k of the expression for specific resistance is found," 
then the independent values of ¢) and ¥) may be computed. For CaO, 
Horton’s data® give b = 21,000 approximately, which taken in conjunction 
with 6 = 40,000 results in the values of @9 and Wo given in the Table. 


DISCUSSION OF LACK OF AGREEMENT FOR Most OXIDES 


In the case of oxides other than those listed (i.e. other than those of the 
alkaline earths) data by Jentsch' show lack of agreement in most in- 
stances between A(calc.) and B(obs.). In every such case B(obs.) is less 
than A(calc.), and this statement holds for any of the above-described 
methods of calculating the constant A. High values of the valence, 
according to Eq. (11), lessen the discrepancy but do not eliminate it 
entirely. It is possible that this lack of agreement may not mean a 
breaking down of the formula, but that it is due to secondary effects. 
Thus, as Richardson has pointed out," in the variation of y with temper- 
ature a term involving the specific heat of electricity was neglected which 
may not be justifiable when dealing with most poor conductors. Further- 
more in the case of emission from a poor conductor like an oxide, the latter 
may find difficulty in replenishing rapidly enough the supply of free elec- 
trons near its surface. Thus the condition for current saturation might be 
hard to fulfill in the case of emission from high resistance material. Diffi- 
culty in attaining saturation is common experience in such investigations. 

The fact that the discordant oxides have melting points close to the 
temperature range of measurement, while the alkaline earth oxides for 
which agreement is excellent have melting points about 1000°C higher 
than the region of observation, suggests the possibility that the process 
of vaporization may affect the magnitude of thermionic currents. Since 
the rate of vaporization follows a law similar to that of thermionic 
emission, no change in form of the thermionic equation might be expected 
if this be true, but its constants would be altered. If vaporization should 
facilitate electron emission this would tend to explain the results ob- 
served. 

The state of poor conductors like the oxides, particularly as concerns 
both conductivity and thermionic investigations, is often apparently an 
extremely sensitive function of temperature, applied potential, and other 
physical conditions, and it is frequently difficult to distinguish between 
reversible and irreversible changes which occur. Under these circum- 

12 Waterman, loc. cit.’ p. 262, Eq. (7) 

13 Horton, Phil. Mag. 11, 505 (1906) 


4 Jentsch, Ann. der Phys. 27, 129 (1908) 
4 Richardson, Electron Theory of Matter, p. 452-454 
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stances it becomes questionable in many cases to expect that any simple 
theoretical expression may be verified until the nature of such peculiari- 
ties is more fully understood. 


CALCULATION OF FREE ELECTRON CONCENTRATIONS WITHIN 
CONDUCTORS 


To summarize the above, in Eqs. (4), (5), and (10) we have different 
expressions for the concentration of free electrons in equilibrium with a 
hot body. Eq. (4) yields a value of A consistent with experiment pro- 
vided n,, the concentration of internal free electrons, is of the order of 
magnitude of the atomic concentration within the substance. Such a 
magnitude is difficult to reconcile with other phenomena, notably the 
contribution of the electrons to the specific heat. According to Eq. (5) 
the concentration of internal electrons does not enter into the problem, 
the work function only being characteristic of the hot body. For the 
reasons above stated this appears merely to be an approximation to the 
truth. It has been shown that Eq. (9) is fully as satisfactory as Eq. (3) in 
accounting for the thermionic current. Assuming the validity of the 


hypothesis underlying Eqs. (8) and (9), it is therefore possible to calculate 


the concentration of free electrons within a conductor from Eq. (7). 
These free electron concentrations are given in Table II, which includes 
the variation of the concentration with temperature for Na. The reason 
that the concentration in metals increases with decreasing temperature 
lies in the fact that conductivity considerations in the case of metals 
require that Wo be slightly greater than @¢o, the energy of dissociation 
(¢)—Wo) of a bound into a free electron being negative.’ For poor con- 
ductors similar considerations require that ¢)—yWo be positive and large. 
. The valence v in the atomic reaction liberating electrons is taken to be 
unity for all substances listed. 


TABLE II 

Conductor N (per cc) T (°K) b (°K) n (per cc) 
Na 2.6107" 273 —26* .0Ox (10)" 
173 .5X (10)? 

100 .6X (10)"” 

50 .6X (10)"" 

20 .4 (10)'8 

10 .5xX(10)'* 

273 .7X(10)" 

273 5X (10)"" 

273 .4X (10)!7 

273 .1X(10)"” 

273 .0xX (10) 
273 .3X (10) 

773 .4X(10)" 

$.3 1273 +21,000* 2.3 (10)* 


*from conductivity data; + assumed negligible. 
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The free electron concentration thus appears lower than most estimates 
given hitherto. Such an order of magnitude however seems demanded by 
the study of the specific heats of good conductors. It will be observed that 
the mean concentration in metals at ordinary temperatures, namely 
about 10" per cc, gives an entirely negligible electronic contribution to 
the specific heat of metals, and in the case of Na, which is typical, the 
specific heat contribution by free electrons remains negligible down to 
very low temperatures, being 0.0024 at 50°K out of a total of 3.50 cal. 
per gm atom, the atomic heat of Na at this temperature. Even at 10°K 
this contribution would only be 0.028. This result appears to be of 
especial interest and significance, in that these concentrations, which are 
computed directly from the theoretical expression, remove what has 
always been a serious obstacle to the progress of a free electron theory. 

Concentrations of this order of magnitude necessitate large mean free 
paths according to the usual simple classical conductivity theory. In 
the case of Na the electronic mean free path, as calculated from its 
specific resistance by this method, would be 0.014 cm at 0°C. On the 
other hand it becomes increasingly probable (cf. results of Ramsauer"®) 
that such a mean free path is possible. Bridgman’ has already concluded 
that the mean free path in metals is at least 100 or 1000 times as long 
as the mean distance between atoms. In fact a long mean free path, 
coupled with the fact pointed out by Eldridge!® that an electric current is 
carried by the free electrons whose thermal motions have a component 
perpendicular to the current flow, may have a share in explaining the 
abnormally high resistance of thin films. 

In conclusion, the hypothesis that in conductors a free electron con- 
centration exists which is determined by chemical equilibrium considera- 
tions, yields a value of this concentration as a function of temperature 
which is competent to explain the observed temperature variation both of 
electrical conductivity and of thermionic saturation current densities. 
Furthermore the computed numerical values of these free electron con- 
centrations are not in conflict with a most important requirement of any 
free electron theory, namely that the equi-partition value for the specific 
heat of the free electrons be negligible. 

SLOANE Puysics LABORATORY, 


YALE UNIVERSITY, 
April 15, 1924 


16 Ramsauer, Ann. der Phys. 66, 546 (1922) 
17 Bridgman, Nat. Acad. Sci. Proc. 7, 299 (1921) 
18 Eldridge, Phys. Rev. 21, 2, 131 (1923) 
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CHANGE OF CONDUCTANCE OF SELENIUM DUE TO 
ELECTRONIC BOMBARDMENT 


By RALPH DE LAER KRONIG 


ABSTRACT 


Effect of electronic bombardment on the conductance of selenium.—A low- 
resistance selenium cell about 1 inch square made by pressing Se between 
twenty double turns of fine Pt wire, was mounted in a vacuum tube with a grid 
and an oxide-coated, equipotential, electrically heated cathode and subjected to 
bombardment by electron currents up to 400u-amp. at potentials up to 95 volts. 
The increase of conductance AC with current J and with voltage V, is found 
to agree well with the theoretical formula AC/C=./1+JV—1, derived on 
the basis of the free electron theory of conduction, assuming that the impinging 
electrons increase the number of free electrons throughout the whole volume 
of the Se by ionizations proportional both to the number of electrons and to 
the energy of each, ic. to JV. The constant & is the ratio of the ions created 


per watt to the number of ions produced spontaneously per sec. and comes 
out 98.5. 





THEORY 


CCORDING to the electron theory of metallic conduction the cur- 

rent is carried by the free electrons in the conductor, and the con- 
ductance is directly proportional to the number of free electrons per 
unit volume. After the conductor has been left sufficiently long under 
unchanged external conditions, a state of equilibrium must exist between 
the un-ionized atoms, the ionized atoms, and the free electrons; i.e. the 
number of free electrons produced per unit volume every second must be 
equal to the number of free electrons recombining per second per unit 
volume. 

In the case of selenium only a very small fraction of the atoms ever 
becomes ionized, and the number of un-ionized atoms per unit volume can 
therefore be regarded as constant. Let ” be the number of ionized atoms 
per unit volume, and therefore the number of free electrons per unit vol- 
ume. The chance of a recombination is proportional to n?, and therefore 
the number of recombinations per unit volume per second may be put 
equal to an’. 

In the case here dealt with, of selenium subjected to bombardment by 
electrons, in addition to the ions produced by thermal agitation and by 
the constant illumination from the filament, whose number per unit 
volume per second we shall denote by az, we may suppose ions are pro- 
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duced by the impinging electrons uniformly throughout a surface layer 
of depth d. Let the conductance of this surface layer be changed by the 
bombardment from Cy to C;. It is reasonable to assume that the number 
of ions produced is proportional to the number of impinging electrons, 
i.e. to the current J falling on the cell, in future simply called the plate 
current, and also proportional to the energy of each electron, i.e. to the 
difference of potential V through which it has fallen and which we shall 
speak of as the plate potential, therefore proportional to the product JV. 
For equilibrium in the surface layer we must have, therefore, 


ayn? = dota3lV 





n=/(deta3IV)/ay . 
Therefore the conductance of this layer is 
Ci=VbitblV 


where 0, and b2 are constants. When no electrons are falling on the cell 
other conditions remaining unaltered, the conductance of this layer is 


Co= Vb1, 
so that 
Ci=VCe+belV . (1) 


For the ratio of the change in conductance AC to the original conductance 
of the whole cell C, we have 


AC Rs Co _ Co 
> + bl “—* Vk+ kel V— hy (2) 


where k, = ‘ans and ke is a constant independent of J and V. 
If the product JV is sufficiently small, we may write, 


AC/C=3(Ro/ki)IV . (3) 
APPARATUS 


The selenium cell used had a dark resistance of about 3500 ohms. It 
was made by winding 20 turns of each of two platinum wires, about .01 
cm in diameter, side by side on a slab of insulating material about 2.5 cm 
square, and pressing the selenium between the wires on one side. A 
shellac coating on the other side prevented the wires from touching. The 
source of bombarding electrons consisted of a heating strip of platinum, 
.025 mm thick, 3 cm long, and .1 cm wide and an outer platinum strip, 
.05 mm thick, insulated from the inner one by a sheet of mica. The outer 
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strip was coated with a half and half mixture of barium oxide and stron- 
tium oxide and connected to one terminal of the inner one, thus forming 
an equipotential emitting surface when a heating current flowed through 
the inner strip. Cell C and filament F together with a grid G were moun- 
ted on the lead-in wires passing through the ground glass stopper of the 
glass tube (Fig. 1). The grid was 3.5 cm on edge, having approximately 
2 wires per cm. The tube was kept evacuated by means of a rotating 
mercury pump and a mercury diffusion pump in series, a mechanical 
pump serving as fore-pump. The mercury vapor was frozen out in a trap 
surrounded by carbon dioxide snow or a mixture of carbon dioxide snow 
and acetone, the last traces being removed by a piece of gold foil inserted 
in the exit of the vacuum tube. Th‘s was found to be very important, 
not only because the vapor would give rise to ionization at high plate 


gold foil 
S 


F 


Fig. 1. Diagram of vacuum tube containing selenium cell C, grid G 
and Wehnelt cathode F. 


potentials, but also because it affects the surface of the selenium, greatly 
increasing its conductivity. A stop cock S was therefore inserted in the 
tube leading from the vacuum chamber and was opened only after the 
cooling mixture had had sufficient time to remove the mercury vapor. 
The pressure in the tube was read on a MacLeod gage; it ranged around 
10-° mm. The outside light wag-swcluded by means of a black cloth so 
that the cell was exposeaOhly: to the constant illumination from the 
filament. The voltage Hetweertcetigand filgwfént could be varied from a 
potentiometer congétted to a stor due attery. An adjustable voltage on 
the grid regulated the plate current without altering the filament current. 


The cell was inserted in one arm of a Wheatstone bridge in parallel 
with a known resistance r, whose midpoint was connected through the 
plate voltage supply to the equipotential filament, a galvanometer or 
microammeter A serving to read the plate current (Fig. 2.) and a volt- 
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meter V indicating the plate potential. The resistance r was necessary in 
order that none of the plate current might flow through the bridge gal- 
vanometer. Furthermore, in order that the electrons striking the cell 
should all have the same velocity, say within .1 volt, the drop across the 
cell due to the bridge current was kept as small as possible, never exceed- 
ing .1 volt. It was on this account that a low-resistance cell was used, 
because with so small a voltage no precision could otherwise have been 


obtained in the resistance measurements. 
— —_ 
4A 
Vi 
| 


Fig. 2. Electrical circuits used in measuring the resistance of the selenium cell, 
and the plate current and the plate and grid voltages. 
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RESULTS 


The resistance of the combination C and r in parallel was measured 
for different constant plate currents, varying the plate voltage up to 100 
volts, and from this resistance X’ the resistance X of the cell was com- 
puted, and then the ratio AC/C plotted as a function of the plate voltage 
V. Fig. 3 shows three characteristic curves for constant plate currents of 
50, 150, and 400 micro-amperes respectively. Similarly Fig. 4 shows three 
characteristic curves obtained by keeping the plate voltage V constant 
and finding the ratio AC/C as a function of the plate current J, the curves 
corresponding to values of V of 20, 50, and 80 volts respectively. From 
the curves it is evident that for sufficiently small values of J and V the 
curves are practically straight lines and can therefore be represented by a 
formula of the type of Eq. (3), 

AC/C=3kIV 
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For larger values of plate currents and plate voltages the values found for 
AC/C are in good agreement with the formula 


AC/C=/1+kIV-1 (4) 


ca 
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Fig. 3. Ratio AC/C as function of the plate voltage for constant plate 
currents of 50(A), 150(B) and 400(C)y-amp. 






































which is identical with Eq. (2), if k:=1. For if the values of AC/C, /, 
and V are used to compute k, it comes out constant over the range of 
observations, (see Table I). We therefore find that k, of the theoretical 
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Fig. 4. Ratio AC/C as function of the plate current for constant plate 
voltages of 20(A), 50(B), 80(C) volts. 
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derivation is equal to unity. Therefore Cy =C, or the original conductance 
of the affected layer is equal to the original conductance of the whole cell. 
In other words ionizations due to the impinging electrons take place 
throughout the entire material at approximately the same rate and do 
not confine themselves to a surface layer. This is in excellent agreement 
with observations by Brown,' which indicate that the change in con- 
ductance due to illumination of selenium spreads from the illuminated 
regions of the material to the unilluminated portions, sometimes as much 
as 10mm. The explanation for this behavior, i.e. the mechanism for the 
necessary transfer of energy will probably have to be sought in the 
production of secondary radiation. ° 


TABLE I 


I k 


(micro-amperes) 
100 


98.7 
200 98.0 
300 97.5 
100 99.3 
200 99.0 
300 97.7 
100 99.3 
200 98.3 
300 98.5 


The writer wishes to express his gratitude to the department of physics 
for many helpful suggestions and especially to Professor A. P. Wills 
for the constant interest taken by him in the progress of this research. 


PHOENIX PHYSICAL LABORATORIES, 
COLUMBIA UNVERSITY, 
May 26, 1924. 


' F.C Brown, Phys. Rev. 5, 404 (1914) 





DIFFUSION OF ALKALI SALT VAPORS 


COEFFICIENTS OF DIFFUSION OF CERTAIN ALKALI 
SALT VAPORS IN THE BUNSEN FLAME 


By GreorceE E. Davis 


ABSTRACT 


The method used is a modification of one suggested by H. A. Wilson. 
The ellipsoidal streak of luminous vapor from a salt bead held in the flame, was 
photographed; then the coefficient of diffusion K was obtained from the equa- 
tion K = 3rv(dx/dr—1), where r is the radial distance from the bead to a point 
on the boundary of the streak, x is the vertical distance of the point above the 
bead, and v is the velocity of the flame gases. A Meeker gas burner of special 
design was used, supplied with air under pressure and with intermittent puffs 
of salt spray to enable v to be determined. The flame temperature, determined 
by Féry’s method, was 1436°C. The means of about ten determinations of K for 
each salt are: lithium chloride, 12.4; lithium sulphate, 9.7; sodium chloride, 
19.0; sodium sulphate, 16.4; potassium sulphate, 11.1; rubidium sulphate, 11.7; 
caesium sulphate, 8.4. These are probably accurate to about 10 per cent. 
The chlorides of both Liand Na diffuse faster than the sulphates. Except in the 
case of the lithium salts, for which the observed values are relatively too low, 
the diffusion coefficients increase roughly as the reciprocals of the square roots 
of the weights of the metal atoms, also roughly as the reciprocals of the square 
roots of the molecular weights. If the free metal atom of the salt is the source 
of luminosity, as some other experimenters have concluded, it is probable that 
it is free only a fraction of the time, the fraction varying with the salt and the 
alkali. 


N 1912, H. A. Wilson published the results of some measurements of 
the coefficients of diffusion of alkali salt vapors in the Bunsen flame.' 
Some time later, while working under the direction of Prof. Wilson, the 
author made some new measurements of the coefficients of diffusion of a 
number of salts, including certain salts of the alkali metals. In these 
determinations, the data obtained for any one salt were not enough to 


make the result very reliable. More recently a new series of measure- 
ments has been carried out. A photographic method suggested by Prof. 
Wilson is used, with certain modifications which simplify the calculations 
and increase the accuracy of the results. 


THEORY 


The streak of luminous vapor extending upward from a bead of salt 
held in the non-luminous Bunsen flame appears to the eye to have a 
rather well defined boundary. This apparent boundary is the locus of 


1 TT A. Wilson, Phil. Mag , July 1912 
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points at which the concentration of salt vapor is equal to the minimum 
value for visible luminosity. Therefore the concentration is the same 
at all points on this boundary. In the paper published by Wilson in 
1912,! it is shown that the concentration C at any given point at a dis- 
tance r from the salt bead is given by 


= q (v/2K) (z—r) 
C= e€ 
4nrKr 


, (1) 


where g is the amount of salt evaporating in unit time, K is the coefficient 
of diffusion, v is the velocity of the flame gases, and x is the distance to 
the given point measured along the axis of the flame. From the graph of 
this equation we find that any surface of constant concentration resembles 
an ellipsoid of revolution with its major axis lying along the axis of the 
flame. Hence the apparent boundary of the luminous streak of salt vapor 
should be such a surface. Expressing r in terms of x and y, in Eq. (1), 
differentiating with respect to x, and setting dy/dx equal to zero, gives 
UF m 
=—(rm—Xm) , (2) 
Xm 

where 7, and x,» are the values of r and x where the horizontal diameter 
of the streak is maximum. Wilson observed x,, and the maximum dia- 
meter of the streak directly and computed K by Eq. (2). The principal 
objection to this method is that the maximum diameter of the luminous 
streak is very difficult to observe and locate definitely by direct observa- 
tion. 

In the experiments performed under the direction of Prof. Wilson, 
referred to above, a different equation was used. From Eq. (1) we may 
obtain 


- v(%2— X41 —retr1) 
2 log (re/r1) 


(3) 





where the subscripts refer to any two given points on the apparent 
boundary of the luminous streak. The values of x and r were obtained by 
photographing the luminous streak and taking measurements on the 
photographic plate. An average value of K was found from the results 
for several pairs of points. 

In the present measurements a still different equation was used. 
Differentiating Eq. (1) with respect to r, we get 


rv {dx 
K=— —-1) : 
2 \dr 
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which holds for points on the apparent boundary of the luminous streak. 
Values of x and r were obtained from the photographic plate, as before. 
These values were plotted and the slope dx /dr obtained from the resulting 
curve for a number of points corresponding to a chosen series of values of 
r. For such a series of points, we obtain from Eq. (4) 


v dx 
cram (C).- 
2(1/r)av OOF «xn. 


Comparing this with Eq. (4), it is seen that for the average slope there is 
a corresponding value of r, say 7, equal to the reciprocal of the average 
value of 1/r and independent of K and v. 7 and the average value of 
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Fic. 1. Gas burner 


dx/dr were substituted in Eq. (4) for r and dx/dr, respectively. This 
method of averaging values makes only one substitution in the equation 
necessary. The same value of 7 can be used for a number of different 
curves. 


EXPERIMENTAL METHOD 


In order to obtain a large gas flame, a Meeker type of burner of special 
construction was used (Fig. 1). A large number of holes were drilled close 
together in an iron disk D. This formed the burner top. B is a vertical 
brass tube, joined to a smaller horizontal tube C. The chamber E con- 
sists of an iron tee connected with the air intake tube A by means of 
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reducing fittings. Compressed air was used. Tube A was arranged as 
shown in order that the air pressure might not force the gas back. The 
long length of tubes B and C allowed the gas and air to mix thoroughly 
before being burned. Shorter lengths, with baffle plates as mixers, were 
tried, but were found impractical because of the tendency of the flame 
to ‘strike back.”” A tube (R, Fig. 3, not shown in Fig. 1) 5 in. long, 
diameter 2.5 in., was supported so as to be coaxial with tube B and extend 
about .75 in. above the burner top. This tube helped to steady the flame. 
The flame appeared to be perfectly uniform and did not taper much for a 
distance of several centimeters above the burner. The small tube T was 
used to introduce puffs of salt spray into the flame in order to measure 
the velocity of the flame gases. 

To obtain these puffs of salt spray, a special type of valve (Fig. 2) was 
designed. C is a cylindrical chamber made by drilling out a solid brass 
cylinder. D is an iron disk screwed onto the end of a shaft S fitted with 
sleeve and set screws (not shown) for connecting with the shaft of an 
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Fic. 2. Valve for giving puffs of salt spray 


electric motor. Face f of disk D is machined smooth and made accurately 
perpendicular to the axis of S. The inner face f of the cylinder end also 
is turned smooth and is held against D by means of spring P which carries 
a small pivot V. This arrangement allows S and D to turn freely without 
imparting much torque to C. Near the circumference of D are drilled 
8 small, evenly spaced holes h, so placed as to come in turn into momen- 
tary alignment with the hole through nipple B, as D rotates. The opposite 
end of chamber C is closed by a brass disk K which carries a nut N and a 
nipple A, all machined from the same piece. When in use, the barrel C 
is supported by a rubber band placed around it, which keeps C from 
turning yet takes up slight movements due to lack of perfect alignment of 
shaft S with the motor shaft. Salt spray entering through A passes out 
at B in a series of puffs as disk D rotates. These puffs are led into the 
flame through tube 7 of Fig. 1. 

The general arrangement of apparatus is shown in Fig. 3. W is a plat- 
inum wire with a very small loop at the end, holding a salt bead B in the 
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center of lame F as seen from above. The section of W heated by the 
flame was adjusted so that only the salt bead showed on the photographic 
plate when an exposure was taken with camera C. S is a vertical steel 
millimeter scale with its edge nearest B located in the plane which passes 
through B perpendicular to the axis of the camera lens. A Bausch and 
Lomb achromatic lens was used, diameter of full aperture 1 in. It gave 
an image free from distortion, using full aperture. Panchromatic plates 
were used. Adjustment was made for nearly equal size of object and 
image. 


70 MERCURY 
MANOMETER 


Fic. 3. General arrangement of apparatus 


Connections were made with the compressed air line as shown. A 
tank 7 of about 12 cu. ft capacity helped to steady the flow. The pressure 
of the air going to the burner was regulated by a special type of escape 
valve, the essential parts of which are shown in Fig. 4. This consisted of a 
brass tube A working easily but snugly inside a second tube B. Melted 
lead poured into A closed the upper end and furnished part of the neces- 
sary weight. Adjustment of weight was effected by adding lead shot. 
Vaseline was used to lubricate the movable tube and prevent the escape 
of air between the tubes. Tube 4A had a number of very small holes 
drilled through it in two vertical rows located opposite each other. The 


action of the valve as a pressure regulator is obvious. It kept the pressure 
constant within a range of 1 or 2 cm of mercury. It was supplemented by 
hand regulation when necessary, using stopcock S; and pinchcock ); 
(Fig. 3). As a rule, the air pressure was adjusted to about 30 cm of 
mercury. Then the gas pressure was adjusted until the inner cones of the 
flame were just below the height at which they flared up. These condi- 
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tions were kept as nearly constant as possible in all the determinations. 
Mercury U-tube manometers were connected as shown. A tilted U-tube 
with water in it indicated the pressure in the gas line (not shown). 
The velocity of the flame gases was determined as follows. Compressed 
air was let into atomizer X by opening stopcock S;. Salt spray, from a 
solution of sodium chloride in X, passed partly to valve V (previously 
described) and partly through a tube leading out of the room. The shaft 
of V was connected directly to the shaft of an electric motor M. The 
motor shaft also carried a 10 in. metal disk D with 8 evenly spaced radial 
slits cut in it. When D and the disk of valve V were rotated, 8 puffs of 
salt spray were introduced into flame F per revolution. Viewed through 
the 8 slits in D, conveniently by means of mirror m, these puffs appeared 
to be stationary in the flame. The distance between puffs was measured 
by use of a cathetometer with telescope ¢t. H is a screen for ¢t. Settings on 
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Fic. 4. Constant pressure air valve Fic. 5 


the puffs were made where their diameters were largest. Three puffs were 
visible at one time. Multiplying the speed of the motor in r.p.s., deter- 
mined with a revolution counter, by 8 times the distance between puffs, 
gave the velocity of the flame gases. 

The temperature of the flame was found by Féry’s method,’ using 
the carbon filament of an incandescent lamp as the body whose tempera- 
ture is adjusted to that of the flame. Adjustment was made for reversal 
of the sodium D-line. The filament temperature was measured by means 
of a Leeds and Northrup optical pyrometer. As a check upon results, 
readings also were taken on the flame from an ordinary Bunsen burner. 


* Féry, Comptes Rendus, 137, 909-912 (Nov. 30, 1903) 
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These latter readings fell above the scale limit of the pyrometer but could 
be estimated without much error. 

All the experiments were performed in a darkroom. The waste flame 
gases and salt vapor were led out of the room by a flue made of 6 in. stove 
pipe, supported with the lower end directly over the flame. 

Several photographic plates (on the average eight) were exposed during 
each run, without allowing the flame adjustment to change. Usually two 
or more different salts were represented in each such group of exposures. 
This made possible a more accurate comparison of the coefficients of 
diffusion of the various salt vapors. 


6 


4 
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Fig. 6. Coefficient of diffusion of salt vapor as function of reciprocal of square root of 
weight of metal atom (I), and of molecular weight (II). 


The method of marking the plates was as follows. Referring to Fig. 5, 
the center of the image O of the salt bead was first marked with a needle 
prick. Then a line ab was cut across the image A of the luminous streak 
as far from O as the distinctness of the apparent boundary EE of the 
streak would allow. This line was cut as nearly perpendicular to the 
axis xx’ as could be judged before the axis was located. The points e and f 
where the line met the boundary EE were then marked, ef carefully 
bisected, and a line xx’ cut through the middle point k and the center of 
the bead. Then the plate was put on the table of a dividing engine and 
from 7 to 15 evenly spaced lines like cd cut 2 or 3 mm apart, perpendicular 
to xx’, above O. The points where these lines met the boundary EE were 
then marked. In many cases the image A faded out quite gradually at 
the edge, so that it was somewhat difficult to locate these points. In 
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locating them, some consideration was given to the fact that the theory 
of diffusion indicates that the points should fall on a smooth curve on 
either side. 

A different method of locating points of equal intensity of blackening 
was tried, using a microphotometer of the form described by Meggers and 
Foote,’ and considering points near the apparent boundary EE. This 
method seemed to promise more accurate results than could be obtained 
by the method of marking just described. However, the results actually 
obtained were often inconsistent and unsatisfactory, so the method was 
not used. 

To measure the plates, they were placed on a small table carried by 
the movable table of a dividing engine and viewed by means of a fixed 
microscope mounted above. Distances could be read accurately to .001 
cm. The distance (as gh, Fig. 5) from edge to edge of the image of the 
luminous streak was measured along each of the cross-lines and the 
corresponding radial distance r calculated. The distance representing 
unit length on the plate was obtained by taking a series of readings on the 
image of a scale appearing on the plate along the image of the luminous 
streak. Incidentally it was found that interposing a sheet of reddish- 
brown glass between illuminator and plate made the images on the plate 
stand out much more clearly than when illuminated by white light. 

Next, x was plotted against 7 on accurately ruled paper of large size. 
The curves were from 45 to 80 cm in length. They were very carefully 
drawn and the plotted points usually fell on, or very close to, a smooth 
curve. The coefficient of diffusion of the salt vapor was obtained from the 
’ The values of the slope 
dx /dr were found by laying a straight-edge tangent to the curve at the 


’ 


curve by the method described under ‘“Theory. 


chosen points. This could be done with a fair degree of accuracy, since 
the curves were very nearly straight in every case. 

It is assumed in the mathematical theory that all the stream lines of 
the moving medium are parallel. Therefore it is important to determine 
whether the tapering of the flame and the distortion of the stream lines 
near the salt bead have an appreciable effect upon the results. As a test, 
a group of 25 plates giving long x, ry curves were chosen at random. Two 
determinations of the diffusion coefficient K were made from each curve, 
one from the upper 5 values of r in the series already chosen, and one 
from the lower 5. Then the average K obtained from the upper points 
was compared with that from the lower. As a further test, two groups of 
4 plates each were chosen, in each group the kind of salt and the flame 


Meggers and Foote, Scientific Papers Bur. Standards, No. 385 (June 30, 1920) 
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velocity being the same for all the plates: From the average dimensions 
in each group, a curve was plotted to represent the apparent boundary 
of the image of the streak of vapor. This curve was compared with a 
theoretical curve which had arbitrarily been made to coincide with it at 
two points. The results are given below. 


EXPERIMENTAL RESULTS 


In Table I is shown a typical set of values of dx/dr obtained from one 
of the x, r curves. In Table II are given the results obtained from the 
individual photographic plates, with their averages. The results obtained 
by H. A. Wilson are shown for comparison. The particular salts used in 
his determinations are not mentioned in his article, so are not specified 
here. 

In Table III are shown the temperature readings at three different 
points in the flame, under different conditions as to pressure of air supply 
to the burner. Since the air pressure determines the velocity of the flame 
gases, the measurements show that the temperature was not appreciably 
affected by velocity variations much greater than those occuring during 
the determinations of the diffusion coefficients. The temperature is the 
same at all three points in the flame, within the limits of accuracy of the 
measurements. These points are so located as to give a fair indication of 
the temperature throughout the region in which the coefficients of diffu- 
sion were determined. Therefore we conclude that the results were not 
appreciably affected by variations in temperature from point to point 
in the flame. The average value, 1436°, may be a little lower than the 
true temperature, since the carbon filament used in the measurements 
was not a perfect black-body radiator. The true temperature may be 
from 45° to 90° higher. However, the value 1850° for a Bunsen flame of 
the ordinary type agrees well with the most reliable values found by 
other observers. 

The gas burned was manufactured from bituminous coal and may be 
described as a mixture of coal gas and “water gas.”’ Its composition 
varied somewhat from day to day, as shown by analyses made by the 
department of chemistry. These analyses showed the following content, 
in volume per cent: CO, 11 to 16; COs, 4 to 7; Oc, 1.5 to 6; unsaturated 
hydrocarbons, 6 to 17; CHa, 11 to 16.5; He, 18 to 21; No, 27 to 32. 

In the first test made to determine whether the values found for the 
coefficients of diffusion K were affected appreciably by lack of parallelism 
of the stream lines of the flame gases, the average value of K for the 


* Burgess and LeChatelier, ‘‘Measurement of High Temperatures’’; observations 
by Burgess, pp. 340, 341. 
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upper 5 values of r came out 4.8 per cent higher than the average for the 


lower 5. 


In the second test, comparing theoretical shapes of streaks of 


salt vapor with shapes experimentally determined, a small departure of 
the latter from the former was found in each case. These tests seemed 
to indicate an appreciable but rather small effect on the values found for 
K, due to lack of parallelism of the stream lines. 
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Reciprocal of average value of 1/r=1. 616. 


TABLE I 


x 
.614 
825 
.035 
243 
.450 
.656 
862 
068 
273 
477 
. 680 
. 884 
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Coefficient of diffusion computed by Eq. (4) =8. 2. 
* These values were obtained from a plotted curve, not computed from the values 


of r and x. 
LiCl Li,SO, 
12.7 12.0 
14.5 12.9 
15.0(3) 10.8 
16.6 10.2 
3.9 11.1 
10.9(4) 6.7 
9.3 8.4 
15.3 9.6 
10.8 8.5 
11.9 6. 5(4) 
7.9 8.9(3) 
Weighted 
average* 12.4 9.7 
Average 
deviation +2.2 +1.6 





Wilson's resultst 


(Li)14.5 


| ‘(Na)11.5. - 


Typical set of results 
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.013 


Average dx/dr=1.031 





Velocity of flame gases=328 cm/sec. 








TABLE II 
Coefficients of diffusion 

NaCl Na:SO, K,SO, Rb.SO, Cs,SO, 
18.8 16.3 12.9(3) 11.7 8.2 
22.6 18.7 17.1 11.1 8.5 
15.8(4) 14.8 14.2 13.7 8.7 
20.8 20.9 12.5 11.4 9.0 
22.9(4) 21.8 8.6 11.3 7.8 
19.5 15.7 9.5 13.1(3) 9.6 
14.0 14.8 8.7(4) 10.6(4) 7.5 
19.2 13.3 9.1(4) 12.0 8.0 
18.6 15.4 6.3(4) 10.3(4) 8.0 
15.8(4) 13.7 
20.4 14.6 

19.0 16.4 11.1 11.7 8.4 
+2.1 +2.3 +2.8 + .8 + .§ 





- (K)4.7. (Rb)4.7. (Cs)4.4 


* Scale of weights, 3 to 5; the weight is 5 unless otherwise indicated by figures in 
brackets; weighted on basis of excellence of photographic plate and apparent reli- 


ability of data. 


t The particular kind of alkali salt examined in each case is not stated by Wilson. 
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TABLE III 
Flame temperatures, by Féry's method 
Position in Pressure of air Temperature 
flame supply 
Lower center of 18.3 cm Hg 1435°C 
flame, just above 30.8 1440 
salt bead 31.1 1433 
.8 cm from 
boundary of flame, 30.9 1426 
just above level of 
salt bead 


Center of flame, 

about 2.5 cm above 30.9 1444 

salt bead 49.0 1440 
Av. temp. = 1436° 

Reading on Bunsen flame of ordinary type, with full air supply = 1850°. 


DISCUSSION 


The coefficients of diffusion found in these measurements (Table II) 
do not agree well with those found by Wilson. It is believed that the 
results obtained in the present measurements are the more reliable for 
reasons given under Theory. A difference in composition of the gases 
burned in the two flames, or a difference in flame temperatures, may 
partly account for the disagreements. The above results are probably 
accurate to about 10 per cent. 

It has been shown by Smithells, Dawson and H. A. Wilson® that the 
light emitted by a flame containing vapor of an alkali salt is probably 
due to atoms of the alkali metal. They pointed out that a metal atom 
might lose one or more electrons under the influence of violent collisions 
with other atoms in the flame, thus becoming a positive ion; and that 
this ion, upon recombination with other free electrons in the flame, might 
be subject to internal vibrations which would give rise to the radiation 
of light. This theory was strengthened by the results of later experiments 
by H. A. Wilson.’ These experiments showed (1) that the luminous 
vapors of salts in a Bunsen flame are not deflected appreciably by an 
electric field; (2) that the positive ions present in the luminous vapor can 
be made to move out of it by an electric field; (3) that the positive ions 
are not luminous, but can form luminous vapor after recombination. 
It has also been shown by Bancroft and Weiser’ that many metallic 
salts are dissociated in the Bunsen flame, the metal atom being set free. 


5 Smithells, Dawson and H. A. Wilson, Phil. Trans.. A 193, 89 (1899) 
6H. A. Wilson Phil. Trans., A 216, 63-90 (1915) 
7 Bancroft and Weiser, Journal of Phys. Chem., 1914. 
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If, as this evidence indicates, the metal atom of the vapor is the source 
of light, we might expect the coefficient of diffusion, determined from the 
shape of the luminous streak of vapor from a salt bead, to be the coeffi- 
cient of diffusion of the free metal atoms. But it is likely that these atoms, 
because of their tendency to combine with other ions (radicals) in the 
flame, are frée only part of the time. Therefore the average rate of 
diffusion for any particular atom probably is less than the rate we should 
expect if the atom were free during the whole of its journey. It therefore 
seems reasonable to assume that our observed values for the coefficients 
of diffusion are lower than the values for free metal atoms. 

It is interesting to examine our results in the light of two different 
assumptions: (1) That the metal atoms were free during the entire time 
spent in diffusing; (2) that the molecules of salt vapor remained undis- 
sociated practically the entire time. In the first case, kinetic theory indi- 
cates that the observed diffusion coefficients should vary directly as the 
reciprocals of the square roots of the weights of the metal atoms. In 
Fig. 6, curve I, these quantities are plotted against each other (points 
marked with circles). In the second case, the observed coefficients should 
vary as the reciprocals of the square roots of the molecular weights of the 
salts. Plotting these quantities, we obtain curve II (points located by 
crosses). In either case, the values for the sodium, potassium, rubidium 
and caesium salts locate, roughly, a straight line passing through the 
origin, as we should expect. On the other hand, lithium chloride and 
lithium sulphate do not fall near the line located by the other salts in 
either case, their observed coefficients of diffusion being only roughly 
1/3 and 1/2 the values indicated by curves I and II, respectively. The 
degree of accuracy with which each curve is located by its set of five points 
is about the same for the two cases. Again, the positions of the two 
sodium salts relative to a straight line drawn between them from the 
origin, depend on whether assumption (1) or assumption (2) is made. 
The same is true of the two lithium salts. These facts suggest that neither 
assumption is correct, but that the metal atoms are free part of the time 
and combined with other particles the rest of the time. As to the two 
lithium salts, their relatively low rates of diffusion may be due to diffusion 
of the particles in groups rather than singly. 

Partial dependence of the rate of diffusion upon the negative radical 
of the salt is shown by the observed differences between the chlorides and 
the sulphates of lithium and sodium. In each case the chloride diffuses 
more rapidly than the sulphate of the same metal. The differences seem 
too large to be attributed to experimental errors alone. Referring to 
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Table II, we see that 4 values obtained for lithium chloride are higher 
than the highest value for lithium sulphate; while 7 values for sodium 
sulphate are lower than the lowest value for sodium chloride except one. 

The various chemical and physical phenomena which occur in the flame 
are not yet well enough understood to warrant drawing more definite 
conclusions. 


=e 


It is a pleasure to the writer to acknowledge indebtedness to Professor 
L. B. Spinney for his kindness in securing the necessary equipment for 
this work; also to Professor G. E. Thompson and especially to Professor 
J. W. Woodrow, for many helpful suggestions. 
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AN ELECTROMETER METHOD FOR MEASURING DIELEC- 
TRIC CONSTANTS OF LIQUIDS 


By A. P. CARMAN 


ABSTRACT 


Differential idiostatic electrometer for the measurement of dielectric 
constants of liquids.—The suspended system includes two cylindrical quad- 
rants C; and C; suspended so as partly to overlap two fixed cylindrical plates 
A and B coaxical with C; and C, and so placed and connected that the torque 
due to the condenser AC, is opposite to that due to the condenser BC;. To 
determine the dielectric constant of a liquid, one of the condensers is im- 
mersed in it and the other immersed in a liquid or gas of known constant, and 
the differences of potentials V4—V, and V:— Vg for zero deflection are read 
by the potentiometer method. The dielectric constants are then proportional 
to (V4—V:)?/(V.— Vpg)?=Ri?/R2?. This method gave for the dielectric 
constant for distilled water at 25°C with 60 cycle e.m.f. the value 78.07. 

HE instrument used is a form of idiostatic electrometer, or electro- 

static voltmeter, part or all of which is immersed in a liquid or liquids. 
The immersed electrometer for measurement of dielectric constants was 
first used by Silow! in 1875, and similar methods have been used by 
several investigators’? since then. Two investigators, Perot® and Heer- 
wagen‘ employed differential instruments, each using a double quadrant 
electrometer with the needles on a common stem, of which the lower set 
of quadrants and needle can be immersed in the liquid. The theory of 
these two methods is similar to that presented below, but the construc- 
tion, adjustments and manipulation of the instrument are very different 
from those of the instrument described in this paper. 

The electrometer used is shown in Fig. 1. The plates C; and C; of the 
needle system are cylindrical and are suspended to rotate about the 
vertical axis. The fixed plates A and B are parts of a coaxial external 
cylinder, placed as shown in the figure to produce a torque on the needle 
when differences of electrical potential between the needle and the plates 
are maintained. The shape of the wire frame carrying the needle plates 
is shown in the figures and, as indicated, either plate C; and its opposite 
fixed plate A or C; and B can be immersed in a given liquid by raising 


1 Silow, Ann. der Phys. 156, 389 (1875) 

*Cohn and Arons, Ann. der Phys. 33, 13 (1888); Tereschin, Ann. der Phys. 36, 
792 (1889); Tomaszewsky, Ann. der Phys. 33, 33 (1888); Rosa, Phil. Mag. 31, 188 (1891) 

* Perot, Jour. de Phys. 10, 149 (1891) 

* Heerwagen, Wied. Ann. 48, 35 (1895) 
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around the plates a beaker containing the liquid. The suspension is by a 
long phosphor-bronze ribbon which also serves for charging the needle to 
the desired potential. The needle is damped by means of an attached 
vane, which dips into a beaker of oil. The deflections of the needle system 
are read by the standard method of mirror, lamp and scale. 

In the actual instrument which was made and used, the plates A, B, C, 
and C2, were each 2.5 cm square and the needle plates were 6.5 cm apart. 
The distance between the needle and fixed plate and also their relative 


’ 



























































Fig. 1A 
Fig. 1A. The electrometer. Fig. 1B. Electrical connections. 
Fig. 1C. Section through the plates. 


positions can be varied. Other sizes can of course be used and a design is 
possible in which the pairs of plates are not the same in size on the two 
sides. Gold plates have been used for liquids which corrode other metals. 
For some dense liquids, the boyancy is compensated by small lead weights 
on the horizontal bar of the needle. 

The electrical connections are shown in Fig. 1B. The plate A is brought 
to the potential V;, the needle plates C,; and C, to V2, and the plate B 
to Vs, the potential V, being higher than V2, and V2 being higher than 
V;. This is conveniently done by the usual potentiometer device of tap- 
ping off connections from a uniform wire which is stretched zig-zag on a 
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board of convenient length. In many of our experiments, the potential 
difference Vi— V3 was 110 volts, and for conducting liquids this was 
alternating. 

The plates are placed so that the charges on the pairs of plates tend to 
rotate the needle in opposite directions. Now for equal arm lengths the 
moments of force are respectively proportional to the squares of the 
difference of potential (Vi— V2)? and (Ve—V3)?, and these in turn are 
proportional to R,? and R,”, where R; and R,» are the resistances of the 
two parts of the potentiometer wire. But the forces are also proportional 
to the dielectric constants K, and Kz of the media between the pairs of 
plates A,C,; and B,C.. Hence for the null position of the needle, we have 


Ky RY 
— Cc —- ° 
Re RR? 


To determine the calibration constant C, the resistances Rio and Reo 
for the zero balance when air is the medium between the plates on both 
sides are found; then if the resistances are R’; and R’s when the first 
pair of plates has the medium air and the second pair of plates has the 
given liquid X as the medium, assuming the dielectric constant of air to 


be unity, we get 
C = (Rao/ Rio)? ; Kx = (R2’/R,’)? XK (Rio/ R20)? - 


By immersing the first pair of plates in a medium Y, and the second 
pair of plates in a medium X, we can calculate the ratio Ky/Ky of the 
dielectric constants of the two media from the ratios of the resistances in 
the calibrating and the comparison balances. This use of a chosen liquid 
instead of air for reference, has at times decided advantages. Thus in 
measurements with a substance such as water having a high dielectric 
constant, we can avoid a large and inconvenient ratio of Re and R, by 
using a reference substance with a fairly high dielectric constant, such as 
acetone or alcohol. The use of a liquid of high dielectric strength for 
reference instead of air also allows the convenient use of higher potential 
differences when the consequent increase of sensitivity is desired. 

The above instrument is simple in construction, with few parts; its 
manipulation and adjustments are easily made, and the plates are readily 
accessible for cleaning, a very important feature in working with most 
liquids. As the plates on either side can be varied in position indepen- 
dently, the sensitivity and range are easily controlled. As the needle 
is brought back to zero each time, the actual readings used are readings 
of resistances which can be made with great accuracy. The immersion 
method has the advantage of allowing the use of alternating electromotive 
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forces of a wide range of frequencies. For alternating currents the resist- 
ance readings must of course be corrected for reactance. The dielectric 
constants of semi-conducting liquids can also be satisfactorily measured 
by this method, particularly if leakage currents are measured and correc- 
tions are made for the potential differences caused by leakage. 

The above method has been used for determining the dielectric con- 
stant of water. The water was doubly distilled in pyrex glass vessels, 
and the first third and the last third of the distillate were discarded, so 
that the water used was probably fairly pure. The value of the dielectric 
constant obtained for this water was 78.07 at 25°C for an alternating 
e.m.f. of 60 cycles, a value in satisfactory agreement with accepted values. 
Most of the work which we have been doing with this method has been 
on the dielectric constants of mixtures of liquids and of emulsions and 
is not yet ready for publication. , 

In developing the above method and instrument, the author ac- 
knowledges the great help received from Mr. F. E. Null, graduate scholar, 
who has carried out the experimental work. 


LABORATORY OF PHysIcs, 
UNIVERSITY OF ILLINOIS 
June 20, 1924. 
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THE ELECTRIC MOMENT OF GASEOUS MOLECULES OF 
HALOGEN HYDRIDES 


By C. T. ZAHN 


ABSTRACT 


Dielectric constant of the hydrogen halides HCl, HBr and HI to 300°C.— 
An improved heterodyne null method of measurement was used in which 
the beats between two high frequency oscillations from separate electron tube 
generators are adjusted to the frequency of a tuning fork, the small change of 
capacity due to the introduction of the gas being compensated by a large change 
in a large capacity in series with the gas condenser. Measurements were made 
at atmospheric pressure, using carefully purified gases, over a range of nearly 
400° from just above the liquefaction point. The results are well represented 
by the Debye equation (e—1)vT=AT+B, where v is the specific volume and 
T the absolute temperature. [The values found for A are .001040 (HCl), 
.001212 (HBr) and .001856 (HI); and for B are .895(HC1), .52(HBr) and .12 
(HI)]. This agreement lends support to the Debye classical theory of electric 
polarization due to fixed moments in the molecule. The values obtained for the 
electric moment of the molecules are, in 107 c.g.s.u., 1.034(HC1), 0.79(HBr) 
and 0.38(HI). The Pauli quantum theory gives an equation of the same 
form as Debye’s but leads to smaller values for the electric moment. The upper 
limit for the moment given by infra-red absorption data for HCI, however, is 
6 times the classical value and 13 times the quantum value and hence does 
not decide between the two theories. 

Dielectric constant for hydrogen, oxygen, nitrogen and air.—The values 
obtained for (e— 1) X10 at 0°C and 760 mm pressure are 265, 518, 581 and 572 
for Ha, Oz, Nz and air, accurate probably to within 1/2 per cent. 


INTRODUCTION 


HE primary object of this investigation was to determine the electric 

moment of the molecules of HCl and HBr in order to compare their 
effective dipole length with their nuclear separation as obtained by 
Kratzer! from infra-red absorption spectra. This comparison would give, 
on certain assumptions, an idea of the distortion of the electronic con- 
figuration of the molecules. Further it was hoped that it might be pos- 
sible by the use of high frequency technique to develop a method of 
measurement of dielectric constants more accurate than had hitherto 
been attainable by the use of electrometer and bridge methods and high 
frequency methods of different character. Similar methods have been 


1 Kratzer, Zeit. f. Phys. 3, 389 (1920) 
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used by Whiddington? for the measurement of small displacements, by 
Herweg’ for studies of the effect of field strength on the dielectric con- 
stant, by Belz‘ for the measurement of the magnetic susceptibilities of 
salts, and by others. An experimental value of the dipole length also 
offered a possibility of choosing between the Debye theory and a more 
recent theory by Pauli applicable to dipole gases and based on Sommer- 
feld’s conditions of quantum mechanics. 
Debye’s equation is as follows: 


=the io —-— = lV — ’ 
+2 3 3RT wil. 


where ¢ is the dielectric constant, N the number of molecules per cubic 
centimeter, \ the constant of elastic binding of the ‘‘optical’’ electrons, 
e the electronic charge, uw the electric moment of the molecule, R the 
Boltzmann gas constant and T the absolute temperature. The Debye 
equation has been found satisfactorily to represent data obtained by 
Jona and by Badeker on H2O, NH3, SOs, CHs;OH, C2H;OH, CSe, CO» 
and CO.®*8 Both Badeker and Jona worked at temperatures in the range 
from room temperature to about 180° C. From the Debye equation it 
will be seen that the greatest contribution of the doublets, and so the 
greatest anomaly in the behavior of the dielectric constant is to be ex- 
pected at the lowest temperatures. A more accurate verification of the 
Debye equation becomes of more particular interest also in view of a 
theory proposed more recently by Lundblad,’ who considers the mutual 
effect of electrons of the same molecule. 


EXPERIMENTAL METHOD. 


The measurement of the dielectric constant of gases is rendered extra- 
ordinarily difficult by the calibration of the apparatus. In fact, a number 
of investigators, some of whose results are quoted in tables of constants, 
have relied upon previous measurements for this calibration. Relative 
values of (e—1) can be obtained much more accurately than absolute 
values. Indeed, different investigators have differed by as much as ten 
per cent in the values of gases such as nitrogen, air and hydrogen, but 


?R Whiddington, Phil. Mag. 40, 634 (1920) 
*J Herweg, Zeit. f. Phys. 3, 36 (1920) 

*V.H Belz, Phil. Mag. 44, 479 (1922) 

5 P, Debye, Phys. Zeit. 13, 97 (1912) 

* K. Badeker, Zeit. f. Phys. Chem. 36, 305 (1901) 
7M. Jona, Phys. Zeit. 20, 14 (1919) 

* J] J. Thomson, Phil. Mag. 27, 757 (1914) 

* R. Lundblad, Zeit. f. Phys. 5, 349 (1921). 
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since the values for the different gases obtained by two observers usually 
bear a nearly constant ratio, the variations in the absolute values can 
be attributed to errors in the calibration of the apparatus. It is believed 
that the method herein described has reduced this difficulty of calibration 
of the apparatus and has accordingly given results less liable to such 
errors. 

(a) The measuring circuits. The method of measurement adopted in- 
volves the use of high frequency oscillations produced by electron tube 
generators. A diagram of the electric circuits is shown in Fig. 1. There 
are two separate generators; one operates steadily at a frequency of the 


+7 





Fig. 1. Diagram of electric circuits. 


order of a million per second; the other, at a frequency differing from this 
by an audible frequency of a thousand per second. These two generators 
are coupled electrostatically to a detector-amplifier circuit consisting of 
three stages as shown. A heterodyne note of frequency one thousand 
per second is heard from the telephone receiver 7. This note is kept 
adjusted by the method of beats with a tuning fork operated by an auxil- 
iary electron tube circuit. The variable generator contains a system of 
condensers one of which can be evacuated or filled with the gas to be 
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studied. The procedure is as follows. With the gas condenser evacuated 
the variable generator is adjusted to the condition of zero beats between 
the notes from the tuning fork and the receiver 7; then the gas is intro- 
duced into the condenser and a measurable compensation is made for 
the resulting change in capacity by means of a device to be explained later. 

The various parts of the electrical apparatus are carefully shielded 
wherever necessary. The two generators and the detector-amplifier are 
separately enclosed in metal boxes. In order to avoid coupling between 
the two generators themselves, and thus to make the frequency of the 
generators independent of each other, the coupling to the amplifier is 
reduced to a minimum. It is found quite sufficient to use a ‘floating 
grid’”’; i. e. a wire connected to the grid of the detector tube and having 
two branches, one passing through a small hole in the shield of each 
generator. This wire receives sufficient charge, mainly by electric in- 
duction, to give a conveniently audible note in the receiver 7. The re- 
ceiver itself is completely shielded electrostatically, the sound being heard 
at the end of a metal tube enclosing it. All the variable condensers are 
operated by shafts running through the shielding boxes, and readings 
are made through small holes in the boxes. 

All five electron tubes are Western Electric 203B tubes operated at 
45 volts on the plate and about 3 volts on the filament. The four induc- 
tance coils of the two generators consist of about 20 turns of 10 cm radius. 
For coupling between the stages of the amplifier, ordinary telephone 
transformers are used. The plate batteries of the generators are shunted 
with telephone condensers in order to reduce the resistance of the circuits 
for high frequency currents. All connecting wires are made of solid brass 
or copper rod. 

(b) The compensation device. The system of condensers used in the 
variable generator consists of the three condensers K, K’, and C shown 
in Fig. 1. K’ isa fixed condenser of about 10000uuf and K is a variable 
precision condenser of about 1500uuf maximum capacity. The precision 
condenser has a slow adjustment screw which permits readings of capacity 
to be taken to within (1/20)uyf. The capacity (K+K’) is connected in 
series with the gas condenser C. Compensations for changes in C by the 
variation of K can be measured accurately since it will be seen that they 
can be made as large as desired according to the value of the shunted 
capacity K’. It would, of course, be quite impossible to measure the 
small capacity changes due to the influx of the gas into the condenser 


by means of an ordinary condenser in parallel with C. The series arrange- 
ment of the condensers magnifies the effect to be measured. Since by the 
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use of zero beats as a criterion of adjustment the total capacity of this 
combination is always kept constant it follows that for a compensation 
AK corresponding to a change AC produced by the influx of gas 

1 1 1 1 


KiK’'C KEK'4AK CHA’ 








C? 
ac=—4 bax 
(K+ K’+C)AK+(K+K’)? 


However, for the purpose of simplifying the calculations, a much more 
convenient form can be obtained for AC by the use of an approximate 
equation, which gives the same result as the latter equation toan accuracy 
within the limits of experimentalerror. For very small changes one could 
use the differential equation dK/dC=d(K+K’)/dC=—(K+K'/C)’, 
which, however, is not a sufficiently close approximation for the require- 
ments of this problem. The next order of approximation consists in as- 
suming that the differential coefficient dK/dC is a linear function of K. 
Then an average value of dK/dC can be used and for finite values of 
AC and AK 


9 


co-{ ve 
K+34K+K’ 


For cases where measurements were to extend overa long period of 
time a special arrangement was used to eliminate errors due to gradual 
changes in the tube circuits external to the condenser system. This 
consists in having a standard condenser so arranged that it can be inter- 
changed for the condenser system shown in Fig. 1. The interchange is 
made by a switch of special construction consisting of three equidistant 
mercury cups. Into one of these is inserted a vertical wire free to turn 
about the vertical axis and having attached a hozirontal wire which makes 
contact in turn with the mercury meniscus of either of the other two cups 
In order to compensate for possible changes in the circuits with time the 
distributed capacity of one of the generator circuits is varied by means 
of a threaded brass rod with a sinail brass disk attached to it on the in- 
side of the shield. Small changes in frequency can be obtained by turning 
this rod and thus varying the distance of the small disk from the generator 
circuit. 

(c) The gas condenser. The condenser C is a cylindrical condenser of 
about 4 cm diameter, 13 cm length and 0.5 mm separation between cyl- 
inders. Fig. 2 shows the main features. The two cylinders are held to- 
gether and insulated by two quartz disks Q, which give a minimum of 
insulating dielectric to be corrected for in the calibration of the apparatus. 
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In order to prevent chemical action the cylinders themselves were made 
of an alloy of gold with small quantities of palladium and platinum. 
A gold plated brass condenser was first tried, but the distillation of zinc 

contaminated the gold plate at the higher temperatures. Around the 

outside cylinder is wound a platinum resistance thermometer P, insulated 

by a sheet of mica M wrapped around the condenser. The condenser is 

mounted in a Pyrex tube and held in position by tungsten wires fastened 

at the top and bottom of the gold cylinder. (These are not shown.) The 

two gold leads LZ and the two platinum leads P are 

first welded to tungsten wires T and then the seals 
in the glass tube are made at the top with the 

tungsten wires. A capillary tube is led from the 
tube to the vacuum and gas system. The Pyrex 

tube is then mounted in a small shielding cylinder 
of brass (plated with nickel to prevent corrosion 

by liquid air). This shielded condenser is then 
mounted in a large shielded box to permit the 
application of cooling baths and of the electric 

furnace without introducing capacity effects with 

the rest of the apparatus. 

‘The resistance thermometer consists of about 
15 ohms of platinum wire wound around the con- 
denser as loosely as is compatible with its staying 
in position. By means of a Wheatstone bridge the 
resistance can be measured to an accuracy of 0.001 
ohm which corresponds to 0.02°C. The ther- 
mometer was calibrated under the experimental 
conditions to be used by means of three fixed 
points; the freezing and boiling points of water and 
the equilibrium point of a mixture of carbon dioxide 
snow and ether. In consideration of possible errors 
due to the calibration it is believed that this 
Fig. 2. Gas condenser. thermometer records temperatures accurate to at 

least 0.1 degree. 

(d) The calibration of the condenser system. It will be noted that this is 
a null method and it depends in no way upon either the electron tube 
constants or the distributed capacity of the coils, etc., provided these do 
not change erratically. Actually, the limit of accuracy of the method 
depends upon the calibration of the condenser system. To make this as 
accurate as possible, the condenser system was so constructed that the 
whole calibration could be made in terms of one single condenser, the 
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precision condenser K. The latter was calibrated by the Bureau of 
Standards in terms of the standard air condensers of the Radio Section. 
The capacity K’ consists of three fixed mica condensers arranged so that 
they can be inserted in parallel with K by means of mercury cup con-. 
tacts. These capacities K’ were such values that they could be calibrated 
in terms of K by the method of substitution. This has the advantage of 
preventing errors which might result from the dependence of the capacity 
of the fixed mica condensers upon the frequency. Here they are calibrated 
at precisely the frequency at which they are to be used. The capacity 
C is determined by shortcircuiting K and K’ by amalgamated contacts 
and tuning the fixed generator to zero beats. Then C is shortcircuited 
and K is substituted and tuned thus giving a measure of C. With this 
arrangement an error of the unit in the calibration of K would produce 
no error in e so long as relative values of K are correct. To check the 
calibration of K a small fixed mica condenser was constructed and dif- 
ferences on K were checked with satisfactory agreement, so that settings 
on K are good to at least 3yyuf. 

In the calibration of the condenser system the greatest possibility of 
error resides in the determination of the portion of the capacity C which 
is not altered by the influx of the gas. This is in part due to the quartz 
disks and in part to the capacity with ground of the lead passing from C 
to (K+XK’). The inductance of this lead, even if appreciable, behaves 
just as though it were outside the system and in series with it and does 
not come into consideration. This lead will possibly have a small specific 
inductive capacity relative to the shield KK’ and also with the grounded 
shield of C and other grounded apparatus such as the large shielding 
box enclosing the whole generator circuit. The former capacity is negli- 
gible as determined by both experiment and approximate calculation. 
The latter capacity is not negligible; it consists chiefly of the capacity 
of two wires about 15 cm long and 2 cm apart. These capacities were 
determined by a rather laborious process to be described below, by con- 
structing a set of leads to imitate the actual leads to be used on the 
condenser tube and making measurements on these placed in position 
within the shields mentioned above. The lead capacity was thus found 
to be 5.5upuf. To obtain the alterable capacity this must be subtracted 
from the total capacity C = 204.5uyf, and in addition there must be sub- 
tracted the capacity due to the portion of the dielectric space occupied 
by the quartz disks. A calculation for this latter capacity gives 0.5yupf. 
Hence the final value for the alterable capacity, which may be; called 
Co = 204.5~—5.5—.5=198.5uuf. The total correction for the lead and the 
quartz disk is less than three per cent of the total capacity C. A fair 
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estimate of the greatest possible error in this determination would be 
about 1pyf or 1/2 per cent. 
The final equation for the dielectric constant is then: 
xe C AK - AK 
6 — ; 
[K+34K+K’}?* 








TOG [K+HAK+R'}? 

The method of determining the lead capacity will now be considered 
in detail. The condensers K’ were disconnected and condenser C was 
not yet mounted in its position in the shield S in Fig. 3a. L represents 
the immitation lead, at the end of which is connected by an amalgamated 
contact A another wire which passes out through the opening made by 
taking off the bottom of the shield S. The contact A is made or broken 
by a sidewise displacement of the lower wire. G is the ground capacity 
of the shield of K; and x and y, the lead capacities mentioned. C; is 
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(a) (b) 


Fig. 3. Connections used in determining lead capacity. 


another precision condenser calibrated in terms of K ; and 7, the terminals 
of the coil of a high frequency generator which produces beats with the 
fixed generator. In Fig. 3b this condenser group is redrawn symbolically 
for the sake of clarity. Now if the circuit is broken at A and a compensa- 
tion A,C is made in C with the capacity K set at a particular value 
AiC=x+G(y+K)/(G+y+K) 

where x, y, and G are unknown and A,C and K are obtained from con- 
denser readings. By taking values of AC for three settings of K there are 
obtained sufficient equations to solve for the three unknowns. A simpler 
and more accurate method is to use one equation of the latter type and 
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two other equations corresponding to the compensations A:C and A,C 
when G and (K+) are shortcircuited respectively. Then the following 
three much simpler equations are obtained 

A,C=x+G(K+y)/(G+K+y) 

A.C=x+K+y 

A;C=x+G. 
The values of ‘K, x, and y were obtained in this way and the method was 
checked by using different settings of K. It was in that way that y was 
found to be negligible. 

(e) The vacuum and gas system. By the use of a diffusion pump the gas 
condenser could be evacuated within three minutes to a pressure as low 
as 10-*mm Hg, which is far beyond the requirements of these experi- 
ments. The use of such a good pump is however important in ridding the 
system of water vapor, which is particularly necessary in the study of 
acids because of catalytic action. Before making a series of measure- 
ments the tube was first baked for several hours at a temperature of 300°C 
and evacuated to a pressure of 10° mm Hg or better. By taking this 
precaution and using P.O; and a liquid air trap no difficulties were en- 
countered such as prevented Bideker from making a satisfactory study 
of HCl gas. Special arrangements were made for the generation and 
purification of each gas. These will be described in connection with the 
data on the individual gases. 


MEASUREMENTS ON AIR, NITROGEN, HYDROGEN, AND OXYGEN 


In order to test the behavior of the circuits, and to obtain data for 
comparison with those of other observers, measurements were made at 
room temperature on air, nitrogen, hydrogen, and oxygen. The air was 
passed through a P.O; tube to absorb water vapor and through a liquid 
air trap to freeze out the COs. The nitrogen was prepared by heating 
sodium nitrite and ammonium chloride solutions together and freed from 
oxygen by storing it in a reservoir containing yellow phosphorous. It 
was then passed through a P2O; tube with a liquid air trap. The same 
process of purification was used for hydrogen, which was prepared 
electrolytically. The oxygen also was prepared electrolytically and then 
passed through a P.O, tube for drying. 

A considerable number of readings was taken for these gases as shown 
in Table I. For all these individual readings the deviation from the mean 
is not greater than 1/2 per cent of the value of (e—1) and the probable 
error of the observations is of the order of 1/10 per cent. The average 
value of the series capacity (K+ XK’) is 10000uuf which gives a magnifica- 
tion of 2500 in AC. The measured capacity AK is about 200uuf. All of 
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the measurements were taken at room temperature and approximately 
room pressure. Then the values were reduced to 0°C and 760 mm Hg 
by the use of the perfect gaslaw pu= RT. This latter procedure is justified 
by the fact that these gases are known to show no anomalous behavior 
of the dielectric constant. 


TABLE I 
Results for air, Nz, Hz and O; 








No. of (e—1) 
Gas readings reduced 





39 .000572 
41 .000581 
19 .000265 

9 .000518 


«4 








It is intended to give in the future a discussion of these data as related 
to the values of other investigators and to data on the refractive index 
in view of the Maxwellian relation. Suffice it to say here that these 
values agree with those of other investigators as well as the latter agree 
among themselves.'® It is, however, believed that, in consideration of 
the improved method previously described, these values are reliable to 
within one per cent. 


MEASUREMENTS ON HCl, HBr, Anp HI 


After the preliminary measurements had been completed a series of 
observations were made on the gases HCI, HBr, and HI. The first method 
of procedure tried with HCI was unsatisfactory. The gas was let in after 
cooling the apparatus almost to the liquefaction point of HCl. Then 
readings were taken on the condenser K as the apparatus approached 
room temperature. This required several hours and it was found after- 


1 In anticipation of later comment, it seems desirable to call attention to a recent 
paper by Fritts (Phys. Rev. 23, 345, 1924) which gives data on the dielectric constants 
of oxygen, nitrogen and air differing by several per cent from those of the present 
paper, although claiming about the same accuracy. Fritts used a very ingenious method 
of recording the beats, so that his actual observations are certainly very accurate. 
His precautions regarding purity of gas and calibration of apparatus do not seem to be 
equally trustworthy, as judged by experience during the present investigation. We 
found the step-by-step calibration, in addition to the Bureau of Standards calibration, 
necessary in order to realize the possibilities of the method. We also found the method 
so sensitive as to detect impurities left in the gases after ordinary chemical 
methods of purification, so that unusual methods were employed to advantage. The 
correction for lead capacity is also important and difficult to be sure of. Since the 
present apparatus was designed primarily to permit convenient temperature adjust- 
ments, a new apparatus is being used to investigate further the sources of error in ab- 
solute values. For this reason further discussion of absolute values is deferred until a 
later communication. 
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ward that the gas had become contaminated, either as a result of chemical 
action or by gases occluded in the condenser tube. It was therefore found 
necessary to use another method of procedure, wherein after each gas 
reading a vacuum reading was taken and impurities in the gas were 
removed by a method described below. Formerly another source of 
error was found in the fact that whenever the temperature of the 
condenser was changing rapidly, the temperature of the inner cylinder 
would lag behind that of the outer, particularly when the condenser was 
evacuated, because of the extremely good heat insulation provided by 
the quartz disks; for this reason it was found necessary to take observa- 
tions after steady temperatures had been reached or under conditions 
of slowly varying,temperature. Further by the first method it was neces- 
sary to make an experimental run with the condenser evacuated in order 
to determine the change in K due to the thermal expansion of the gas 
condenser. By the second method of procedure the importance of this 
temperature effect is eliminated since now the effect appears only as a 
change in C and Co, whereas before it appeared also as a correction in 
AK, a quantity of the same order of magnitude as the correction. In C 
there is produced thereby a change of the order of 1/10 per cent, which is 
neglected in the computations. 

The preparation and purification of the three gases was achieved in 
the following manner. HCl was generated by dropping concentrated 
H.SO, slowly onto concentrated HCl. The gas thus generated was 
passed consecutively through concentrated H2SO,, a P.O; tube, and a 
liquid air trap where it was frozen. Then the liquid air trap was evacuated 
over the frozen HCI to remove any non-condensible gases such as air and 
hydrogen. Before letting the gas into the condenser tube a portion of it 
was boiled off and removed in order to eliminate possible impurities of 
lower boiling than that of HCl. Then the HCI was allowed to boil and 
to pass into the condenser until atmospheric pressure was reached, the 
stopcock to the generator was opened, and the surplus gas was allowed to 
bubble away through a trap of concentrated H2SO,. Other possible 
impurities would boil off last and have little chance of diffusing through 
the capillary tube running to the condenser. In this way the same gas 
could be used again and again. It may be noted that on renewing the 
gas no change in € occurred. 

For the generation of HBr bromine was dropped slowly into a moist 
mixture of sand and red phosphorous, and the gas evolved was passed 
through a tube of red phosphorous to remove free bromine vapor and to 
the P.O; tube and the liquid air trap where it was distilled as was the 
HCl. The surplus gas was bubbled through concentrated HBr. 





ELECTRIC MOMENT OF HALOGEN HYDRIDE MOLECULES 411 


HI was generated similarly by dropping a thin mixture of water and 
red phosphorous onto moist iodine. It was bubbled through concentrated 
HI to remove free iodine vapor and then dried and frozen as was done for 
the other two gases. 

Table II shows the experimental and calculated data for these three 
gases. Each value (e—1) represents the mean of from 5 to 10 or more 
readings taken at the given temperature. In the third column is given 
the value 7/273, which is the specific volume relative to 0°C and 760mm 
Hg and calculated according to the perfect gas law. The fourth column 
gives the correction to be applied to this value of v as obtained from the 
constants of the van der Waals equation; and the fifth column, the cor- 
rected value of v. In the sixth column is given the value of («—1) reduced 
by a slight correction for pressure made on the simple gas law. The next 
column shows (e—1)v, which is the value of (e—1) reduced to constant 
density by use of the Clausius-Mossotti relation («—1)v=const. for a 
given temperature. These values of (e«—1) for constant pressure and for 
constant density are shown in Fig. 4 to give an idea of the order of mag- 
nitude of the variations occurring throughout the range of temperature 
studied. ; : 


TABLE II 








Correction 
tov j (e— 1) 760 (e—1)v 





60% i . 007452 . 005462 
07 . 004716 .004494 
08 : . 003792 . 004092 
.22 .002672 . 003524 
23 ‘ .001948 . 003104 
. 23 . 001526 . 002823 
.27 ‘ .001182 .002555 





.004545 .003602 
. 003209 . 003160 
.002796 . 003020 
. 002228 .002773 
. 001303 . 002295 
. 000943 .002081 


mNORON| COONN +! 
++t+i+]++4++4/ | 





. 002687 . 002402 

.002123 .002296 

.001760 .002231 ; 

. 001239 .002153 1.022 
.0009706 = .002177 k. 


wNAawoewN 


612. 








* No values of the van der Waals constants are available for this correction. 
t The uncorrected value of v is used here. 
The values in the last column of Table II serve as a test of the Debye 
equation as follows. This equation can be written as 
(e-1)vT=AT+B 
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which shows that if the values of the last column be plotted against the 
absolute temperature there should be obtained a straight line with 
intercept B determined by yu, the electric moment, and with slope A 
depending upon the contribution of the elastic electrons to the polariza- 
tion. 
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Fig. 4 


These points are plotted in Fig. 5. For HCl and HBr straight lines 
could be drawn so that none of the points deviate by more than 1/2 
per cent of (e—1), even at temperatures near the liquefaction point. 
This agreement cannot be showm adequately by the necessarily small 
scale drawing of Fig. 5. For HI two sources of difficulty arise. At high 
temperatures this gas dissociates to a marked degree. In Gmelin and 
Kraut" the equilibrium degree of dissociation of HI is given as about 
3 per cent at 0°C and 18 per cent at 300°C. In the experiments described 
this degree of dissociation was probably not reached at the moderately 
high temperatures on account of the slowness of the reaction velocity 
at such temperatures. (Any dissociated gas was removed after each 
reading by the purifying process so that.no accumulation of dissociated 
gases could take place.) Since the number of molecules does not change 
in the dissociation, except possibly by a slight amount due to condensa- 


‘Gmelin und Kraut, Handbuch der anorganischen Chemie 1, (2), p. 321. 
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tion of the iodine on the walls of the tube, this consideration becomes 
less important. Also from data on the refraction indices of HI, He, and 
I it is seen that the refractivity of HI is about ten per cent greater than 
would be expected on the principle of additivity of refractivities. Hence 
in dissociation the effect of the destruction of the moment would in part 
be compensated for by the increase in the contribution from the elastic 
electrons. However, the experimental values of (e—1) give a much better 
straight line than might have been hoped for. With the above considera- 
tions in view a straight line was drawn through the lower three points, 
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Fig. 5 


where very probably the above mentioned errors are only very small; 
the fourth point is fairly good but the last point is considerably off the 
straight line. Another difficulty in the case of HI is that no data are 
available on the constants of the van der Waals equation. The un- 
corrected values of v were used here. This introduces a possible error of 
about 1/2 per cent. 

Table III gives the results obtained from the lines of Fig. 5. The 
second column under (¢— 1)» gives the contribution of the elastic electrons 
to (e—1); and the third column, values of (n?—1) for the Na D line, since 
(n?—1) for long wave-lengths is not known. These two columns should 
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be identical if the values of the refractive index were given for long wave- 
length instead of for the Na D line. The last column gives the values 
of the electric moment as calculated from the intercept of the lines of 
Fig. 5. 

B= (4r/3)(Nou?/3R), (using e+2=3) 

a *=(3X1.36 10/42 X2.71 X10") B=1.198 X10 B 


TABLE III 
Results for HCl, HBr and HI, 








Gas (e—1) (n? —1) u X10" 





HCl .001040 .000888 1.034 
HBr .001212 .001140 . 788 
HI 001856 .001812 . 382 








DISCUSSION OF DATA 


Of the three dipole gases here studied only HCI has been worked upon 
before so far as is known by the author. Badeker® attempted to study HCI 
over a considerable range of temperature but was prevented by chemical 
action which he attributed to moisture. However, his values of («—1) 
for 84.4° and 106.6°C, .000269 and .000256 show very good agreement 
with the corresponding values from Fig. 4, .000270 and .000247, better 
than might be expected considering the fact that Badeker did not put 
much confidence in his values. These values for HCI do not extend over 
a sufficiently large temperature range to allow the calculation of the 
electric moment with any degree of certainty. 

It should be mentioned here that recently Falkenhagen” has published 
a value of the electric moment of HCI, »=2.1510-"* c.g.s.u., obtained 
from dielectric constant measurements. This is about twice the value 
givenin Table III. Since no values of ¢ are given and only a brief descrip- 
tion of the method, it is impossible to arrive at an explanation of the dis- 
crepancy. 'A calculation has been made of the contribution of the doublets 
alone to («—1) using the value given by Falkenhagen for the moment. 
This gives B=y?X 10*°/1.198 =3.86; hence B/T =3.86/373 =.0104 for 
100°C. Thus it is seen that Falkenhagen’s value of the moment leads to a 
value of electric polarization about four times the total value observed 
by the author and also by Biadeker. Hence the value given by Falken- 
hagen would seem to be impossibly large. 

Frivold and Hassel'* have made a determination of the dipole length 
of the HCI molecule by means of electrostriction. They give for the value 
of the dipole length 0.31108 cm. This corresponds to a value of 


'? H. Falkenhagen, Phys. Zeit. 23, 87 (1922). 
QO. E. Frivold und Odd Hassel, Phys. Zeit. 24, 82 (1923) 
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u=1.48X10"* c.g.s.u., which is considerably less than the value of 
Falkenhagen. These electrostriction measurements are of very limited 
accuracy. 

In support of the value obtained by this investigation Smyth" has 
calculated the electric moment of HCI by a simplified use ” of an equation 
developed by Gans.” This equation applies to both the gaseous and the 
liquid states and is developed on the dipole assumption taking into ac- 
count the intermolecular forces. By Smyth’s method the electric moment 
can be calculated approximately from a knowledge of the indices of re- 
fraction at two different wave-lengths, the corresponding density, the 
dielectric constant at one temperature, the density at this temperature, 
the molecular diameter, and certain universal constants. For the value 
of the dielectric constant he used Badeker’s two values for the gaseous 
state and Schlundt’s value for the liquid state at 27.7° C and obtained 
for uX10'* at these three values, 1.096, 1.054, and 1.075 c.g.s.u. respec- 
tively. The mean of these three vatues, 1.075, agrees with the value of 
Table III, 1.034, better than could be hoped for taking into consideration 
possible errors in the calculations and in the experimental values in- 
volved. 

A further argument in favor of the value 1.034 is the fact that («—1)>, 
in Table IIT agrees fairly closely with the refractive index data. 


RELATION TO INFRA-RED SPECTRA 


Values of the electric moment of the halogen hydrides are of particular 
interest for two reasons; first their structure is perhaps the simplest 
dipole structure because of the simple hydrogen atoms, and second their 
nuclear separation (except for HI) has been determined from infra-red 
absorption data.’’ Now it will be noted that once the nuclear separation 
is known a knowledge of the average position of the electronic configura- 
tiun obtained from the value of the electric moment woufd be of consider- 
able importance in the study of the structure of the halogen hydrides. 
Suppose that in the formation of the HCI molecule, got example, the 
electron from the hydrogen atom combines with the seven outer electrons 


of the chlorine atom to form a group similar to the outer group of the 
argon atom. If this argon group be assumed to be symmetrical and not 
appreciably distorted by the presence of the hydrogen nucleus, then the 


4 This work is not yet published. 
% C, P. Smyth, Phil, Mag. 45, 850 (1923). 
* R. Gans, Ann. der Phys. (4) 44, 481 (1921); 
H. Isnardi and R. Gans, Phys. Zeit. 22, 230 (1921). 
17 See A. Kratzer, Zeit. f. Phys., 3, 305 (1920); 
Herbert Bell, Phil. Mag. p. 549, March 1924 
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chlorine atom plus its acquired electron would be equivalent to a single 
negative charge at the chlorine nucleus. This latter charge combined 
with the hydrogen nucleus gives an effective moment of HCI such that 
its dipole-length is equal to its nuclear separation. But one would expect 
the hydrogen nucleus to distort the electronic orbits in such a way as to 
decrease the dipole-length. Hence it seems reasonable to regard the 
nuclear separation as an upper limit to the dipole-length as pointed out 
by Pauli. C. P. Smyth’® has made a calculation for the effect on the 
electric moment of the HCI molecule due to the distortion by the hydrogen 
nucleus. For this purpose a statical structure of the Langmuir type was 
assumed and the distortion of the individual electrons was calculated 
from the coefficient of elastic binding given by the refractive index of 
the substance. These calculations give the proper order of magnitude 
for this distortion. From thé point of view of the Bohr theory one would 
not expect a calculation based upon the above assumptions to give a 
better agreement. 

Such a consideration is of interest in connection with the quantum 
theory of dielectric polarization developed recently by Pauli.1® Whereas 
Debye and Thomson used in the development of their formulas classical 
mechanics, Pauli has derived for the case of dipole gases a formula based 
on Sommerfeld’s quantum conditions. The latter theory leads to an 
equation of the same form as that of Debye, differing only in the 
numerical factor which appears in the doublet term after integration. 


e—1 





= (40/3) N{de?+1.5367 w?/RT} . 


e+2 





Thus so far as moments are concerned the two theories would give dif- 
ferent values for the electric moment calculated from data on the dielec- 
tric constant, such that: 
u (classical) = 2.15 w (quantum). 

The measurement of the dielectric constant of dipole gases at different 
temperatures offered a possibility of distinguishing between the two 
theories, should the classical theory lead to a value of the dipole-length 
greater than the nuclear separation. 

In Table IV are given the results of some calculations relating to the 
above mentioned quantities. d isthe dipole-length calculated on the clas- 
sical theory, and ro is the nuclear separation obtained from the infra-red 
data. D is the distance of the center of gravity of the electronic configura- 


‘8 C, P. Smyth, Phil. Mag., 47, March 1924 
'° W. Pauli, Jr., Zeits. f. Phys. 4, 319 (1921) 
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tion from the halogen nugleus; this may be regarded as a measure of the 
distortion caused by the hydrogen nucleus. A simple calculation gives 
D=(ro—d)/(N+1), where N is the atomic number of the halogen atom. 
Since the value of the dipole-length is about one-sixth the value of the 
nuclear separation it is impossible to choose between the quantum theory 
and the classical theory. 


TABLE IV 


Comparison of calculated dipole length d with nuclear 
separation ro from infra-red data. 





Gas X10" d=* ro X108 





HCl 1.034 217. 1.2 
13 


HBr . 788 . 165 


7 
7 
HI . 382 . 080 : 








In conclusion the writer wishes to acknowledge a grant from the Joseph 
Henry Fund of the National Academy of Sciences, which has in part made 
possible the construction of special apparatus for this problem. The 
writer is also deeply indebted to Professor K. T. Compton whose stimulat- 
ing interest and whose most generous help and advice have contributed 
largely to the success of this investigation. 


PALMER PHysICAL LABORATORY, 
PRINCETON UNIVERSITY, 
May 8 1924 
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THE MAGNETIC SUSCEPTIBILITY OF HELIUM, 
NEON, ARGON, AND NITROGEN 
By L. G. Hector 
ABSTRACT 

The apparatus described in a recent paper by Wills and Hector! for the 
measurement of the magnetic susceptibilities of gases, has been rebuilt so as 
to enable the temperature to be controlled to .02° and the sensitivity of the 
method increased. The gases used were all carefully purified, and spectroscopic 
examination showed no oxygen lines. The results obtained for the volume 
susceptibilities at 20°C and 760 mm pressure for helium, neon, argon and 
nitrogen are —0.780, —2.77, —7.52, and —4.91, respectively, all times 10—°. 
They are probably accurate to within 2 per cent. These results are of the 
order of magnitude required by Pauli’s modification of the Langevin theory 
of diamagnetism and also agree with values predicted by Joos from the electron 
arrangements. 


INTRODUCTION 


N a previous paper by A. P. Wills' and the writer a method was devel- 

oped for determining the susceptibility of fairly strong magnetic gases 
by balancing them magnetically against a water solution of nickel 
chloride. It was further shown that the method could be extended to the 
case of the nearly neutral diamagnetic gases by balancing the gases 
magnetically at various densities (pressures) against the same solution 
at slightly different temperatures, since the magnetic temperature 
coefficient of the solutions was large and that of these gases relatively 
negligible. The method was carried out for the cases of hydrogen and 
helium, but only tentative results were obtained for the latter gas be- 
cause of its extremely small value. The present research involves per- 
fecting the details of the method in order to be able to determine such 
small susceptibilities as that of helium, and applying the method to the 
measurement of the other gases, argon, neon, and nitrogen. 

It was shown in the previous paper that the working formula for this 
method of determining the susceptibilities of diamagnetic gases is: 


PA PwXw dRy Aé 


ee _ (1) 
6, 1+R-—K do AP 





gq 


where xo is the volume susceptibility at the absolute temperature 4, 
and standard pressure Py) in mm mercury, and Rp is the ratio of mass of 


1A. P. Wills and L. G. Hector, Phys. Rev. 23, 209 (Feb. 1924) 
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water to mass of nickel chloride for a magnetically neutral solution at 
the absolute temperature @ at which the experiment is performed. p, is 
the density and yx, the volume susceptibility of water at the tempera- 
ture 6. 

R is the actual ratio of the mass of water to mass of nickel chloride in 
the solution used, and may be replaced by Ro without appreciable error 
for these nearly neutral solutions. K is a constant depending solely on 
the chemical concentration and density of the mother solution from which 
all the solutions used were made by dilution by weighing. 

The mother solution used was the same as that employed in the 
previous research, and hence we have the same value of K, 0.88102. We 
have also the same Curie constant (C= R,@) for this solution of nickel 
chloride. As before specified, there was about 0.08 per cent cobalt chloride 
in this solution, and the Curie constant obtained was C=13890. The 
value of x, is taken to be —0.72010~°. If we insert these values in 
Eq. (1) we have for standard pressure (760 mm Hg) and a temperature of 
20°C 
Se Aé 


K o> aaa 0.02594 oa 
0.1190 6+ 13890 AP 


where A@ AP is to be determined by experiment and expressed in degrees 
per mm Hg pressure. The temperature @ must be observed at the same 
time and the corresponding value of p, taken from a table of water 
density. 

The method employed was to determine a fairly large number of 
values of 6 and the corresponding values of P for a given solution. These 
values were then plotted, and the quantity 4@/AP determined from the 
graph by an approximate least squares solution. 


PROCEDURE 


Oxygen coming from the air dissolved in the solution formed a principal 
source of trouble at the beginning of our attempts to measure the diamag- 
netic gases. However, by means of a vacuum pump a pressure of about 
1 cm of mercury was held on the balance for about one half hour at the 
beginning of a run. This process hastened the removal of dissolved air 
from the solution. The apparatus was then flushed with hydrogen for 
from two to three hours, and then with the gas to be measured for about 
one half hour. 

The apparatus was similar to that used in the previous research with 
the addition of a double-walled house about the entire magnet. The 
temperature of the air within the house was controlled to within plus or 
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minus 0.02°C by means of a bi-metallic strip which applied potential to 
the grid of a three electrode tube, in the plate circuit of which was a high 
resistance relay. Temperature adjustment was effected by means of a 
bank of lamps controlled by the relay, and by means of pails filled with 
ice and placed in openings in the top of the house. An electric fan kept all 
the air in the house in circulation. 

The gases (except neon) were supplied to us under high pressure in 
steel cylinders. Pressures, except for neon, were read by means of an 
ordinary gauge calibrated against the mercury manometer described in 
the previous paper. Absolute values of temperature were read from a 
mercury thermometer calibrated to 0.02°C, placed between the poles of 
the magnet, very near to the balance. The small variations in the tem- 
perature were read from the constant deflections produced on a galva- 
ometer by means of a copper-constantin thermopile. One set of junctions 
was kept in'a Dewar bottle filled with ice and water, and the other set 
was wrapped about the horizontal tube of the balance and protected from 
the air by cotton wool and tape, so that variations in temperature would 
come from the liquid within. The temperature variations in the liquid 
were produced by slightly varying the temperature of all the air inside 
the double-walled house. These variations were small, seldom exceeding 
0.2 or 0.3° at a time. 

All the gases measured were examined spectroscopically for oxygen, 
but in no case were any oxygen lines observed. 


HELIUM 


The helium supplied to us contained about 8 per cent nitrogen. It was 
purified by passing it through a long tube filled with magnesium granules 
and kept at a dull red heat by means of an electric furnace. This proce- 
dure served to remove any traces of oxygen as well as the nitrogen. 

The following tables give the data from which the graphs and calcu- 
lations were made. Galvanometer readings and gauge pressure readings 
are omitted and only the corresponding columns of corrected pressures 
and temperatures are given. The numbers in the temperature column 
are in centigrade degrees, but are referred to an arbitrary zero which is 
different in the different tables. Pressures are in mm Hg at 0°C. 

Fig. 1 shows the (@—P) graph for only one set of data for each gas. All 
the values of A@/AP given have, however, been taken from similar graphs 
of the various sets of data given in this paper. 
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NITROG 
(TABLE 


ARGON 
(TABLE 


HELIUM 
(TABLE 








. 2345 678 4 WH 12 
PRESSURE (IN METERS HG.) 
Fig. 1. 
TABLE | TABLE II 
Results for helium Results for helium 
Pressure Temperature Pressure Temperature 

1230mm 94° 1180mm 3.17° 
6000 a2 4950 3.31 
9000 1.26 8450 3.46 
11440 1.36 10880 3.70 
9280 1.27 10875 3.64 
7870 1.25 10570 3.58 
6200 1.15 7770 3.51 
2930 .99 5950 3.46 
1180 .96 5200 3.43 

3820 3.34 

3220 3.32 

2400 3.31 

1230 3.26 

930 3.25 

From the data of Tables I and II with the aid of Eq. (2) we find the 





following values 


6 Pw A6/AP Ky 
(1) 292.2°K 9985 4.12«K10° — 0.766 &107'° 
(Il) 294.1° 9980 4.28 —(.795 


Mean —0.780 x 10~-'° 
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ARGON 


The argon as supplied to us contained about 4.5 percent nitrogen and 
0.5 per cent oxygen. The purifying train consisted of heated tubes, one 
filled with magnesium granules, the other with copper filings which had 
previously been oxidized and reduced. 


TABLE III TABLE IV 
Results for argon Results for argon 
Pressure Temperature Pressure Temperature 
1680mm 1. 86° 1375mm 63” 
2250 2.4% 4740 2.33 
3320 2.47 4920 2.38 
3875 2.69 4590 2.42 
3950 2.82 4180 2.03 
3580 2.67 3500 1.99 
3080 2.41 3120 1.79 
2750 2.34 2850 1.63 
1980 2.07 2600 1.46 
1380 1.80 2140 1.26 
990) 1.62 1680 1.07 
1375 1 
1200 . 86 
1080 74 
900 48 


4 


From the data of Tables III and IV, with the aid of Eq. (2), we obtain 
the following values 


6 Du Aé/AP Ko 
(111) 291.7°K 9985 4.10107! —7.63 X10-" 
(IV) 290.1 ORY 3.98 —7.41 


Mean —7.52x10-" 
NITROGEN 


The nitrogen as supplied to us was supposed to be better than 99.9 
percent pure, but as a final precaution against oxygen we passed it 
through a heated tube filled with oxidized and reduced copper filings. 


TABLE V TABLE VI TABLE VII 

Results for nitrogen Results for nitrogen Results for nitrogen 
Pressure Temperature Pressure Temperature Pressure Temperature 
4720mm 2.68° 3710mm 1. 58° 7250 3.44° 
7140 3.43 6810 2.15 9950 4.33 
10570 4, 38 8220 2.63 11340 4.75 
11340 4.59 9890 3.16 8030 3.93 
11700 4.63 11340 3.48 5350 3.20 
9330 4.02 11640 3.56 2700 2.49 
7970 3.64 10410 3.28 1320 2.08 
5400 2.98 8890 2:93 

2700 2.23 7250 2.48 

1030 1.80 5750 2.02 

3380 2.40 3250 1.45 

1860 a3 


1230 .88 
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From. the data of Tables V, VI, and VII, with the aid of Eq. (2) we 
tind the following values 


6 Pu Ad/AP Ko 
(V)  291.2°K ~—_.9986 2.68 X10-' ~ 4.98 x 10-” 
(VI) 290.6° .9987 2.57 —4.78 
(VII) 291.4° .9986 2.67 — 4,97 


Mean —4.91 x 107, 
NEON 


The neon was purified by passing it through a trap packed with char- 
coal and immersed in liquid air. Quantitative chemical analysis showed 
no oxygen present, but there was a small amount of residual nitrogen and 
some helium. These impurities were estimated at 2 per cent or less 
nitrogen and 3 per cent helium. Such a mixture would have a suscepti- 
bility so near that obtained for neon that any errors arising from con- 
sidering the mixture to be 100 per cent neon would be within the experi- 
mental limits of precision. As a final precaution against oxygen a tube of 


hot reduced copper filings was kept in the gas train throughout the 
measurements. 


as 

|! MANOMETER TUBE. 

‘| (4% M, HIGH) 
a 
TO VACUUM 
TO PURIFYING TUBE PUMP 
AND BALANCE 
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Fig. 2. Apparatus used in compressing neon. 


The neon was supplied to us in 350 cc glass flasks at atmospheric pres- 
sure. Fig. 2 shows the method of compressing this gas in the balance 
without any opportunity of introducing impurities. 








(a) 


Pressure Temp. 
2017mm 2.53° 
1349 2.44 
910 2.30 
1379 2.47 
1967 2.52 
3572 _ 2.82 
3213 2.76 
2655 2.72 
2117 2.63 

(a) 
Pressure 
808mm 

908 
1906 
2964 
3313 
2405 
1576 
848 

(a) 

Pressure Temp. 
836mm 1. 86° 
1892 1.99 
2919 2.13 
3288 2.27 
2092 2.06 
796 1.89 


Temp. 


—me bho hte ee 
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TABLE VIII 


Results for neon 


(b) (c) 
Pressure Temp. Pressure Temp. 
1867mm 2.65° 3532mm 2.95° 
1290 . waa 2697 2.84 
896 2.57 1808 2.75 
2236 a.42 871 2.60 
2824 2.82 
3074 2.90 
3562 2.95 
TABLE 1X 


Results for neon 
(b) 


From these data with the aid of Eq. (2) we find 


7] 
VIII (a) 291.0°K 
(b) 291.3° 
(c) 291.4° 
IX (a) 292.3° 
(b) 292.4° 
. | (a) 292. 8° 
(b) 293 .0° 
(c) 293.0° 


Pressure remp 
a 848mm 1.91° 
.80 1526 2.02 
91 2305 2.14 
11 3053 2.26 
.19 3323 2.29 
.06 2305 2.16 
.95 987 1.95 
. 84 
TABLE X 
Results for neon 
(b) (c) 
Pressure Temp. Pressure Temp. 
796mm 1.99° 816mm 2.03° 
1922 2.32 1484 2.43 
3248 2.34 2760 2.29 
3059 2.30 3284 2.41 
2172 2.19 2042 2.20 
816 2.02 1055 2.09 
p,, A6/AP Ko 
. 9987 1.72xX10~ —3.20 107" 
. 9986 1.56 —2.90 
. 9986 1.26 —2.4¢ 
.9985 1.31 —2.81 
. 9984 1.53 —2.85 
. 9983 1.48 —2.75 
.9983 1.35 —2.51 
.9983 1.47 —2.73 
Weighted mean —277. X10 


In obtaining the mean, the mean values from each of the three tables were given 


equal weight. 


TABLE XI 
Volume susceptibilities at 20° C and 760mm pressure 
Ko «10 
Helium —0.780 
Neon —2.77 
Argon —7.52 
Nitrogen —4.91 
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The neon, argon, and nitrogen were made by fractional reduction of 
liquid air by the Air Reduction Company of New York City. I am in- 
debted to Dr. F. J. Metzger of this company for providing me with 
exceptionally pure samples for our measurements. In the case of argon 
and nitrogen, which were furnished us under high pressure, water lubri- 
cated compressors were used to avoid contamination with oil vapors. The 
helium was provided by the Bureau of Aeronautics, Navy Department. 


SENSITIVITY OF APPARATUS AND PRECISION OF MEASUREMENTS 


A variation of plus or minus about 0.03° from the ideal zero setting 
corresponded to the smallest magnetic difference for which a lack of mag- 
netic balance could be observed. The corresponding change in the sus- 
ceptibility of the solution is about 0.7 x 10-'°. With the field used (25,000 
gauss) and a tube of diameter 4 mm, this gives a force of about 1/300 
dyne as the limit of the sensitivity of the device. These figures give 
some idea of the dependability of the individual points determined. 
But the slope of a line on the (6,P) graph is, of course, determined from 
a number of points. Judging from the agreement of the independent 
determinations, the mean values are probably accurate within +2 
per cent. The accuracy of the results depends also onthe value of the 
susceptibility of water, which I have assumed to be —0.72 X10-* plus 
or minus one half of one per cent. 

DISCUSSION OF RESULTS 

Pauli? has shown, by a modification of the Langevin theory of dia- 
magnetism, that the order of magnitude of the susceptibility of helium 
and argon as based on the knowledge of atomic magnitudes from other 
sources should be less than —2X10-!° for helium at 0°C and between 
—2.7 and —8.1X10-" for argon at 0°C. By a comparison of data on the 
susceptibilities of ions having similar electron arrangements, Joos* has 
also predicted values for these gases, and obtained — 1.3 X 10~"° for helium 
at 0°C and —7.8X10~-"° for argon at 0°C. An extension of his method 
and data would give —3.9X10-'® for neon at 0°C. The agreement be- 
tween the experimental values and these predicted orders of magnitude 
is striking. 

The writer is indebted to Professor A. P. Wills for advice and sugges- 
tions throughout the research. 


DEPARTMENT OF PHYSICS, 
CoLuMBIA UNIVERSITY, 
June 10, 1924. 


* W. Pauli, Jr. Zeit. f. Physik 2, 201 (1920) 
3 George Joos, Zeit. f. Physik 19, 347 (1923) 
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CRYSTAL STRUCTURE OF HARD STEEL 


By K, HEINDLHOFER 


ABSTRACT 


Using an x-ray spectrometer with ionization chamber and reflecting 
Mo Kaand K§ rays from small blocks of steel, spectra were determined for 0.80 
and 1.31 per cent carbon steels both in coarse and in fine martensitic condition. 
The same lines were found as in the spectrum of carbon free iron, but they 
were only one-third as intense and were somewhat broader and slightly shifted. 
Lattice of martensite. These results confirm the previously reported conclusion 
that martensite has a body centered lattice like that of a iron but show that 
the lattice is contracted about 0.4 per cent on a side and is also distorted. This 
distortion, which changes the shape of the lines, disappears largely however 
after heating the martensite at 260°C for 30 min., although the hardness is 
not much decreased. Nature of martensite. From a comparison of the densities 
of the specimens with those of a iron and of cementite, it is concluded that 
martensite is a mixture (1) of a solid solution of carbon in iron, the C atoms 
replacing the Fe atoms, and (2) of finely dispersed cementite. The solid 
solution is unstable even at room temperature. Both constituents seem to 
contribute to the hardness, the solid solution by distorting the lattice and the 
cementite by being so fine grained. 


THEORIES OF THE HARDNESS OF METALS 


HE knowledge of the mechanism of hardening of steel after being 

quenched from a certain definite temperature range is both of 
theoretical and practical importance. We shall recapitulate two im- 
portant theories which seem best to explain hardness of alloys in general, 
very likely including steel as a special case. 

(a) The hardness of solid solutions explained by Rosenhain.' It is a 
well known fact that metals which are mutually soluble in any propor- 
tions do not form alloys of great hardness as compared with that of the 
component metals. Rosenhain points out that the crystalline habit of 
such component metals is very similar. For instance gold and silver are 
soluble in any proportion, their alloys are soft, and both crystallize in 
face centered cube of almost identical size. In this and similar cases the 
atoms of the one metal indiscriminately replace the atoms of the other, 
without materially disturbing the space lattice. 

There are, however, a great number of metals, which show a limited 
solubility. This is because when some of.the atoms of the metal A are 
replaced by the atoms of the metal B, the lattice is disturbed due to 


1 Walter Rosenhain, Proc. Roy. So . A 99, 196 (June 1921) 
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disturbed symmetry. Distortion affects the internal energy of the system, 
which effect increases as there are more disturbing B atoms dissolved. 
A point is reached when the energy contents of the system as a whole 
will be less if the atoms of B are separated and arranged on a space lattice 
of their own. The solid solubility of metal B in A will be as a first approxi- 
mation inversely proportional to the amount of distortion of the normal 
space lattice of A caused by the substitution of an atom of B for an atom 
of A in that lattice. It is the power of the crystal of resisting slip which 
determines the hardness or power of resisting permanent deformation. 
“If this view is correct, then anything which tends to hinder the free 
occurrence of slip on the principal planes of a crystal will increase the 
hardness and will lower the ductility of the material.’’? Since a distortion 
of the lattice as mentioned above must serve as a hindrance to slip on the 
crystal planes, we should expect to find that the hardening effect of one 
metal upon another in the form of a solid solution, should be inversely pro- 
portional to its solid solubility. An example for low solubility is carbon in 
gamma iron (limit about 2 percent by weight). Accordingly the effect 
of carbon on the hardness of iron is great as compared with the effect of 
the highly soluble elements such as nickel or silicon. 

(b) The slip interference theory of Jeffries and Archer*. In addition to 
a theory of hardness of solid solutions which is similar to that given by 
Rosenhain, the authors offer a theory to explain the hardness in general. 
They suggest that any structural condition which interferes with slip 
along the cleavage planes of the crystals of the metal increases its hard- 
ness. 

(1) In a pure metal the most simple source of increased hardness is 
grain refinement, which introduces slip interference at the grain boun- 
daries due to different orientations of the adjacent grains, and especially 
in fine grained metal between the grains. 

(2) Cold working introduces slip interference by the fragmentation of 
the grains. 

(3) Slip within grains is opposed by presence of hard grains uniformly 
distributed as fine particles. 

(4) The effect of a given amount of hard constituent increases with the 
fineness of subdivision reaching a maximum at an average particle size 
denoted by the term critical dispersion (order of magnitude 10-7 cm). A 
higher degree of dispersion, particularly the atomic dispersion of solid 
solutions, is less conducive to hardness. Regarding the special case of 

* Ewing and Rosenbain, Phil. Trans. A 249, (1899) 


’ Zay Jefferies and R.S. Archer, Chem. Met. Eng. p. 1057, June 15, 1921 and p. 242, 
Feb. 8, 1922 
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steel, the authors conclude that martensite (structure of hard steel) 
consists of alpha iron which may hold carbon in solid solution. They claim 
that hardness is due chiefly to the grain refinement of the ferrite (alpha 
iron), but partly to the carbon in solution. 


OBJECT OF THE PRESENT INVESTIGATION 


It seems impossible to furnish by the microscope alone sufficient evi- 
dence to ultimately decide upon the various theories on the structure of 
martensite. The microscope fails to show the distribution of the carbon, 
the presence of which is essential to give iron the quality which allows it 
to be hardened. To bring the problem nearer to its final solution x-ray 
crystal analysis was resorted to by A. Westgren,‘ who in a pioneer investi- 
gation showed that martensite has the same body centered cubic struc- 
ture as the alpha modification of pure iron. He found even the lattice 
parameters of the two to agree within the limit of observation. The only 
difference was that the x-ray spectral lines of the martensite were broader 
and less intense. Since very fine grain colloids give wider spectrum lines 
Westgren concluded that martensite may consist of alpha iron grains 
which are of about the same size as the grains of gold colloids. Undoubt- 
edly such a fine grain size would also tend to explain the hardness of mar- 
tensite. 

Since, however, no martensite can be produced in absence of carbon, 
it seems that carbon must play an essential role in rendering steel hard. 
It was, therefore, of importance to trace the carbon in the iron lattice of 
hard steel. Due to both its low atomic number and the relatively small 
quantity present, it seemed to be hopeless to do this directly. It was, 
however, thought possible that from small changes found in the x-ray 
spectrum indirect conclusions as to the distribution of cai bon might be 
drawn. 

Evidently a method which would show the distribution of intensity of 
the x-ray spectrum, such as the ionization chamber, would be more prom- 
ising for this purpose than the photographic method. It was, therefore, 
decided to study the structure of martensite under different conditions by 
means of an x-ray spectrometer employing the ionization chamber as an 
indicator. 

The brick shaped specimen was placed on a small spectrometer, 
with its polished and slightly etched plane through the axis of rotation of 
the tables. The smallest interval of reading on the vernier is 5’ of arc on 
the chamber slit, and 10’ on the specimen table. The arrangement of 


4A. Westgren, J. Iron and Steel Inst., p. 303, 1921; and London Engineering, May 
19, 1922 











CRYSTAL STRUCTURE OF HARD STEFL 429 


slits and target is shown in Fig. 1. The slits, axis of spectrometer and 
target were lined up with a silk thread. All the slits are cut into sheet 
lead of 3 mm thickness and are non-adjustable. The ionization slit, 
1 mm wide and 5 mm high, is mounted on an arm which is rigidly con- 
nected with the table. 

The ionization chamber is mounted on a support concentric with the 
spectrometer and resting on a reinforced wooden board to relieve the 
spectrometer from carrying its weight. This chamberconsists of an 
earthed brass tank 40 cm long by 10 cm in diameter filled with methy! 
bromide vapor. All joints and pet cocks were sealed with shellac. The 
leak was checked by ionization measurements and was found to be 
negligible even after a months time. The methyl bromide was specially 
prepared and was passed over phosphorous pentoxide and soda lime before 
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Fig. 1. Diagram showing path of x-rays through apparatus 
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being admitted in the tank. Inside this tank was a 3 mm diameter brass 
rod electrode insulated by amber, and surrounded by another electrode 
which is a brass cylinder of oval section. The saturation voltage on the 
plate was found to be about 40 volts at maximum vapor density and 
strongest ionization occurring during the experiments, hence 72 volts was 
considered safe to maintain saturation under all conditions. There are 
two openings 1 cm in diameter sealed by thin mica windows, one at each 
end of the chamber, enabling a convenient lining up of the chamber with 
the slit. Care was taken to prevent the beam from striking the electrodes. 
The wire connecting the chamber with the electrometer is shielded by a 
sealed brass tube of 25 mm diameter. 

The electrometer is of the Dolezalek type, with a sensitivity of 3500 mm 
per volt. On the advice of Dr. K. T. Compton a platinized mica needle 
of only 22 mm diameter was used, which reduced the period of swing of 
the needle at a given sensitivity. This was suspended by a fine Wollaston 
wire. 














430 K. HEINDLHOFER 


The angular speed of the needle at different intensities of ionization 
was observed at the different parts of the scale and at each intensity it 
was found to be constant within one or two per cent at any reading beyond 
40 mm. Readings were thus started at 50 mm from zero. 

Considerable trouble was experienced with the earthing key and all 
the different types tried failed after a certain number of breaks. A clean- 
ing and sharpening of the contacts removed the trouble only for a limited 
time. 

The x-ray tube was of the well known water-cooled molybdenum 
target Coolidge type. It was enclosed in a lead box with walls 1/4 in. 
thick. The x-ray generating outfit and the measuring systems were 
installed in adjoining rooms, the systems being separated by a lead door 
1/8 in. thick. Measurements comparing ionization with the x-ray tube 
turned off and then on, but with specimen removed from the spectro- 
meter showed that shielding from stray rays was perfect. 

The x-ray tube was directly connected with the secondary of the 
transformer, the tube rectifying its current. The tube voltage was kept 
constant by a small variable resistance and a Weston A.C. voltmeter, 
both inserted in the primary of the transformer. The absolute peak value 
of the secondary voltage was determined by the vanishing point of the 
x-ray spectrum at the short wave end and was found to be 42,000 volts 
at 25.3 m-amp. tube load. 


THE SPECIMENS 


The size of all hardened specimens was 1.2 X1X0.4 cm. 

Specimen 1 (carbon free iron). Iron melting bar. This sample was 
almost free from carbon, but contained a considerable amount of silicious 
slag as seen under the microscope. It was not analyzed. To refine the 
grain, it was heated to 900°C, held for six minutes, quenched in oil. 
Treated on March Ist, 1923. 

Specimen 2. Analysis: C =.80%, Cr=.14%, Mn=.35%, Si=.19%. 
Other ingredients were not determined. Treatment 2f (fine martensite). 
Quenched three times from 850°C in oil to dissolve the cementite (Fe3C). 
It was then heated to 765°C to give fine martensitic structure; quenched 
in ice water to avoid troostite. The microscope at 1000 diameters magni- 
fication shows absence of cementite. The martensite needles are very 
fine. No austenite or troostite could be seen. Treatment 2c (coarse marten- 
site). Quenched three times from 850°C in oil; final quench in ice water 
from 1030°C to give coarse structure. Microscope at 1000 diameters 
shows absence of cementite, coarse martensite needles and no austenite 
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or troostite. Specimen 2a (annealed cementite) was cut from the original 
annealed bar containing globular cementite. 

Specimen 4. Analysis: C=1.31%, Cr=.12%, Mn=.24%, Si=.17%. 
Other ingredients were not determined. Treatment 4 f (fine martensite). 
Quenched three times from 850°C and twice from 1035°C in oil to dissolve 
cementite; final quench from 850°C in ice water. No austenite or troostite 
could be seen under the microscope, and but little undissolved cementite ; 
the martensite was fairly fine. Treatment 4c (coarse martensite). Quenched 
three times from 850°C and twice from 1035°C in oil; final quench from 
1035° in ice water. No austenite or troostite could be seen under the 
microscope, and only a few fine grains of cementite; the martensite was 
coarse. In order to avoid decarburization the specimens were soaked for 
only about two minutes after their surface assumed final temperature. 
Specimens 2f, 2c, 4f, 4c, were heat treated at end of Aug. 1923. Specimen 
4a (annealed cementite) was cut from the original annealed bar containing 
globular cementite. 


TABLE I 
ITardness and density of the specimens. 

Specimen Brinell hardness Density 
2a(0.8% C) annealed 140 kg/mm 7.823 g/cm’ 
2f fine martensite 660 7.757 
2c coarse martensite 640 7.737 


2 fine martensite 
drawn to 260° C 590 7.781 


for 30 minutes 


4a (1.31% C) annealed 170 7.793 
4f fine martensite a= 7.729 
4c coarse martensite 640 7.729 


DATA OBTAINED WITH THE X-RAY SPECTROMETER 


Readings were taken in at least 30’ intervals over a range of 36° to 50°. 
Fig. 2 shows the ionization current as a function of the slit angle for 
specimens 1 and 4f. The reading corresponding to the primary beam 
was at 45°. 

The lines with the corresponding atomic spacings of the carbon free 
iron are compiled in Table II. The location of lines check with those found 
by other investigators within close limits and the space centered lattice 
of alpha iron is once more verified. 
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Fig. 2. X-ray spectra for specimens #1 and #4F; ionization as a function of angular 
position of the slit. Peak marked 220a should be 110a(2). 
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TABLE II 
X-ray spectrum of carbon free iron 
Reading 8 Line* Spacing a Deviation 
62° 50’ g° $5’ 1108 2.02914 2.870A +0.21% 
65° 10’ es 110a 2.0305 2.872 +0.28 
70° 25’ 12° 42’ 1008 1.4296 2.859 —0.14 
73° 45’ 14° 22’ 100a 1.4319 2.864 0 
76° 20’ 15° 40’ 21186 1. 1646 2.854 —0.35 
80° 30’ 17° 45’ 21la 1.1661 2.858 —0.21 
86° 5’ 20° 32’ 110a(2) 1.0131 2.865 +0.03 
91° 12’ 23° 6’ 310a 0.9061 2.866 +0.07 
93° 50’ 24° 25’ 3218 0.7608 2.847 —0.59 
95° 40’ 25° 20’ llla 0. 8308 2.880 +0.56 
100° 20’ 27° 40’ 321a 0.7656 2.865 +0.03 
© 27 o4cr 4lla 0.6755 2.865 +0.03 
ae’ oS {tt08(3) Mean 2.864} 





* The line whose peak is measured is designated by the plane and the wave-length 
a=Ka=.711A; B=KB=.629A 

+ The density of extremely pure vacuum fused iron was recently (April 1924) deter- 
mined by the Bureau of Standards to be 7.8685 at 20°C. This gives for a the value 
2 8691A, 0.17 per cent greater than the above mean. The metal had been rolled and 
contained: C 0.001%, S .011%, Si .004%, Mn none, P not determined. These results 
are published with the kind permission of the director Dr. George K. Burgess. 


The main object of taking the spectrum of a iron was, however, for 
comparison with that of martensite. It was found that martensite shows 
the same lines in number and location as soft iron. The most outstanding 
difference is the decrease of intensity of the martensite peaks to about one 
third of the intensity of the iron peaks, and their slight broadening which 
is more pronounced at the right end of the spectrum. Also, particularly 
at the right end of the spectrum, a shift of the martensite lines to the 
right was noticed, indicating a contraction of the martensite lattice com- 
pared with that of pure iron. The formula for the relative value of con- 
traction may be obtained by differentiating Bragg’s formula 2a sin @= , 
giving for small increments, 


Aa/a=—A6\/(2a/d)?—1 


To determine the shift more accurately, peak 310a was gone over in 
detail taking many readings at close intervals, and repeating each reading 
two to four times (see Fig. 2). The martensites with different carbon 
analyses seem to show approximately the same amount of shift. This 
indicates that only a limited amount of carbon can be dissolved (which 
is in accordance with Rosenhain’s view quoted above). The average 
value of the shift of six determinations was found to be 2A@=14’. The 
relative linear contraction is therefore Aa/a = —.00469; and the relative 
volume contraction is Av/y=3Aa/a=0.0141. 211a was treated in a 
similar way and gave a shift of only about 5’ and in the opposite direction 
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to that for 310a, which the writer is unable to interpret at the present 
time. 

Outside of the a iron lines one feeble austenite line at 67°50’, @= 11°25’, 
could be seen on specimen 4c (1.31% carbon, coarse martensite). 

The spectrum of the 1.31% carbon steel was also determined in the 
annealed condition. It is interesting to note that the intensity of the 
lines is somewhat lower than in case of pure iron. This is due to the large 
amount (about 20%) cementite present, which displaces the iron, and 
consequently reduces the quantity of scattering substance. A few cemen- 
tite lines while weak were readily identified. Spectra of coarse and fine 
martensite of medium and of high carbon steel showed no marked differ- 
ences. 


STRUCTURE OF MARTENSITE DERIVED FROM THE FOREGOING DATA 


The x-ray spectrum of the martensite determines its space lattice 
to be a body centered cube as in the case of a iron. The shift which was 
determined on one line only, shows that the martensite lattice is spaced 
approximately 0.4 per cent closer than the a iron lattice. There also is a 
decrease of intensity of the martensite lines to about 0.4 that of pure 
iron lines, accompanied by slight broadening. 

This evidence together with the densities should be considered in con- 
nection with the following various possible assumptions: 

(a) All carbon is present as extremely fine grains of cementite (FesC), 
too small to be visible under the microscope. In this case no distortion of 
the a lattice would be expected. 

(b) The lattice is substitutional, the carbon replacing the iron atoms. 
This would cause a drop of intensity of the spectrum lines due to the 
distortion of space lattice by the comparatively small carbon atoms. 

(c) The carbon is within the lattice. This also may cause distortion 
and might also contract the lattice. 

In order to be able to choose between the assumptions enumerated 
the densities to be expected in the case of each of these assumptions will 
be determined. : 

(A) The density of a mixture of cementitie and a tron. If the mass of 
carbon combined as Fe;C is x, and the ratio of the atomic weights is 
p=Fe/C it can readily be shown that, taking the density of cementite’ 
to be 7.59 the density of steel containing all the carbon as cementite is 

x 


p 7.59 
» Levin and Dornhecker, Ferrum 10, 321 (1914) 


1 
(3p+1) + ——(1-2(39+1)), 
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(B) The density of a substitutional solid solution of the body centered 
lattice type. In the case of a body centered lattice, there are two atoms per 
elementary cube. Since some of the atoms of one kind are substituted 
by the other kind we may speak of an average number 1, mp, of atoms 
of each kind per cube; then 


Neotup=2. 


If the proportional weight of carbon present is y, then the weight ratio 
of carbon to iron is 


y/(A— y)=ncC/neFe . 
Hence the average number of iron atoms per elementary cube is 


2(1i— y) 
nt pe =————- 


“a(p—1)+1 | 


Let p denote the density of the solid solution of carbon in iron, mC the 
absolute mass of a carbon atom and a the lattice parameter, then 


npele+ (2—ne.)C 2 pmC 1 
| — eee * ———— 


a’ a’ x(p—1) 


(C) The density of the martensitic steel under the assumption that all 
carbon is inserted and forms no part of the structure, is 
p=7.86 (1+ y) 


provided the lattice is unaltered. Since the x-ray spectrum shows that 
the lattice contracts, the density would be about 1.014 times greater. 

The densities calculated under the various assumptions are as follows 
for 0.8 and 1.31 per cent C: 


80% C 1.31%C 
(A) Martensite is a mixture 
of submicroscopic Fe;C 7.825 7.803 
and Fe 


(B) Solid solution, all car- 
bon pone the 
iron atoms, lattice 7.637 
assumed to be uncon- 
tracted. 


(B’) Solid solution, all car- 
bon substituting the 7.74 7.60 
iron atoms, lattice con- 
tracted 0.4% 


(C) Solid solution, carbon 
between the lattice, 
lattice assumed the 
same as that of pure 
iron. 


. 490 


~r 


.923 7.963 


~I 
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In Fig. 3 the densities calculated under these various assumptions are 
shown graphically as a function of the relative weight content of carbon. 
The actual densities of the different specimens as determined by weighing 
are also entered in the diagram. These data show that a considerable 
part of the carbon present in martensite replaces iron atoms in the space 
lattice. Due to the fact that the carbon atoms are in the proportion 5:9 
smaller in diameter than the iron atoms which they replace, there must be 
considerable lattice distortion. But a distorted lattice must scatter x-rays 
very like a lattice subjected to temperature agitation. I. Backhurst*® 
found that the peak intensity due to a single crystal of aluminum heated 
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Fig. 3. Density as a function of per cent carbon for various assumptions. 

All carbon in solution and assumed within the lattice; (CJ) 

All carbon present as cementite; (A) 
: Carbon replacing iron atoms in lattice, causing contraction; } (ge) 
: Carbon replacing iron atoms in lattice, without contraction. 
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to several hundred degrees Celsius is considerably less than the cold value, 
the relation of the spectra of cold and hot aluminum thus being similar 
to the relation between the spectra of pure iron and martensite. It is, 
therefore, probable that martensite is a distorted lattice of iron. This 
would also explain its hardness. 

The assumption of colloidal crystal grains would also account for the 
drop of intensity of the martensite lines. In this case the lines correspond- 
ing to the planes of lowest atomic population (fewest number of lines in 
the grating) should drop to the greatest extent in intensity and vice versa. 
This is contrary to what was observed, the drop of line (100) being the 





greatest. 


6]. Backhurst, Proc. Roy. Soc. A 102, 340 (1922) 
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Summarizing the observed facts it seems probable that martensite is a 
mixture, consisting of two constituents which may be present in varying 
proportions, one a solid solution of carbon in a iron, the carbon atoms re- 
placing the iron atoms, the other finely dispersed cementite. 

The solution is unstable and some of the carbon gradually precipitates 
from the lattice even at room temperature; hence it is difficult to keep more 
than part of the carbon in solution. 

The gradual unloading of carbon from the lattice at room temperature 
is shown in Fig. 4 where the 310a line of two specimens is compared, both 
made of the same steel, treated alike, one about five months, the second 
only three or four days old. The intensity of the lines of the former is 





*| 
2F OLD 


2F FRESH 
310% 











90 qr = gi Qe 
Fig. 4. Influence of seasoning on the line 310a. 


slightly greater. To prevent possible error due to slight change of sensi- 
tivity of the electrometer, readings were taken in succession on the three 
specimens shown in Fig. 4, before changing the setting of the spectro- 
meter table. 

Fig. 3 shows that there is a limit of the amount of carbon held in this 
unstable solution. Specimens 4f and 4c have considerably higher density 
than the total carbon in solution would require. Hence the balance is 
likely to be present as dispersed submicroscopic cementite. 

Fig. 5 shows a few drawing experiments. The degree of lattice distor- 
tion indicated by the intensity of lines is small after drawing to a tempera- 
ture of 260°C applied for thirty minutes. The density of this specimen 
increased and approached the all cementite line (Fig. 3). The fact that 
the cementite cannot yet be seen under the microscope may be due to the 
submicroscopic size of the Fe;C grains, but their presence was proved by 
the well known darkening of the polished surface by etching in weak 
acids. 
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It is interesting to note that at the drawing temperature of 260°C the 
lattice distortion has to a great extent disappeared; on the other hand 
the density in Fig. 3 shows that only a part of the carbon has been pre- 
cipitated. The hardness only dropped from 660 to 590 kg/mm*, which 
is far greater than the hardness of the annealed state, 140 kg/mm’. It 
may be that in the undrawn martensite hardness is chiefly due to lattice 
distortion as required by Rosenhain’s theory of solid solutions, but almost 
the same hardness may be maintained by the carbon which is precipi- 
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Fig. 5. Influence of drawing on the line 211a. 


tated at drawing to a moderate temperature, say at 260°C and is likely 
to be present as cementite in fine dispersion, when the distortion of the 
lattice is negligible (Jeffries and Archers slip interference theory). 

In concluding, I wish to express my indebtedness to the S K F Indus- 
tries Inc. for giving me the privilege of devoting part of my time to this 
research, to Prof. C. B. Bazzoni, Dr. Lucian of the University of Pennsyl- 
vania, to Dr. H. Styri and Mr. Frank Wright of the S K F Industries, 
to Prof. C. W. Hewlett and to Prof. K. T. Compton. 

RESEARCH SECTION OF THE RANDAL MorGAN LABORATORY. 


UNIVERSITY OF PENNSYLVANIA, 
May 18, 1924. 
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THE BARKHAUSEN EFFECT 
By E. P. T. TYNDALL 


ABSTRACT 


silicon (4.2%) steel, after having been carried through several hysteresis 
cycles and brought to a steep part of the B-H curve, was subjected to a mag- 
netic field which was increased in a continuous manner by 0.13 gauss in 2 sec. 
To obtain a record of the discontinuities in magnetization, the specimen was 
surrounded by a small search coil connected through an amplifier to a moving 
coil oscillograph. The oscillograph records show many almost instantaneous 
deflections with a random distribution both as to time of occurrence and as 
to magnitude, each followed by an approximately exponential return to zero. 
The duration of an impulse depends probably on eddy currents in the specimen. 
A rough time constant of 3X10— sec. is computed for a 4% silicon steel wire, 
1 mm in diameter. The apparatus was calibrated by means of artificial 
impulses of the same type. Assuming each impulse is due to a sudden satura- 
tion of a small portion of the material, the change in magnetic moment of this 

portion is found to vary from .001 to .008 e.m.u.; the average change .003 is | 
sufficient to saturate a volume of 1.7X10— cc. This volume, while com- j 
parable with that of a crystal grain in the particular case described, was found | 
not to depend upon grain size in any marked degree in other experiments. | 
The results favor the suggestion of Barkhausen that magnetic materials 

magnetize discontinuously, but leave open the question as to what determines 

the size and shape of the portions which suddenly change. 


| 
Barkhausen effect for silicon steel—In a typical experiment a strip of | 
| 
i 
i 





INTRODUCTION 


HE discontinuities in the magnetization of iron and other magnetic 

materials have been investigated by Barkhausen,' Van der Pol,? 
Gerlach and Lertes,* Weiss and Ribaud,‘ Zschiesche,’ and S. R. Williams.°® 
Barkhausen and Van der Pol consider the effect as due to the sudden re- 
orientation of groups or chains of molecular magnets, possibly all the 
magnets in a single crystal turning simultaneously. The latter has pub- 
lished pictures of the discontinuities as recorded by a fluxmeter and by a 
string galvanometer. Gerlach and Lertes, and Zschiesche placed the 
specimen in a solenoid which was then rotated at an arbitrary rate 
between the poles of an electromagnet. The amplified and rectified 





1H. Barkhausen, Phys. Zeits. 20, 401 (1919) 
* B. Van der Pol, Proc. Acad. Amst. 23, 637, (1921); 23, 980 (1922) | 
’ W. Gerlach and P. Lertes, Zeits. f. Phys. 4, 383 (1921) 

‘ P, Weiss and G. Ribaud, Jour. de Phys. (6) 3, 74 (1922) 
* K. Zschiesche, Zeits. f. Phys. 11, 201 (1922) 

*S. R. Williams, Phys. Rev. (2) 22, 526 (1923) 
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current from the solenoid as measured by a galvanometer is defined by 
them as the “effect.” Various experiments led them to the conclusion 
that the effect is dependent on magneto-striction and not primarily on 
discontinuities in magnetization. Williams has recently obtained oscillo- 
grams of the effect using a telephone receiver as an oscillograph. All 
observers seem agreed that the effect takes place only where the B-H 
curve for the material is steep. Preliminary work in this laboratory 
showed that the effect as measured by the noise in a telephone receiver 
is roughly proportional to the slope of the B-H curve. Iron dust, loose 
or compressed, gave a negligible effect. The work reported here was 
undertaken to test the hypothesis that the phenomenon is due to the 
magnetization of single crystals and to obtain a more detailed picture of 
it than has yet been presented. Most of the work has been done on one 
specimen, though enough observations have been made on others to 
give confidence in the generality of the results. 

The discontinuities in magnetization are best measured by the induced 
voltages in electrical circuits, magnetometer methods being unsuited to 
the proper analysis of the effects. In an attempt to record the discon- 
tinuities separately and with as little distortion as possible various 
methods have been tried. The first instrument used was the oscillograph 


tube developed in this laboratory by J. B. Johnson.’ This gave visual 
evidence of the discontinuities but it was at once apparent that photo- 
graphs could not be taken quickly enough to follow the changes in a 
single event. Later a moving coil oscillograph of the ordinary type® was 
used and found more satisfactory. 


EXPERIMENTAL METHOD 


The circuit used in obtaining oscillograms of the effect is shown in 
Fig. 1. Changes in magnetization are produced by passing current 
through either of two single-layer solenoids S; and S2, each 80cm long and 
2.5 cm in diameter, wound in parallel on a glass tube. The search coil 
S3 is a multiple-layer solenoid of 10,000 turns, 2.5 cm long and 2.2 cm in 
outside diameter. The specimen is a straight wire or flat strip longer than 
the search coil and supported at the center of the magnetizing coils S; 
and S:. The voltage generated in S; by sudden changes in magnetization 
of the specimen is amplified by the amplifier A, the output of which 
goes to the oscillograph (Vibrator No. 2). The circuit connected to 
S; is used to carry the specimen through several magnetizing cycles 


1 J. B. Johnson, Jour. Opt. Soc. Am. and R.S.I. 6, 701 (1922) 
* General Electric Company Oscillograph, Type EM Form C. 
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before taking an exposure. The specimen can thus be brought to any 


desired point in the cycle. The experimental change in magnetization is 
then effected by closing the filament circuit of the vacuum tube V. As 
the temperature of the filament rises the current through S. grows con- 
tinuously and almost linearly with time, thus giving a smooth variation 
of magnetizing field. This current is recorded on the oscillograph (Vibra- 
tor No. 1). A time scale is given by a 60-cvcle alternating current (Vibra- 
tor No. 3). No transformers are used in the amplifier circuit, the coupling 
being by means of resistance (100,000 ohms) and capacity (0.1 mf). A 
“V”’ tube® is used in each of the first three stages; in the fourth are two 
“L” tubes,'® connected in parallel. The space current (20) 10-* amp.) 
from the last stage passes directly through the oscillograph vibrator. An 











TO VIBRATOR *| 








TO 
VIBRATOR *2 








Fig. 1. Diagram of apparatus for recording Barkhausen effect. 


amplifier of this type is comparatively free from distortion. The effect 
on the records due to amplifier and oscillograph characteristics is dis- 
cussed below. 

Before taking an oscillogram the specimen is taken about ten times 
around a hysteresis loop wide enough to ensure nearly complete satura- 
tion at either end, and so brought to some desired point, usually that 
corresponding to zero current in solenoid S;. The magnetizing field thus 
produced is not quite zero since the vertical component of the earth’s 
field is not neutralized. The light is then screened from vibrators No. 2 
and No. 3 and the zero position of vibrator No. 1 marked on the film 
during one rotation of the drum. After these preliminaries the oscillogram 
is taken at once to minimize the disturbing effects of stray magnetizing 


® Western Electric Company, 102-DW 
10 Western Electric Company, 216-A 
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fields. The exposure is made just after closing the filament circuit of the 
auxilliary vacuum tube 1’. To simplify the comparison of records two 
standard rates of the oscillograph drum were chosen. In the slower case 
the closing of the key in the filament circuit of V also opens the shutter 
of the oscillograph by means of a relay. The picture thus obtained shows 
the process from the beginning nearly to the end, lasting about two 
seconds. In the faster records the vacuum tube filament is lighted and 
after a suitable short interval the shutter is opened. A record at this 
rate covers half a second or less and a more detailed picture is given of a 


small part of the process. 
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Fig, 2. A record taken at the slower rate. 


DATA AND DISCUSSION 


A portion of a typical film, taken at the slower rate, is shown in Fig. 2. 
This record was obtained with a strip of silicon steel (4.2%) 0.043 cm X 
0.178 cm X17.8 cm which had been annealed in vacuum for 11 hours at 
900°C and had then had a further heat treatment, also in vacuum, of 
30 min. at 1200°C. The average linear dimension of the grains as meas- 
ured on etched areas is about 1.5X10-* cm. The record of Fig. 2 shows 
well the discontinuities in the process of magnetization. The amplified 
curve is approximately the derivative of the B-H curve for the specimen, 
It shows that a portion of the B-// curve consists of short straight lines 
parallel to the B axis and // axis alternately. 
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It might be expected that consecutive experiments under identical 
conditions would give superposable oscillograms, the magnitude and 
order of the impulses being completely determined by the initial state of 
magnetization. In no case, however, was it possible to obtain two identi- 
cal records, even though many successive films were compared, in which 
the specimen had been left undisturbed in position and had been brought 
to the same state through identical preliminary magnetizing cycles. 

A study of the occurrence and magnitude of the discontinuities was 
made. Five successive records like that of Fig. 2 were analyzed. The 
change in the magnetizing field was 0.130 gauss and the change in the 
induction corresponding thereto was about 2000. This increment was 
divided into intervals 0.005 gauss long and the number of impulses in 
each interval recorded. A count was then made of the number of inter- 
vals containing 0, 1, 2, 3 etc. impulses. These data are graphically pre- 
sented in Fig. 3. The ordinate of any point is the number of field intervals 
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Fig. 3. Distribution among equal field intervals of discontinuities in magnetization. 


containing the number of impulses shown by the abscissa. The graph 
shows that this distribution is a random one, except for the excess of 
intervals containing no impulse, most of which occur at the beginning of 
the record. It is thought that this results from the unavoidable fluctua- 
tions in magnetizing field between the end of the preliminary treatments 
and the beginning of the exposure. In Fig. 4 is plotted the distribution 
of the impulses as to magnitude. The ordinate of any point is the number 
of impulses having a magnitude (simply the deflection in cm measured 
on the film) of the size given by the corresponding abscissa. This distri- 
bution also appears to be a random one. 
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Fig. 5. Records taken at the faster rate. 
A. Change in H 0.022 gauss. Measured impulses ¢,f,g. 
B. Artificial impulse, 7 =0.110— sec. 
C. Artificial impulse, 7 =1.45 X10-* sec. 
D, E, F. Measured impulses a,),c,d. 


In a record of the second kind described above each impulse consists of 
i = 
an almost instantaneous deflection followed by an approximately, ex- 


onential return to zero. Part of one such record is shown in Fig. 5A. 
£ 


Oscillograms taken with fresh damping fluid show higher and sharper 
peaks, corresponding to the lower value of damping factor thus ob- 
tained. 
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To calibrate the amplifier artificial impulses were similarly recorded by 
the arrangement shown in Fig. 6. The search coil S;, and the amplifier 
circuit are as before, there being, however, no specimen in S;. In the 
circuit with S, are now an inductometer (Ls, r3), resistance rs, and milli- 


Ww oa 
o o 


af 
oOo 


NUMBER OF IMPULSES 
Ss Ss 





0 2 4 6 68 0 2 14 16 18 20 22 
MAGNITUDE 


Fig. 4. Distribution in magnitude of discontinuities in magnetization, 


ammeter M. By adjusting the ratio of rs to R and the value of Ls it is 
possible to establish a current of a chosen magnitude in the circuit in 
which S: is included and to adjust the time constant of this circuit as 
desired. The equation for this current is given by 


i=] (1—e€*/*) (1) 

in which 
t= (Lot+Ls)/(ratrst+r,) . (2) 
In these expressions J is the steady state current with the key K closed 
and 7 is the time constant. The voltage induced in the search coil Ss, is 


e= Mdi/dt=(MI/r)e-*'" (3) 


where M is the mutual inductance of S; and S;. It isan amplified picture 
of this voltage which is recorded on the film as a function of the time. 
Its initial value and the area under the curve are given by 


én=MI/r and fedt=MI. 
0 


Oscillograms showing two such impulses are shown in Fig. 5B, 5C. 
Table I gives results of measurements on such artificial impulses, the 
deflection h, and the area under the curve A, being the quantities meas- 
ured from the oscillograms. 
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The values of the ratio ¢,/h are interesting in that they show how the 
amplifier and oscillograph distort a very sudden voltage impulse, the 
maximum value being cut down while the rate of return is diminished. 
If there were no distortion the ratio e»/h would be constant. In sp.te of 
the distortion, however, the system is a competent integrator, for, as the 
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l‘ig.6. Diagram of apparatus for producing and recording artificial impulses. 


last column shows, fedt/A is approximately constant for values of 
7 down to 10~* seconds. It is moreover considered constant for values of 
t less than 10~ seconds, as it was later checked experimentally to r=10~ 


TABLE I 

















T em Sedt h A €m/h fedt/A 
(X10-*) = (X 107%) (X10-) (X10-*) (X10-3) (107) 
0.1 sec. 2.93Volts 0.293 0.75 cm 0.64 3.9 0.46 
0.1 6.30 0.630 0.78 0.97 8.1 0.65 
0.2 3.23 0.648 0.90 0.65 3.6 1.00 
0.2 6.25 1.25 1.33 1.26 4.7 0.99 
0.2 1.53 0.31 0.55 0.635 2.8 0.49 
0.2 3.97 0.80 1.02 1.30 3.9 0.62 
0.28 8.16 2.28 2.13 3.55 3.83 0.64 
0.3 3.11 0.923 1.17 1.64 2.66 0.56 
0.41 2.85 ‘27 1.58 1.94 1.8 0.60 
0.60 1.80 1.08 ae 1.64 1.5 0.66 
0.62 3.68 2.28 1.86 3.12 2.0 0.73 
0.63 1.79 1.13 1.60 2.42 1.1 0.47 
1.45 0.764 1.11 1.01 1.78 0.76 0.62 


Average 0.652 














seconds, there being no considerable change in fedt/A between r= 10~* 
and r=10~- seconds. These data are not included in the table as they 
were taken after replacement of the damping fluid in the oscillograph, 
which changed the actual value of fedt/A for any value of r. 
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From oscillograms of the effect taken at the faster rate some specula- 
tions concerning the magnitude of the change in magnetic moment pro- 
ducing a sudden voltage impulse may be made. Impulses of various 
sizes, all from the same specimen previously described and with zero 
current in S;, Fig. 1, were studied. Some of these are shown in Fig. 5 
already referred to. The time constant of any of these impulses is by 
inspection not more than 10~* seconds and may be very much less. Since 
it has been shown that the area under the curve of the amplified voltage 
is practically independent of the time constant of the voltage impulse, an 
estimate may be made of the total change in magnetization correspond- 
ing to any impulse without knowing the time required for the change. 
It will merely be assumed that a change in magnetization of a small 
portion of the material takes place in a short space of time. The portion 
of material will be considered as a dipole whose magnetic moment changes 
suddenly from one value to another. Let M be the variable magnetic 
moment of the dipole and consider it situated at the center of the short 
solenoid $3, Fig. 1. It is then very easy to show that the area under the 
curve for the voltage induced in such a solenoid when M changes from 
0 to M is 








2nM un : 
edit= — —~— (volts~—sec.) 
10°V a?+R® 


where is the number of turns on the solenoid S;, a is half its length in 
cm, and R is the mean radius of the multiple layer winding. Some objec- 
tion may be made to the assumptions that have been made but until 
more is known of the mechanism of magnetization they do not seem 
unreasonable. The dipole has been considered as surrounded by non- 
magnetic material. This assumption, however, does not seem so arbi- 
trary when it is remembered that it is evident from the records that the 
rest of the material does not greatly change its magnetization when one 
small portion does. That a sudden change in magnetization might set 
up eddy currents in the surrounding material which would modify the 
observable effect, is quite possible, and will be discussed later. 

For a given Barkhausen impulse the edt may be determined from 
the area by the calibration factor derived in Table I. The experimental 
data on a series of impulses, shown in Table II, have been interpreted 
on the basis of the analysis given above. The letters a to g designate the 
impulses shown on Figs. 5. For simplicity the area under each deflection 
has been approximated to that of a triangle whose base b is the time 
during which the vibrator is deflected from the zero position. 
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TABLE II 
Impulse b h Area Sedt M 
(X10-%) (10-5) 
a 2.16 sec. 1.00 cm 1.08 0.71 1.5 e.m.u 
b 2.43 2.15 2.61 1.70 3.6 
c 1.69 0.95 0.80 0.52 & 
d 2.60 1.62 2.10 1.37 2.9 
€ 3.04 1.28 1.95 1.27 2.7 
f 1.59 0.65 0.52 0.34 0.72 
g 1.85 1.10 1.02 0.67 1.4 
h 2.58 1.20 1.55 1.01 2.1 
i 1.98 0.72 0.71 0.46 0.97 
j 1.52 0.70 0.53 0.35 7.3 
k 1.06 1.05 0.56 0.37 » Pe 
Average 2.9 








DISCUSSION OF RESULTS 


Assuming that a small volume v of the material is magnetized to 
saturation to produce each of the changes in magnetic moment given in 
Table II, the average value of v is given by 


v=M/c. 


Putting o, the intensity of magnetization at saturation, equal to 1700 
gives 


v=1.7X10-*cm’ , 


and 
v/3=1.210-*cm . 

Although this is not far from the average linear dimension of the grains 
in the specimen, it is not to be concluded that the grains become magnet- 
ized as units. Experiments to be quoted below will make this clear. 
However the change in magnetization arises, the order of magnitude of 
the change has been estimated on the basis of somewhat speculative 
assumptions. 

It will be interesting to compare the result with that obtained by 
another experimenter. Van der Pol assumes groups of elementary mag- 
nets in long filaments and computes that the change in flux corresponding 
to one discontinuity is 0.33 lines. On the same assumption the average 
magnetic moment shown in Table II would correspond to a change in 
flux of 0.14107? lines. Van der Pol used a different material and a 
different method and it is not surprising that the results should differ. 
An experiment was performed to determine whether the filaments which 
he describes existed in the specimen here studied. Two identical search 
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coils, one of which was that previously used, were placed end to end over 
the specimen and substituted in the measuring circuit for the single coil 
S;, Fig. 1. Connections were such that a change in magnetization of given 
sign would cause oppositely directed voltage impulses in the two coils. 
The record showed deflections on both sides of the zero line, similar 
in type and approximately equal in number and magnitude to those 
obtained with a single search coil, clearly indicating that few changes in 
magnetic moment affected both coils. This was taken to mean that this 
sample differed from that studied by Van der Pol in having few long 
filaments or at least in having few that became magnetized as units. The 
result of this experiment justifies the assumption that the magnetic 
element may be treated as a dipole. ; 

To test the effect of grain size on the phenomenon a strip of the same 
silicon (4.2%) steel sheet was vacuum annealed at about 1300°C for 
30 minutes. This produced large grains, extending through the thickness 
of the strip. A specimen 4.3 cm X 0.238 cm X 0.043 cm was sawed from the 
portion where the grains were largest. About half this strip consists of 
one large grain. Another grain is about half as large, and altogether there 
cannot be more than twenty-five grains. That such large grains really 
represent single crystals was shown by examining the etching pits, which 
were alike over any one grain, and was checked by making a Laue photo- 
graph." The records taken with this specimen do not differ more from 
those taken with the specimen first described than is to be expected in 
view of its shortness and consequent smaller change of induction for a 
given increment of magnetizing field. The impulses occur in the same 
random manner and are of about the same general intensity and number 
as before. A single grain, therefore, can hardly be conceived of as chang- 
ing its magnetization as a unit. 

If a dipole be instantaneously created in the interior of a conducting 
mass, the sudden growth of the magnetic field will set up eddy currents. 
The field of these currents will oppose the growing field of the dipole. A 
solenoid surrounding such a system would have a voltage pulse generated 
in it, consisting of an instantaneous rise to a finite maximum (initial) 
value and a gradual decay. The time of decay would depend on the time 
constant of the eddy current circuit. A paper by Wwedensky” enables a 
rough approximation of the time of the decay of the Barkhausen im- 
pulses to be made. Wwedensky considers an infinite cylinder of magnetic 
material, whose permeability is constant, subjected to the influence of a 


4 By Dr. R. M. Bozorth of these laboratories. 
1 B. Wwedensky, Ann. der Phys. (4) 64, 609 (1921) 
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longitudinal magnetic field whose intensity changes instantaneously to a 
higher or lower value. He obtains an expression for the reduced mean 
flux through the cylinder as a function of time. Defining a sort of rough 
time constant as the time, r’, for the growing flux to reach (1—1/e)th of 
its final value, and using Wwedensky’s expression, a silicon steel wire of 
0.05 cm radius (permeability 400, conductivity 0.0210-* e.m.u.) has 
r’=3 10> sec. A lower value of permeability or smaller radius would 
give a smaller value of 7’. 

The impulses of Table II were produced in a strip only 0.043 cm thick 
and the corresponding value of 7’ although not directly calculable may 
be 10-° sec. or even less. An attempt was made to detect the effect of 
eddy currents by comparing records taken with two specimens of silicon 
(3%) steel using fresh damping fluid in the oscillograph. One specimen 
was a wire 0.16 cm in diameter; the other was a tape 0.30 cm wide and 
0.048 cm thick. Under these experimental conditions it was observed 
that the impulses were somewhat broadened in the case of the wire. The 
time constant of the deflection as measured on the oscillograph did not in 
either case differ much from 10~ sec. The energy dissipated by such 
transient eddy currents is to be considered inseparable from the process 
of magnetization and not, like that dissipated by the eddy currents 
ordinarily considered, a function of the frequency of alternating magnet- 
izing fields. There should, therefore, be a difference between the energy 
loss in low frequency magnetization as measured by calorimetric methods 
and that measured by plotting stationary values of induction and magnet- 
izing field, in addition to the difference ascribable to eddy currents of the 
same frequency as the magnetizing field. 


CONCLUSION 


The above results favor the theory advanced by Barkhausen and 
upheld by Van der Pol that magnetic materials magnetize discontinu- 
ously, but do not substantiate the belief that each discontinuity is due to 
a single crystal. What determines the selection, number, and grouping of 
atoms which may at any time form a simultaneously magnetizable unit, 
and what determines the order of response of such units to an externally 
applied field, are not apparent, nor can any conclusion be drawn as to 
whether a given unit magnetizes to saturation in one jump or in several. 
That there should be any continuous variation of magnetization between 
the discontinuities seems unlikely. The apparently smooth parts of the 
curve observed by Van der Pol with the fluxmeter can be considered as 
‘regions in which a large number of small discontinuities occur, that is, the 
phenomenon can be regarded as too fine-grained for the method of obser- 
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vation. Little justification has been found for the theory proposed by 
Gerlach and Lertes and supported by Zschiesche. Their adoption of a 
complicated and somewhat obscure explanation depends mainly on thei: 
failure to detect an effect in iron powder and divided iron bars. With iron 
powder of the sort here examined even saturation of individual grains 
would give discontinuities too small and too numerous to give effects 
measurable either by their method or by that here described. It seems 
best to adopt the simpler notion that for some reason, limited portions of 
the material magnetize by jumps, either partially or to saturation, and to 
leave open the question of what determines the size and shape of these 
portions in various materials and in various specimens of the same 
material which have undergone different treatments. The identification 
of these individuals and a better understanding of their behavior would 
give microscopic detail to pictures of the mechanism of magnetization. 
RESEARCH LABORATORIES OF THE 
AMERICAN TELEPHONE AND TELEGRAPH COMPANY 


AND THE WESTERN ELEcTRIC COMPANY, INCORPORATED. 
April 21, 1924. 
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ACCURATE DETERMINATIONS OF THE SPEED OF 
SOUND IN SEA WATER* 


By E. A. ECKHARDT 


ABSTRACT 


Sound and radio signals were transmitted simultaneously to receiving 
apparatus at a distance of over 60 miles. The distances were determined to 
within 6 ft and the times to within .001 sec. so that the mean of 8 independent 
values obtained within 10 minutes is accurate to about 1 in 50,000. With 
practically no tide, the mean result was 4896.1 ft/sec. or 1492.3 m/sec. The 
method is sufficiently accurate to determine the variation in the speed of sound 
with temperature or salinity or to measure the component of tidal flow paralle! 
to the sound path. The practical limitations imposed on distance determina- 
tions by the use of sound signals in sea water are thus determinable. 


HE values of the speed of sound in sea water to be reported here are 

believed to be unique as to accuracy as well as interesting from other 
points of view. They were obtained incidentally during a program of 
work having other objectives. This fact will explain the limited number of 
the observations, all of which were made on a single afternoon. 

The method was essentially that described by Stephenson! except for 
modifications in the distance determinations which were made necessary 
by the much greater distances involved. A sound and a radio signal were 
transmitted simultaneously at a point approximately 53 nautical miles 
(61 miles) from the receiving apparatus. The latter served to measure 
the time interval between the arrival of the radio and sound signals. 
This time interval was of the order of 60 sec. and was measurable to the 
nearest .001 sec. The distance between the sound transmitter and the 
sound receivers was something over 321,000 ft. The end points of this 
distance were located on the map to within 6 ft. The resulting speed 
determinations should therefore bé accurate to better than one part in 
25,000. 

The reception and recording of the signals was accomplished by the 
Subaqueous Sound Ranging Section of the Coast Artillery Corps at 
Fort H. G. Wright who also furnished the map coordinates of the hydro- 
phone positions. The sound and radio signals were transmitted from the 
steamer Guide of the U.S. Coast and Geodetic Survey. The sound signals 


* Published by permission of the Director of the U. S. Coast and Geodetic Survey, 
the Chief of Coast Artillery and the Director of the U. S. Bureau of Standards. 
' E. B. Stephenson, Phys. Rev. 21, 181 (1923) 
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were made by the explosion of bombs which were suspended at a depth of 
42 feet vertically below the center of a pyramidal target of a type com- 
monly used by the Sound Ranging Section in similar work. The bombs 
contained about 314 pounds of T.N.T. and were fired electrically from 
the 110 volt supply circuit of the ship. The operation of a single switch 
accomplished both the firing of the bomb and the transmission of the 
radio signal. It had been previously demonstrated by experiment that 
with the arrangement used the closing of the circuit, the detonation of 
the bomb, and the transmission of the radio signal are simultaneous to 
within .001 sec. 

The position of the target was obtained at the instant of firing. Observ- 
ing parties with theodolites were detached from the Guide and occupied 
stations adjacent to two lighthouses the positions of which were known 
from previous triangulations. They sighted on the bomb target and saw 
it tilt when the bomb exploded. The distances “lighthouse to ship” were 
of the order of five miles. As the lighthouses were occupied eccentrically 
the observations were reduced to center. The resultant positions of target 
are comparable in accuracy with that of intersection stations in tertiary 
triangulation. The distance from each position of the target to each 
hydrophone was computed by geodetic methods. 

The Guide was anchored in about 100 feet of water and the water 
intervening between it and the receiving hydrophones ranged from 90 to 
120 feet in depth. 

The data follow. 

TABLE I 


Measurements of speed of sound in sea water, made Nov. 1, 1923. 
Temperature of water at position of boat at depth 42 ft =13°C 











Bomb Time Hydrophone Distance Time Speed Speed 
#1 1.55p.m. H; 321,900 ft 65.659 sec. 4902.6 ft/sec. 1494.3m/sec. 
H, 323,973 66.090 4902.0 1494.1 
#2 3.57p.m. H, 294,483 60.146 4896.12 1492.34 
Hz 297,879 60.840 4896.12 1492.33 
H; 321,741 65.714 4896.06 1492.32 
H, 323,808 66.135 4896.18 1492.36 
#3 4.06p.m. H, 294,498 60.148 4896.21 1492.37 
He 297,894 60.843 4896.12 1492.34 
H; 321,756 65.719 4896.00 1492.30 
Hg 323,826 66.138 4896.21 1492.37 
Average velocity for Bomb 41 = 4902.3 ft/sec. (1494.2 m/sec.) 


Average velocity for Bombs #2 and 3 = 4896.12 “ (1492.34 * ) 


It is quite obvious that the data for distance and time given in the 
above tables permit the computation of the sound speed between the 
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hydrophones, leaving out of consideration entirely the radio signal accom- 
panying the detonation of the bomb. In this way the speed over approxi- 
mately the last 30,000 ft of the 300,000 ft path is obtained. The results 
are assembled in the following table. 


TABLE II 
Inter-hydrophone sound speeds 








Bomb Hydrophones Distance Time Speed Speed 

#2 Hi—Hs 27,258 ft 5.568 sec 4895.4 ft/sec. 1492.1 m/sec. 
H,—H, 29,328 5.989 4893.9 1491.7 
H.-H; 23,862 4.874 4895.4 1492.1 
H.—H, 25,929 5.295 4897.2 1492.7 

43 H:—H; 27,258 5.571 4893.0 1491.4 
Hi—H, 29,328 5.990 4896.0 1492.3 
H.—H; 23,862 4.876 4893.6 1491.6 
H.—H, 25,932 5.295 4897.2 1492.7 





Average speed from inter-hydrophone measurements 4895.1 ft/sec. or 1492.1 m/sec. 


While the speeds obtained from bombs No. 2 and No. 3 are in satis- 
factory agreement, that obtained from bomb No. 1 differs from them by 
approximately 6 ft/sec. This difference can be due only to a real change 
and not to experimental errors. By referring to a tide table it was found 
that at the time bomb No. 1 was fired there was a fairly strong tidal cur- 
rent running in the direction of sound propagation, while several hours 
later when bombs No. 2 and No. 3 were fired the tide was practically at 
a standstill. An interval of only nine minutes separated the firing of 
bombs No. 2 and 3. Tidal conditions for them were therefore substan- 
tially the same. The speed obtained for bomb No. 1 might therefore be 
taken as the speed of sound with respect to the medium plus the speed of 
flow of the water in the direction of propagation of the sound. The excess 
of 6 ft/sec. would be accounted for by a tidal current of 4.1 miles per 
hour, if it is assumed that the speeds are simply additive. 

The Division of Tides and Currents of the U. S. Coast and Geodetic 
Survey has determined from data at hand what the current speeds at 
the time might have been and finds for the time of bomb No. 1 an average 
current of 1.8 miles per hour, for bomb No.2 0.6 miles per hour and for 
bomb No. 3 0.36 miles per hour, all in the direction of sound propagation. 
While these figures are deduced from an average of actual observations, 
the deviations from the mean are all too small to make a current of 4.1 
miles per hour probable. A large part of the 6 ft per second discrepancy 
between the two speed values is undoubtedly attributable to the water 


currents, but further work to clear up the whole discrepancy seems 
desirable. 
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The revived interest in the speed of sound is largely due to recent 
activity in applying it to the measurement of distances. When sea water 
is the medium involved, the currents, as well as inhomogeneities due to 
temperature and other causes limit the accuracy attainable. If speeds 
which are independent of the effects of flow of the water are wanted this 
may be accomplished by having duplicate equipment at each station and 
observing the time of travel for the sound in each direction. By sending 
the signals in both directions in rapid succession the effects of motion 
of the water may be completely eliminated. 


DIVISION OF HYDROGRAPHY, SUBAQUEOUS SOUND RANGING SEC1 
U.S. Coast AND GEODETIC SURVEY COAST ARTILLERY COoRPs, 
Fort H. G. WRIGHT 
SouND LABORATORY 
BUREAU OF STANDARDS 
June 25, 1924 
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ON THE PARTIALS OF A PIANOFORTE STRING STRUCK 
BY AN ELASTIC HAMMER 


By R. N. GHosH 


ABSTRACT 


Dynamics of a stretched string struck by an elastic hammer.— Neglecting 
small terms in the partial differential equation of motion, a solution is obtained 
involving the mass, elasticity and impact speed of the hammer and the distance 
of impact a from the nearer fixed end. Also approximate expressions are ob- 
tained for the pressure exerted on the string and for the duration of impact. 
The final solution for the displacement contains besides the Helmholtz term, 
three others. For a given a, the Helmholtz term disappears for partials with 
periods 2, 2/3, 2/5, etc., but not the other terms; this explains the observed 
existence of these partials with feeble intensity. The equation shows how the 
relative strengths of the partials depend on the various factors. If the duration 
of impact is half the period of vibration of a component, that component will 
become very prominent. For a given velocity of impact, the resulting am- 
plitude increases with a. 


INTRODUCTION 


ELMHOLTZ built up his classical theory! supposing that the dura- 
tion of impact was small compared to the period of vibration of the 
string, that it was not displaced appreciably from its equilibrium position, 


and that the hammer was elastic. Consequently he found that the dura- 
tion of impact was independent of the striking point. Kaufmann in the 
year 1895 showed however that the duration of impact was not indepen- 
dent of the striking point, and gave an approximate theory for the calcu- 
lation of the duration of impact when the striking distance was small, on 
the supposition (1) that the string was appreciably displaced from its 
equilibrium position, and (2) that the hammer was perfectly hard. The 
theory has since been extended by Raman and Banerji,’ and also by 
P. Das* for the calculation of the duration of impact at any point of the 
string. 

Recently it has been shown’ that the observed values of the duration 
of impact do not agree with those calculated from Kaufmann’s theory, 
at least in the middle octaves of the piano where the hammer is elastic. 
In the higher octaves where the hammer is almo-t hard Kaufmann’s 
theory is partly true. Further it has been shown that his theory can be 
easily extended to the case of an elastic hammer. 

' Helmholtz, Sensation of Tone (Ellis trans.) page 380 

* Kaufmann, Ann. der Physik, 54, 675 (1895) 

* Raman and Banerji, Proc. Roy Soc. 97, 99 


‘ P. Das, Proc. Indian Asso. Calcutta 7, 13 
* Ghosh, Indian Science Congress Jan. 1922 
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In the present paper an attempt is made to point out the differences 
in the relative strength of the partials in the cases stated above. The 
main results are, (1) that the partials whose periods are 2/3, 2/5, 2/7, 
etc. of the duration of impact do not completely disappear; (2) that the 
strengths depend upon the striking point; and (3) that certain undesirable 
partials may become prominent under some circumstances—depending 
upon the striking point, the elasticity and the mass of the hammer. 


Law OF MOTION OF THE STRING 


Let: /=length of the string; 
a =the striking distance from the nearer fixed end; 
M =the effective mass of hammer; 
=the strength of the elastic spring of the hammer felt; 
c=/T/p =velocity of transverse wave in the string; 
y=the displacement of the string; 
yo= the displacement of the center of mass of the hammer: 
=the compression of the hammer felt. 

The equation of motion of the hammer is given by 


Md’ y)/dt?= —p 
where p represents the pressure upon the hammer. 


youyté, and p=pé 


1 d da? 
and r\2 + zI- p=? "ss (approx. ). 


a * dt dt? 

Kq. (3) is derived on the assumption that the small part a of the string 
included between the striking point and the fixed end, follows the motion 
of the hammer, i.e. the wave takes very little time to reach the end from 


the striking point. Eliminating p from (1) (2) and (3) we obtain 


Mpés d? T ad d* MT 
«el elas eee {a+ — pis 
m pe dt at* pa 


Tdy T 
tities. w 
c dt a 
A complete solution of (4) is difficult and even if it were obtained, it 
would be too complicated to be of any practical value; hence we shall 
neglect the first two terms which are very small compared to the others; 
then (4) reduces to 
T \d@? T d T 
u( / : 


14+-— J—+ — —+— = 0. 5) 
er ca { 
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Since 7°/c is small, the solution of (5) is obtained in the form 








v~ Ae! sin qt (6) 
7 pa 
bua (— ) (7) 
Mc\T+ua 
T ma 
aM \T+ua 
If v is the velocity of impact of the hammer at ¢=0, then A is given by 
v( Tf1 k\ )7 
ati ———} } . (9) 
q\ pra 2cF | 


When the hammer is perfectly hard, » becomes infinite; then k becomes 
independent of the striking distance and g? comes out 


gi2=T/aM —3(T/Mo)?. (10) 


The duration of impact is approximatly equal to 7/g. It is observed from 
(8) and (10) that the duration of impact in the case of an elastic hammer is 
greater than with a hard hammer. In both cases however it increases with 
increasing distance of the striking point from an end. 

The amplitude of the resulting motion of the string is given by Eq. (9). 
It is proportional to the impinging velocity of the hammer and increases 
as a increases. 


THE FORCE ON THE STRING 
From (3) and (6) p is easily found to be 


p=ATe*'|singt/a+(R/c)cos(qt+¢)] (11) 
where R= \/q?+k?/4 and tan ¢=k/2g. 


From (11) it is observed that the pressure on the string is not zero at the 
beginning of the impact but is equal to TAg/c, and at the end of the 
interval r/q its magnitude is again the same; it becomes zero after a short 
interval 6¢ obtained by putting t= (a/q+6t) in Eq. (11). The fact that 
the pressure is not zero at the beginning of the impact is also confirmed 
by Raman’s curves® where the hammer is considered to be perfectly 
hard. The interval 6¢ however is very small; it is approximately given 
by a/c. Hence we shall for simplicity consider the interval +/q to be 
the duration of impact. This amounts to neglecting the second term in 
Eq. (11). Then # is given by . 

p=(AT, a)e*'sin gt. (12) 


®* Raman and Banerji, loc. cit.*, Figs. 1 and 2 
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Physically the neglect of the term (R/c)cos(gt+@) amounts to the 
assumption that the disturbance has traveled a long distance along the 
string but has not yet been reflected from the farther end, i.e. the duration 
of impact is smaller than the time taken by the wave to reach the farther 
end of the string; so that if we continually increase a, our approximations 
will be sufficient as long as the above condition is satisfied. 


CALCULATION OF THE PARTIALS 


The equation of motion of a point on the string in terms of normal 


d*9, src\* 2 
—— + [ —- ¢,= —?q, (1.3) 
dt’ l Ip 


where ® is the generalised component of force, and 


co-ordinates is given by’ 


®, = Asin(sra/l)e* "sin gt’ . 


If @, and d¢,/dt are both zero initially then for {> q 


w/q 


24 _ msa , . 
¢,=—- sin —— | e7}*” sin n(t—?’) sin gt'dt 
nlp l . 
0 
2A . smaf2q cos (ng/2z) sin n(t— 2/2) 
peed epee: fc. be Berni sete len baw. (14) 
nlp gn 


bk sin (uw; 2g) cosn(t+n/2q) — kw sin n(t+7/q) 


(q—n)2 4q(q—n) 





kr sin n(t—72/q) , 
le ——-|(approx.) (15) 
4q(q+n) 





where n=src, I. 
The final solution for the displacement y becomes, when we substitute 


the value of , to the same order of approximation, 


2Ar’ i sin (sra/l) sin (smc/l) { 2g cos (sr*c/2lq) sin (src/l)(t—/2q) 
y= —_ -_ d — omer 





prc s(lq—smc) | = (lg+src) 
>k sin (sr*c/2lq) cos { sxc, I(t+2/2q) } rk 


+ —— sin “(+4/q) 
— —_—_— 1 —— 
(iq—sxc) al lacie 





kr sin (src/l) (t—x/q) 


_ (16) 
4gl(lq+smxc)/(lq—smc) } 





* Rayleigh, Theory of Sound, Vol 1, Art. 128 
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Equation (15) shows besides the Helmholtz term three other terms 
which appear on account of the finite displacement of the string, but on 
account of the presence of k their magnitudes are small. The convergence 
of the first two terms is of the order s~*, but the last two terms converge 
slowly. They increase if k increases, that is when the striking distance 
increases. This partly explains the choice of the striking point which is 
about 1/7th the length of the string. Whenever the striking distance is 
large, the higher partials become prominent and disturb the musical 
quality of the note. 

For a given striking distance the first term disappears for those partials 
whose periods are 2/1, 2/3, 2/5, 2/7 etc. of the duration of impac:, but 
the other three terms do not, hence their existence with feeble intexsity 
is explained. When sin (s7a/l) =0 (i.e. the point struck is the first node 
of the sth component) then a=//s, and the component will be absent; 
but if at the same time the duration of impact 2/qg=I/cs, then 
sin (sr’c/2lg=1 and cos (sx*c/2lq) =0, and the first term and the last two 
terms disappear, but the second term does not and its magnitude is then 
approximately given by k/2s which is small and decreases with the order 
of the component. So we find that a component having a node at the 
point struck is not entirely absent if the duration: of impact is half the 
period of vibration of the component. Hipkins and others have found 
this to be the case experimentally.*® 

A given component will converge more rapidly the greater the value 
of g, which increases as a decreases (see Eq. 8), hence a must be small so 
that the higher undesirable partial may be very feeble. It is further 
observed that the magnitude of the second term is inversely proportional 
to (lg—smc)* hence when g=szc/l, that is when the duration of impact 
is equal to half the period of vibration of that component, it will become 
very prominent. 

A is not constant throughout all the length of the string for the same 
velocity of impact of the hammer. It changes with g and therefore in- 
creases as a the striking distance increases. 

Puysics DEPARTMENT, 

ALLAHABAD UNIVERSITY, 


ALLAHABAD, INDIA. 
December 7, 1923. 


§ Helmholtz, loc. cit.! pp. 76, 77 
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BOOK REVIEWS 


Color and Methods of Color Reproduction. By MArtTIN anD GAMBLE.—This book 
contains a great deal of interesting and useful information on color. The exposition of 
the more elementary portions is clear, and it is evident from the frequent references to 
easily observable phenomena in nature dependent on color, and to works of art, that 
the writer uses and enjoys his eyes. The chapters on colored light and coloring materials 
provide useful, practical information on methods of producing color, both in the labora- 
tory and in the arts. Detailed information on the chemical characteristics and the 
permanence of various coloring materials is something that the physicist working with 
color often wants and usually has difficulty in finding. 

Chapters on color vision, color blindness and color in human experience give an 
interesting general review of the physiological and psychological aspects of color, treated 
in a fairly popular manner. Photometry and color measurement are considered at some 
length, but not entirely satisfactorily, since a number of erroneous statements, partic- - 
ularly of experimental findings in the field of intermittent vision, indicate that the 
author’s acquaintance with this part of the subject has been gathered from rather hasty 
and undigested reading. The method given of calculating hue, luminosity and purity 
from spectrophotometric data does not show acquaintance with the more recent work 
on color measurement transformations. In particular the “saturation” arrived at is not 
the saturation with respect to the spectrum, but to a non-existent “two sensation” 
spectrum. 

Two chapters of a very popular and non-technical character by the second author, 
on ‘Color printing” and ‘‘Color photography” would have been of greater assistance 
to a clear presentation of three-color printing processes if placed in inverted order. 
The scientific basis of these processes and of three-color photography are not clearly 
grasped by the authors. “Red, yellow and blue” are given as the primary printing 
colors,and the author of the first part of the book, while in one place correctlystating that 
the ideal subtractive primaries would be “peacock blue,” ‘‘chrome-yellow” and ‘‘cerise,” 
does not stand by this consistently. His colored plates illustrative of additive and sub- 
tractive color mixture (p. 17) show in each case the complementary of red as blue. The 
foot-note to page 166 which endeavors to make clear the three color additive process of 
color photography by optical superposition of the three color images is unfortunate in 
the parenthesis ‘“‘both for taking and viewing.”” The simple principle that the trans- 
mission of the “taking” screens should copy the color mixture curves of the spectrum 
while the viewing screens should be the monochromatic primaries used for determining 
the color mixture curves was clearly enunciated at least thirty years ago. 

These illustrations are sufficient to show that, interesting and valuable as this book 
is in many ways, it nevertheless cannot satisfy the student who wants sound detailed 
information. It is only fair to state however that a comparatively small number of text 
changes would suffice to correct the errors which have been noted. Such changes would 
make a second edition a much more valuable book.—Pp. 187, 73 fig., Van Nostrand, 
1923. HERBERT E, IvVEs 


The Elementary Principles of Lighting and Photometry. By J. W. T. WaLtsn.—To 
quote the author, this book “‘is mainly intended to give a description of the nature and 
amount of the illumination required for different purposes, and of the way in which the 
desired result may be attained and its attainment checked by photometric measure- 
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ment.” Thus while a large part of the subject matter is divided into ‘‘Indoor Lighting,”’ 
“Industrial and School Lighting,”’ ‘Outdoor Lighting,’’ ‘Light Projection’’ and such 
practical illumination matters, these are well introduced by several chapters dealing 
with “Light, Vision and the Eye,” and an outline of photometry. 

The chapters dealing with photometry present a well balanced review of the elemen- 
tary aspects of the subject. Good judgment is shown in the relative space given to 
subjects of academic rather than practical interest, such for instance as that of funda- 
mental photometric standards. The treatment of color and color difference photometry 
is not so satisfying. The rating of illuminants by average candle power rather than 
lumens follows an English practice which is some years behind our own. Omission of 
any discussion of luminous efficiencies or any adequate treatment of the physics of light 
production is to be regretted. A carefully chosen bibliography at the end, together with 
appendices giving the most authoritative definitions of photometric quantities will assist 
the reader to fill in these gaps. 

The guidance offered to the practical problems of lighting is good. ‘‘The important 
factor in seeing is the brightness of the thing looked at and this is the product of the 
illumination and the reflection ratio,” but the qualifications to this introduced by the 
effects of bright lights in the field of view, and glare in general, are well covered. Specific 
treatments of lighting schemes for individual rooms in the house should prove of real 
help to people of the class to whom this book is addressed, namely those possessed of an 
elementary scientific and technical education, who are, let us say, not lighting engineers 
but nevertheless are desirous of lighting their living rooms safely and sanely.— 

HERBERT E. Ives 


Light and Color in Advertising and Merchandising. By M. LuckKe1sH.—This book 
does not purport to be a scientific treatise on light and color. To quote the author “‘it 
seems unwise to go into these tedious details here.” The chief interest of the book to 
the physicist will be the revelation it gives of how far reaching may be the reactions, on 
the human complex, of certain rather simple physical variables. Thus, starting with 
radiant energy which can be varied physically merely in intensity, frequency and direc- 
tion of propagation, our visual apparatus attributes it to this or that localized material 
object, and adds the quality of color. Our mental experiences differentiate between 
pleasant and unpleasant arrangements of the visual field of view, and sets of color associa- 
tions arise, harmonious and unharmonious. These become in time part of our language 
and literature. Ultimately a study of how these factors may be intelligently applied by 
the advertising specialist, the interior decorator, and the window trimmer, comes to be 
written, showing how these practical men who know nothing of lumens or Angstroms, 
can make real money by proper use of light and color. 

The book consists essentially of a popular review of what is known of the psychology 
of light and color, with numerous specific commercial applications. A large number of 
colored illustrations make the book a very attractive one, even to the reader whose 
interest in the subject matter may be slight.—Pp. xvi+268, 38 illustrations, Van Nos- 
trand, 1923. HERBERT E. Ives 








